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We  recently  showed  that  the  spectrum  of  light  emitted  by  a  source  depends  not  only  on  the  spectrum  of  the  source  distribution 
but  also  on  the  degree  of  spectral  coherence  of  the  source  fluctuations.  In  this  note  we  show  that  with  a  degree  of  spectral  coherence 
of  certain  kind,  specified  by  two  parameters,  the  spectrum  of  the  emitted  light  will  be  displaced  relative  to  the  source  spectrum. 
The  displacement  will  be  either  toward  the  lower  or  toward  the  higher  frequencies,  depending  on  the  choice  of  the  parameters. 


1.  Introduction 

It  has  been  known  for  some  time  that  the  spectrum 
of  light  generally  changes  on  propagation,  even  in  free 
space  [1,2],  Such  changes  are  basically  due  to  corre¬ 
lation  properties  of  the  source.  Recently  we  derived 
a  condition  for  the  normalized  spectrum  of  light  gen¬ 
erated  by  a  planar,  secondary,  quasi-homogeneous 
source  to  be  the  same  throughout  the  far  zone  and  in 
the  source  plane  [3]  .  We  referred  to  this  condition, 
which  is  a  requirement  on  the  functional  form  of  its 
degree  of  spectral  coherence,  as  the  scaling  law  and 
we  noted  that  all  quasi-homogeneous  lambertian 
sources  satisfy  this  law.  We  have  also  shown  that  when 
the  scaling  law  is  not  satisfied  the  spectrum  of  the 
emitted  light  will,  in  general,  no  longer  be  invariant 
on  propagation.  These  theoretical  predictions  have 
been  recently  verified  by  experiments  [4  j. 

In  another  recent  paper  (5]  we  considered  radia¬ 
tion  from  three-dimensional,  quasi-homogeneous 
sources  and  we  showed  that  if  the  source  spectrum 
consists  of  a  line  with  a  gaussian  profile  and  if  the 
degree  of  spectral  coherence  of  the  source  is  appro¬ 
priately  chosen,  the  spectrum  of  the  emitted  light  will 
also  consist  of  a  line  with  gaussian  profile,  but  this 
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line  will  be  redshifted  with  respect  to  the  spectral  line 
of  the  source  distribution.  The  amount  of  the  red- 
shift  depends  on  the  spectral  correlation  length  of  the 
source.  This  result  has  important  implications  for 
astrophysics,  some  of  which  were  briefly  mentioned 
in  ref.  [  5  ] . 

In  the  present  note  we  again  consider  a  source 
whose  spectrum  consists  of  a  single  line  with  a  gaus¬ 
sian  profile  but  we  assume  somewhat  different  cor¬ 
relation  properties  of  the  source.  More  specifically 
we  choose  a  degree  of  spectral  coherence  of  the  source 
distribution  which  depends  on  two  parameters  rather 
than  on  a  single  parameter  as  we  have  done  previ¬ 
ously.  The  spectrum  of  the  emitted  light  is  again 
found  to  be  a  line  with  gaussian  profile,  but  this  line 
may  be  redshifted  or  blueshifted  relative  to  the  spec¬ 
tral  line  of  the  source  distribution,  depending  on  the 
choice  of  the  parameters. 

2.  The  spectrum  of  light  produced  by  a  three- 
dimensional  quasi-homogeneous  source 

Let  us  consider  a  fluctuating  source-distribution 
Q(r,  t)  occupying  a  finite  domain  of  volume  D  in 
free  space  and  let  T(r,  t)  denote  the  field  generated 
by  the  source.  Here  r  denotes  the  position  vector  of  a 
typical  point  and  l  the  time.  Both  l'(r.t)  and  Q{  r.  I) 
are  taken  to  be  analytic  signals  [6],  They  are  related 
by  the  inhomogeneous  wave  equation 


12  0  030-4018/87/$03.50  <fc)  Elsevier  Science  Publishers  B.V. 

The  U.S.  Government  Is  authorlaed  to  reproduce  and  sell  this  report.  ( North-Holland  Physics  Publishing  Division) 

Permission  for  further  reproduction  by  others  must  be  obtained  from 
the  copyright  owner. 


Volume  62,  number  I 


OPTICS  COMMUNICATIONS 


I  April  1987 


F2V(r,t)-c~2(d2/dt2)  V(r,t)  =  -4n  Q(r,t)  . 

(2.1) 

We  will  assume  that  the  statistical  ensembles  that 
characterize  the  source  fluctuations  are  stationary. 
Let  fVa(rh  r2,  co)  and  W/,  (*V  r2,  co)  be  the  cross- 
spectral  densities  of  the  source  distribution  and  of 
the  field  distribution  respectively.  They  may  be  rep¬ 
resented  in  the  form  [7] 

WQ(r,,r2,ai)  =  (UZ(r,,a>)  UQ(r2,  co))  ,  (2.2a) 

W'.  (r, ,  r2,  co)  =  {  ,  co)  Uy(r2,  co)  >  ,  (2.2b) 

where  {V^ r,  co))  and  (LV(r,  o>)}  are  ensembles  of 
suitably  chosen  realizations,  angular  brackets  denote 
averages  taken  over  these  ensembles  and  the  asterisk 
denotes  the  complex  conjugate.  As  consequence  of 
the  wave  equation  (2.1)  the  two  cross-spectral  dens¬ 
ities  may  be  shown  to  be  related  by  the  equation  n 
(ref.  (7],  eq.  (3.10);  ref.  (8a],  eq.  (2.1 1 )] 

(F? +*2)  (F\+k2)  W>{rur2,co) 

=  (4n)2  WQ(r„r2>co)  ,  (2.3) 

where  Ff  and  F \  arc  the  laplacian  operators  acting 
with  respect  to  the  coordinates  of  the  points  r ,  and  r2 
respectively  and 

k  =  co/c  (2.4) 

is  the  wave  number  associated  with  the  frequency  co, 
c  being  the  speed  of  light  in  vacuo. 

Using  eq.  (2.3)  one  can  show  that  the  radiant 
intensity  J,„(u)  generated  by  the  source,  i.e.  the  rate 
at  which  energy  is  radiated  at  frequency  co  per  unit 
solid  angle  around  a  direction  specified  by  a  unit 
vector  « is  given  by  ( ref.  ( 8a  ] ,  eq.  ( 3.9 ) ] 

•/„(«)  =  (2ar)6  tf\j{-ku,  ku,  co)  ,  (2.5) 

where 

froUC, .  K2,  co)  =  (2*) J  jVtf(r, ,  r2,  co) 

/>D 

xexp(  —  i ( #Ti  t,  -t-iTiT;)]  d’r,  d3r2  (2.6) 

ri  The  derinition  of  the  cross-spcctral  densities  employed  in  refs. 
[7|  and  [8]  differ  by  complex  conjugation.  Throughout  this 
note  we  employ  those  of  ref.  ( 7  ] ;  hence  some  of  the  formulas 
we  now  use  |e.g.  eq.  ( 2.5)  below]  differ  trivially  from  the  cor¬ 
responding  formulas  of  refs.  ( 8 1 . 


is  the  six-dimensional  Fourier  transform  of  ll ’0. 

We  will  restrict  our  attention  to  quasi-homogene- 
ous  sources.  For  such  sources  one  has,  to  a  good 
approximation, 

^0(r,,r2,  co) 

=  S0[(r,  +  r2)/2,  cu]  H(j(r2  -r,.co)  ,  (2.7) 

where 

SQ(r,co)  =  WQ(r,  r,  co) 

=  <UJ(r,w)  U0(r,co))  (2.8) 

is  the  source  spectrum  and 

Ve(r2-rl,io)s  WQ(rly  r2,  co) 

X  (S^fr, ,  (u)]  ~ 1/2  (S0(r2, <u )] " 1/2  (2.9) 

is  the  degree  of  spectra!  coherence  of  the  source  dis¬ 
tribution.  Moreover,  for  each  effective  frequency  co 
contained  in  the  source  spectrum,  Sy(r,  co)  varies 
much  more  slowly  with  r  than  co)  varies  with 
r.  With  sources  of  this  class  eq.  (2.5)  takes  the  form 
(ref.  (8b], eq.  (3.11)] 

J<0(u)  =  (2ji)6Sv(0,co)fiQ(ku,co)  ,  (2.10) 

where  the  tilde  now  denotes  three-dimensional  Four¬ 
ier  transforms. 

Let  us  next  assume  that  the  source  spectrum  is  the 
same  at  each  source  point.  We  will  then  write  Sy(o>) 
in  place  of  S0(r,co).  In  this  case  £y(0,  co) 
=  DSq(co)/( 27t)Jand  the  formula  (2.10)  becomes  ( D 
again  denoting  the  source  volume ) 

J,0(u)  =  {2n)2DSQ(co)fi0(ku,co)  .  (2.11) 

Now  the  radiant  intensity  J„,(u)  is  trivially  related 
to  the  spectrum  S^'iRu,  co)  =  W\^'(Ru,  Ru,  co)  of 
the  far  field  by  the  formula  [9] 
SI"  '(Ru,  Ru,  co)  ~J,„(u)/R2  as  kR-*cc ,  with  the 
unit  vector  u  fixed.  Hence  we  obtain  at  once  from  eq. 
(2.11)  the  following  expression  for  the  spectrum  of 
the  emitted  light  in  the  far  zone: 

Sr  ’(«.«)  =  ( 2n)'{D/R 2 )Sv{(o)fi(J{kii,co)  . 

(2.12) 

This  formula  shows  that  the  spectrum  Si"  ’(u.co)  of 
the  emitted  light  in  the  far  zone  depends,  in  general, 
not  only  on  the  source  spectrum  Sy(w)  but  also  on 


i  I  .  «  ,  i’  »i.  it ,  ii  •.!, 
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the  degree  of  spectral  coherence  of  the  source  distri¬ 
bution.  It  seems  worthwhile  to  note  that  the  dimen¬ 
sions  of  JC>  and  of  Sy  are  different.  Since  ftQ  is  the 
three-dimensional  Fourier  transform  of //y,  [/i^J  =  /.  •* 
(brackets  denoting  dimensions  and  L  denotes 
length).  Hence  eq.  (2.12)  implies  that 
[S','  ’]  =  [Sy]L4,  in  agreement  with  eq.  (2.3). 


3.  A  class  of  source  correlations  that  generate 
lineshifts 


In  ref.  [5]  we  considered  quasi-homogeneous 
sources  whose  spectrum  was  a  line  of  gaussian  profile. 


Sy(w)=Aexp[-(w-a»0)2/2<5,2)]  ,  (<5()/w„«l ), 


and  whose  degree  of  spectral  coherence  was  also 
gaussian  viz.. 


Hu(r ,  «)=exp[  -  r'2/2<x2((o)]  , 


where  r'  =  |r'|.  The  three-dimensional  Fourier 
transform  of  nv  is  then  given  by 


fi(j( K.  aj)  =  [a(to)/K  27t]3exp[  -  }A2<r2(w)], 


(A  =  |A| ).  In  particular  we  showed  that  if  <x(co)  is 
constant  ((  say)  such  a  source  will  emit  light  whose 
spectrum  in  the  far  zone  is  rcdshifted  with  respect  to 
the  source  spectrum,  the  amount  of  the  shift  depend¬ 
ing  on  the  effective  source  correlation  length  £. 

The  degrees  of  spectral  coherence  of  the  form  (3.2), 
with  <?(«)=£  (constant)  form  a  one-parameter 
family.  In  this  note  we  will  consider  quasi-homoge¬ 
neous  sources  whose  degrees  of  coherence  are  of  a 
somewhat  more  general  form.  Specifically  we  assume 
that  for  these  sources 


fiQ(K)  =  Bcxpl  -  {( A  -  A,  )2£2]  , 


where  B,  A',  and  C  ate  positive  constants.  We  have 
written  fiv(K)  rather  than  A.  to)  on  the  left-hand 
side  of  eq.  (3.4),  because  ftv  is  now  independent  of 
to.  Only  two  of  the  three  constants  in  the  expression 
(3.4)  are  independent,  because  the  Fourier  trans¬ 
form  /iy(r')  of  /iy(  A)  satisfies  the  requirement  that 
/iv(0)  =  |,  which  is  a  necessary  condition  for  fiyir  ) 
to  be  a  correlation  coefficient. 


(units  :  10  cm  ) 


Fig.  1.  The  behaviour  of  the  correlation  coefficient  py(r)  =  ((sin 
A,r' )/K,r  ]  exp(  -r'1/ 2(!).  with  K,  =  1 .07*  10'  cm  ',  Z  =  1.5  cm 
(associated  with  curve  (d >  in  fig.  2). 


It  can  be  shown  by  a  long  but  straightforward  cal¬ 
culation  (which  we  omit  because  of  limitation  of 
space)  that  if 


*iC»l 


the  degree  of  spectral  coherence,  whose  Fourier 
transform  is  given  by  eq.  ( 3.4),  is 


Mv(r')  =  l( sinAT, r' )/AT, r' ]  exp( -r'2/2C2)  (3.6) 


and  that  the  constant  B  in  eq.  ( 3.4)  is  given  in  terms 
of  the  two  other  parameters  by  the  formula 


B=C/2(2n)V2At  . 


From  now  on  we  will  only  consider  situations  for 
which  the  constraint  (3.5)  holds.  Eq.  (3.6)  then 
shows  that  the  degree  of  spectral  coherence  has  the 
form  of  the  sine  function  (sin  A, r')/A,r', modulated 
by  the  gaussian  function  exp(  -r'2/ 2£;).  The  behav¬ 
iour  of  such  a  two-parameter  correlation  coefficient 
is  shown  in  fig.  1 . 

It  follows  on  substituting  from  eqs.  (3. 1 )  and  (3.4) 
into  eq.  (2.12)  that  the  spectrum  of  the  light  in  the 
far  zone,  generated  by  such  a  source,  is  given  by 


,SV  ’(«)  =  ( 2it )'(/)//? 2 ).4Bexp|  -  ( to -<o0):/26;,] 


Xexp{  -  (r u-<Oi 


1-wUUimilW 
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where  to = kc  as  before  and 

w,=A'|C\  <$,  =c/y  .  (3.9) 

We  have  written  5}*' (w)  rather  than  w) 

on  the  left-hand  side  of  eq.  (3.8)  since,  because  of 
the  assumed  isotropy  of  the  source,  S\^  '  is  now 
independent  of  u.  In  terms  of  the  parameters  w  ,  and 
<5,  the  factor  B,  given  by  eq.  (3.7),  becomes 

B=ci/2(2it)iaoj}dl  .  (3.10) 

Let  us  now  consider  the  expression  (3.8)  more 
closely.  For  this  purpose  it  is  convenient  to  set 

au  =  1/2  Jo,  a,  =  l/2<*i.  (3.11) 

One  then  finds  after  a  straightforward  calculation  that 
eq.  ( 3.8)  may  be  expressed  in  the  form 

.S',3'  l(w)=/lCexp[  -(w-w0,)2/2<5s,]  ,  (3.12) 

where 

<o,„  =(anco„  +  alojl)/(a0  +  al)  ,  (3.13) 

I/*,’,,  =2(a0 +  «, )  =  ( l/*o)  +  ( 1/*? ) ,  (3.14) 

and 

C=(2n)'(DB/R2) 

X  exp{  -  [a0O!|/(ao  +  «,  )](w,  -w0);)  • 

(3.15) 

The  formula  (3.12)  shows  that  the  spectrum  of  the 
emitted  light  in  the  far  zone  is  also  a  line  with  gaus- 
sian  profile,  but  it  is  not  centered  on  the  frequency 
a»<>  of  the  source  spectrum  [ cf.  eq.  (3.1)]  but  rather 
on  the  frequency  to,,,.  given  by  eq.  (3.13).  Since 
according  to  eqs.  (3.1 1 )  nr0  and  «,  are  positive  con¬ 
stants  one  can  readily  deduce  from  the  expression 
(3.13)  that 

Wo,  <w„  when  to,  <w„  , 
and  that 

Win  xo,i  when  to,  >  w„  . 

Since  according  to  eq.  (3.9)  w,~A',c,  this  result 
implies  that  if  the  parameter  K,  of  the  degree  of  spec¬ 
tral  coherence  ( 3.6)  is  smaller  than  the  wavenumber 


k„— to  0/c  associated  with  the  source  spectrum  SQ  (oj), 
the  spectrum  S,'1  (to)  of  the  emitted  light  is  red- 
shifted  with  respect  to  S0(to);  and  that  ifK,  is  greater 
than  k„  it  is  blueshifted  with  respect  to  it.  We  also  sec 
from  eq.  (3.14)  that  l/*<2„  >  l/*2,  i.e.  that  *0i<*u- 
Hence  in  either  case  the  spectral  line  of  the  emitted 
light  is  narrower  than  the  spectral  line  of  the  source 
distribution. 


4.  Examples 

To  illustrate  the  preceding  analysis  we  consider  a 
few  examples.  For  simplicity  we  will  choose 


*i=*n-  (4.1) 

Then,  according  to  eq.  (3.11),  £V|=a0  and  the 
expression  (3.12)  becomes 

S'i^’tw)  =.4f  exp[ - ( w  —  w)2/<5i;i]  ,  (4.2) 

where 

w=](w0  +  W|),  (4.3) 

C=  ( 2tt)3( DB/B) exp [  -  (w,  -wn)2/4*o]-  (4.4) 


We  see  that  the  spectral  line  of  the  emitted  light  is 
now  centered  on  the  average  value  w  of  the  frequen¬ 
cies  w(>  and  w,. 

Let  us  consider  the  normalized  spectrum 

_  x. 

s\*'((0)=SlS'(Qj) I  dw  (4.5) 

of  the  emitted  light.  On  substituting  from  eq.  (4.2) 
into  cq.  (4.5)  and  on  using  eq.  (4.3)  we  obtain  the 
following  expression  for  s\  r  ’  ( w): 

4*  ,(w)  =  (l/*„v  n) 

xexp-  -  [w-  J(w„  +w,  )]-/*r,}  .  (4.6) 

In  fig.  2  curves  arc  plotted  showing  the  normalized 
source  spectrum 

st,(w)  --  (l/(5„v  '2n  )  cxp|  -  (w-w,i)’/2*r,)  . 


(4.7) 
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( units’  I0'5  sec  ') 


Fig.  2.  Redshifts  and  blueshifts  of  spectral  lines  caused  by  source 
correlations.  The  normalized  spectrum  3y(<u)  of  the  source  dis¬ 
tribution  is  a  line  of  gaussian  profile  [given  by  eq.  (4.7)],  with 
<o„=3.l97049x  10"  s  1  (sodium  line  of  wavelength 
A„=5895.924  A)  and  rms  width  S  =2x10"’  s  Curves 
( a )  -  ( d )  show  the  normalized  spectra  „f  the  emitted  light  [  lines 
with  gaussian  profiles  given  by  eq.  (4.6)],  generated  by  the  source 
distribution,  each  with  <5,=<5„  (f  =  c/<5,  =  1.5  cm)  and  with 
m,  =<u(l-  1.4x  10"  s  '  (a),  0.7x10"  s  1  (b), 

m,=(u„  +  0.7xl0"s  '(c)  and <o,  =  tog  +  t.4x  10"  s ~ 1  (d). 

taken  to  be  one  of  the  sodium  lines,  as  well  as  a  num¬ 
ber  of  emitted  lines  for  different  values  of  the  param¬ 
eter  uj,=K,c ,  of  the  degree  of  spectral  coherence  of 
the  source;  the  other  parameter,  £,  is  kept  fixed  and 


chosen  so  that  S,=c/{  is  equal  to  60 .  It  is  seen  that 
with  increasing  values  of  the  difference  |  w0-cu,  |  the 
shift  of  the  emitted  spectral  line  also  increases.  This, 
of  course,  is  to  be  expected  since  when  <5,  =<50,  the 
shift  is  given  by  1 to  -  to,,  \  =  1 1  tu0  -  w ,  | . 
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A  general  expression  is  derived  for  the  ratio  of  the  radiated  power  and  the  source-integrated  intensity  for  any  planar  gaussian 
Schell-model  source.  The  behavior  of  this  quantity,  known  as  the  radiation  efficiency  of  the  source,  is  displayed  graphically  as  a 
function  of  the  rms  width  of  the  intensity  profile  and  the  spatial  coherence  length  of  the  light  distribution  across  the  source.  Some 
limiting  cases  are  discussed  and  il  is  shown  that  a  gaussian-correlaled  quasi-homogeneous  source  may  have  higher  radiation 
efficiency  than  a  fully  coherent  Schell-model  source  with  a  gaussian  intensity  profile  (e  g.  a  single  mode  laser). 


1.  Introduction 

In  the  last  few  years  there  has  been  considerable 
interest  in  radiation  produced  by  partially  coherent 
sources.  In  particular,  the  radiation  efficiency  of 
sources  of  different  states  of  coherence  have  been 
investigated.  The  radiation  efficiency  of  a  source  is 
defined  as  the  ratio  of  the  total  outgoing  flux  to  the 
source-integrated  intensity.  The  radiation  efficiency 
of  planar  quasi-homogeneous  sources  was  studied  by 
Carter  and  Wolf  [1,2].  More  recently  the  radiation 
efficiency  of  three-dimensional  gaussian  Schell-model 
sources  was  calculated,  and  was  compared  with  the 
radiation  efficiency  of  a  corresponding  coherent 
source  [  3  ] . 

In  the  present  paper  we  extend  the  analysis  of 
Carter  and  Wolf  to  the  important  class  of  planar 
gaussian  Schell-model  sources.  We  derive  an  explicit 
expression  for  the  radiation  efficiency  of  sources  of 
this  class  and  present  diagrams  which  show  its 
dependence  on  the  rms  widths  of  the  intensity  pro¬ 
file  and  its  degree  of  coherence.  We  also  examine 
homogeneous  sources,  completely  coherent  sources, 
and  quasi-homogeneous  sources  as  limiting  cases  of 
gaussian  Schell-model  sources.  Finally  we  derive 
conditions  under  which  a  planar  gaussian-correlated 
Schell-model  source  is  more  efficient  than  a  com- 

*  Research  supported  by  the  Air  Force  Geophysics  Laboratory 

under  AFOSR  Task  23IOGI  and  by  the  Army  Research  Office. 


Fig.  1.  Illustration  of  the  notation.  P  represents  a  field  point  in 
the  far  zone. 

pletely  coherent  source  with  a  gaussian  intensity  pro¬ 
file  (e.g.  a  single  mode  laser). 


2,  The  radiation  efficiency  of  a  planar  gaussian 
Schell-model  source 

Consider  a  planar  secondary  Schell-model  source 
occupying  a  domain  D  in  the  plane  z  =  0  and  radiat¬ 
ing  into  the  half-space  z>  0  (see  fig.  1 ).  Such  a  source 
is  characterized  by  a  cross  spectral  density  function 
of  the  form  [  4  ] 

fF(r, ,  r2,  to) 

=  [/(r,,  w)  /(r2,  to)]'/J£(r,  -r2,  to)  ,  (1 ) 

where  /(r,  w)  is  the  intensity  profile  and  g(r,-r2. 
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i o )  is  the  complex  degree  of  spatial  coherence,  both 
taken  at  frequency  co.  The  symbols  r,  and  r2  are  posi¬ 
tion  vectors  of  typical  points  in  the  source  region  D. 

It  is  known  that  the  radiant  intensity  produced  by 
a  secondary  planar  source  in  the  direction  specified 
by  a  unit  vector  s  is  given  by  [  5  ] 

J<u(s)  =  (2nk)2cos2d  W(ks  L  ,-ksx  ,  co)  .  (2) 

Here  k=co/c  ( c=  speed  of  light  in  vacuum)  is  the 
wave  number  associated  with  frequency  co,  6  is  the 
angle  between  the  s  direction  and  the  normal  to  the 
source  plane,  and 

•JL 

^(/■|./2.w)  =  (2k)-4  ||  H'(r,,r„e>) 


Here  /0,  a,,  and  ax  are  positive  quantities  depending 
only  on  the  frequency  co  (dependence  not  dis¬ 
played).  Such  sources  are  known  as  gaussian  Schell- 
model  sources. 

On  substituting  from  eqs.  ( 6 )  into  eq.  ( l )  and  tak¬ 
ing  the  four-dimensional  spatial  Fourier  transform 
one  can  show  after  a  lengthy  calculation  that  [cf.  refs. 
6  and  7] 

tV(ksx  .  —  ks,  .  co) 

/t) _ /  Ar2sin  -  67  \ 

4(27t)2a2(a2 +  2/?2)eXf\  2(a2  +  2p2))'  *  ^ 

where 

a2  =  l/(4o;)  .  p2  —  l/(2o2)  .  (8) 


xexp[-i(/  r,  +/2r>)]  d2r,  d2r2  (3) 

is  the  four-dimensional  spatial  Fourier  transform  of 
the  cross-spectral  density  of  the  light  distribution  in 
the  source  plane,  with  /,  and  f2  representing  two- 
dimensional  spatial-frequency  vectors. 

The  total  flux  emitted  by  the  source  into  the  half¬ 
space  z>  0  is  given  by  the  expression 

<P,„=  |  J„As)dS2,  (4) 

(i/ti 

where  the  symbol  (2n)  under  the  integral  sign  indi¬ 
cates  that  the  integration  is  taken  over  the  solid  angle 
subtended  by  a  hemisphere  in  the  half-space  z>0, 
centered  at  the  origin. 

We  define  the  radiation  efficiency  of  a  source  [cf. 
ref.  2.  eq.  (3. 1 1 )]  by  the  formula 


(5) 


The  integration  in  the  denominator  of  eq.  (5)  is  taken 
over  the  source  domain  D.  We  show  in  the  Appendix 
that  <  1  for  any  planar  source. 

We  will  now  consider  planar  Schell-model  sources 
for  which  both  the  intensity  distribution  and  the 
degree  of  spatial  coherence  are  gaussian,  i.e.  they  have 
the  form 

Hr.  co)  =  /()exp( -r2/2aj)  ,  (6a) 

and 


g(r,-r2.  co)  =exp[  -  (r,  -r2)2/2tx2]  . 


Next  if  we  substitute  !rom  eq.  (7)  into  eq.  (2)  we 
obtain  the  following  expression  for  the  radiant  inten¬ 
sity  generated  by  a  source  of  the  type  we  are 
considering: 


./,„(*)  = 


k2L 


4or(«2  +2  p:) 


cos 2  6 


xexp 


( 


k2sin2e  \ 
2(a2  +2P2) ) 


(9) 


It  follows  on  substituting  this  expression  into  eq.  (4), 
that  the  total  flux  at  frequency  co  radiated  by  a  planar 
gaussian  Schell-model  source  into  the  half-space  z>0 
is  given  by 


to  ~ 


k2I0 

4«:(ft 2  +  2  P2) 


|  cos 2  6 
(2*1 


:“p(~:7SW)(l~cos!6,)dfi  "0I 


After  some  algebraic  manipulation  this  expression 
can  be  reduced  to 


<ft,u=2;tff;/„|l— CXp(g  *•  )  |  exp(r)  dt| .  (11) 


where 

i2  =  [M2(ka,)2  +  2/{ko,)2\  (12) 


The  denominator  ineq.  (5)  with/(r,  co)  given  by  eq. 
(6a)  can  also  be  readily  evaluated  and  we  find  that 


(6b) 
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j"  I(r,  (o)  d2r=2rtl0a}  .  (13) 

On  substituting  from  eqs.  (11)  and  (13)  into  eq.  ( 5 ) 
we  finally  obtain  the  following  expression  for  the 
radiation  efliciency  of  a  planar  gaussian  Schell-model 
source: 

«»  =  1 -*>({)/«,  04) 

where 

c 

£>({)=exp(-£2)  J  exp(r2)df  (15) 

o 

is  the  Dawson  integral  [8], 

Fig.  2  shows  a  three-dimensional  plot  of  the  radia¬ 
tion  efficiency  tlo  as  a  function  of  ka„  and  ka,  calcu¬ 
lated  from  eqs.  (14)  and  (12).  Fig.  3(a)  shows  the 
behavior  of  the  radiation  efficiency  as  a  function  of 
kaK  and  fig.  3(b)  shows  its  behavior  as  a  function  of 
ka,  for  some  selected  values  of  the  other  parameter. 


3.  Physical  interpretation 

As  can  be  seen  from  eq.  (14)  the  radiation  effi¬ 
ciency  t„,  depends  on  the  rms  widths  of  the  intensity 
profile  and  of  the  degree  of  spatial  coherence  only 
through  the  parameter  defined  by  eq.  (12).  A  con¬ 
sequence  of  this  fact  is  an  equivalence  theorem  for 
the  radiation  efficiency:  there  exist  an  infinite  num¬ 
ber  of  planar  gaussian  Schell-model  sources  of  differ¬ 
ent  rms  intensity  width  0\  and  different  spectral 
coherence  lengths  at  which  have  the  same  radiation 


Fig.  3.  Radiation  efficiency  «„  as  a  function  of  ka,  for  selected 
values  of  ka ,  (a)  and  as  a  function  of  ka,  for  selected  values  of 
*<r,(b) 


XlTg 


Fig.  4.  Contours  of  equal  radiation  efficiency  as  a  function  of  Iter, 
and  ka,. 
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Fig.  5.  Graphical  representation  of  £)(£)/{,  where  D(4)  is  the 
Dawson  integral  (15). 


efficiency  t,u.  For  a  given  value  of  a  class  of  equiv¬ 
alent  sources  is  represented  by  a  single  curve  in  fig. 
4. 

We  shall  now  consider  a  number  of  limiting  cases 
that  are  of  special  interest. 


W(ru  r2,  w)=/0exp[  -(r,  —  r2) 2/2ey|]  .  (17) 


Since  W(r„  r2,  to)  now  depends  on  r,  and  r2  only 
through  the  difference  r,  -  r2,  the  source  is  homoge¬ 
neous.  It  follows  from  eq.  (12)  that  as 
t-tkOg/'Jl  and  hence  the  expression  (14)  for  the 
radiation  efficiency  now  becomes 


D(koglJl) 

kojjl 


(18) 


Furthermore,  on  inspecting  eq.  (18)  and  using  the 
fact  that  the  function  D(£)/£  monotonically  decreases 
from  unity  to  zero  as  £  increases  from  zero  to  infinity 
( see  fig.  5 )  we  see  that  for  radiation  from  a  homoge¬ 
neous  gaussian  Schell-model  source,  tVJ  increases 
monotonically  with  increasing  kog  and  asymptoti¬ 
cally  approaches  the  '  alu1'  unity  as  kog-*oo.  This 
limiting  case  corresponds  to  the  situation  where  the 
field  generated  by  the  source  coincides  with  a  wave- 
front  of  a  plane-wave  field  that  propagates  in  the 
positive  z-direction. 


3.3.  The  quasi-homogeneous  limit  (o,>og) 


3.1.  The  coherent  limit  (kag- >oo) 

When  the  source  of  the  class  that  we  are  consider¬ 
ing  is  completely  coherent,  kag—ao  and  eq.  (12) 
implies  that  {-»v^2/c(t,.  The  expression  (14)  for  the 
radiation  efficiency  then  becomes 

t,„  =  \-D(j2kol)lJlkol  .  (16) 

Since  the  second  term  on  the  right  of  eq.  (16) 
approaches  zero  as  ko/-+oo  (see  fig.  5)  we  see  that 
the  radiation  efficiency  of  a  coherent  source  then 
approaches  the  value  unity.  The  formula  (16)  applies 
to  certain  types  of  lasers  operating  in  their  lowest- 
order  mode. 

3. 2.  A  homogeneous  Schell-model  source  (ka,-*<x>) 


When  ko,$>kog  a  gaussian  Schell-model  source 
reduces  to  a  gaussian  correlated  quasi-homogeneous 
source  with  a  gaussian  intensity  profile.  The  radia¬ 
tion  efficiency  of  such  sources  was  shown  by  Carter 
and  Wolf  [  1 )  to  be  given  by 


D(kog 


(19) 


It  is  clear  that  our  expression  (14),  together  with  eq. 
(12),  indeed  reduce  to  eq.  ( 1 9 )  in  this  limiting  case. 

We  may  also  consider  the  limiting  case  of  a  com¬ 
pletely  coherent  quasi-homogeneous  source  by  let¬ 
ting  kog->co,  k<Ji~*cc  with  ko, / kog= const. $>  1.  Since 
£>(£)/<!; -* o  as  £— oo  it  follows  from  eq.  (19)  that  in 
this  limit 


Another  interesting  limiting  case  is  obtained  by 
letting  A<7,-»ao  (with  k  being  fixed),  and  kog  having 
an  arbitrary  but  fixed  value.  Eq.  (6a )  reduces  to  /(r, 
fu)  =  /0,  and  if  we  also  make  use  of  eq.  (6b)  the 
expression  (1 )  for  the  spectral  density  of  the  source 
becomes 


(,u  - 1  •  (20) 

Hence  the  radiation  efficiency  of  a  coherent  quasi- 
homogeneous  source  is  unity. 

Finally  we  deduce  from  eqs.  (16)  and  (19),  if  we 
recall  once  again  that  D(£)/£  decreases  monotoni¬ 
cally  with  increasing  £,  that  when 
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(ff*)qh  >2(®r)coh  (21) 

the  radiation  efficiency  fat  frequency  to)  of  a  gauss- 
ian  correlated  quasi-homogeneous  source  will  be 
greater  than  that  of  a  completely  coherent  Schell- 
model  with  gaussian  intensity  profile.  As  we  have 
mentioned  above,  certain  types  of  lasers  operating  in 
their  lowest  order  mode  correspond  to  a  coherent 
Schell-model  source  with  a  gaussian  intensity  pro¬ 
file.  It  is  therefore  clear  from  eq.  (21)  that  a  gauss- 
ian-correlated  quasi-homogeneous  source  may  have 
higher  radiation  efficiency  than  a  coherent  laser 
source  emitting  radiation  of  a  gaussian  intensity 
profile. 


Appendix 

Proof  that  i,0  <  I  for  planar  Schell-model  sources 

It  was  shown  in  ref.  [  1  ]  that  for  quasi-homogene¬ 
ous  sources  <*^1.  We  will  now  show  that  this 
inequality  holds,  in  fact,  for  all  planar  sources. 

We  start  by  showing  that  W{f,  -/,  w)  >  0  for  all  real 
two-dimensional  vectors  /(0<  \f\  <oo).  The  cross 
spectral  density  Wfr,,  r>,  at)  is  known  to  be  non-neg¬ 
ative  definite  [9]  i.e. 

J  j  Wfr,,  r2,  w)/(r,  )/*(r2)  d3r,  d2r:  ^0  ,  (Al) 

with  any  arbitrary  function  f(  r)  for  which  the  double 
integral  converges.  Let  us  choose  /( r)  =  exp(  -ifr). 
The  inequality  (Al )  then  gives 

jj  fV(r,,r2,(o) 

Xexp[ -i/(r, -r2)l  d2r,  d2r2  2.0 ,  (A2) 

which  implies  at  once  [cf.  eq.  (3)]  that 

,-ksL ,  w)s0  (OSk)s,  |  Soo)  .  (A3) 

If  we  substitute  eq.  (4)  into  eq.  (5)  and  use  the 
expression  (2)  for7,„(s)  we  obtain  the  following  for¬ 
mula  for  the  radiation  efficiency  . 

e,„  =  (2n/c)2  |  cos20#'(fcj,  ,-fca,  ,o;)dfl 
( 2») 


j  /(r,ct>)d2rj  ( A4) 

If  we  next  make  use  of  the  relation  cos20 
di2=;  ( 1  -  s2  -  s2) 1/2  dsr  ds,.  in  eq.  ( A4)  and  recall  that 
the  intensity  /(r,  a>)  =  W(r,  r,  to),  we  find  that 

t,uS(2nk)2  Jj  W(ksl,-ksi.,(o)  ds,  ds, 

s]  +52  i,  l 

X^  j  W(r,  r,  o>)  d2r^  (A5) 

In  view  of  the  inequality  (A3)  we  may  replace  the 
integration  over  the  unit  circle  in  the  numerator  of 
eq.  (AS)  by  integration  over  the  whole  sx ,  s,, -plane. 
After  doing  so  we  substitute  for  W(ks , ,  —  ks , ,  w) 
from  eq.  (3)  and  interchange  the  orders  of  integra¬ 
tions.  We  then  obtain  the  inequality 

00 

<„,S  Jj  H''(r1,r2,w)  <5(r,  -r2)  d2r,  d2r2 

-  «> 

*(  1  r,  w)d2r^  ,  ( A6) 

where  5  is  the  Dirac  delta  function.  On  carrying  out 
the  trivial  integration  with  respect  to  r2,  we  finally 
obtain  the  inequality 

<<„S1  ,  (A7) 

valid  for  all  planar  sources. 
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Mueller  matrices  and  depolarization  criteria 

U.  SIMON 

The  Institute  of  Mathematical  Sciences,  Madras  600  1 13,  India 
( Rt  i  eii  i  d  l(i  Scploithi-r  I  V.X6) 

Abstract.  The  question  of  whether  a  given  Mueller  matrix  represents  a 
determinism  or  a  non-deterministic  system  is  analysed  by  means  of  a  matrix 
condition.  The  possibility  of  replacing  this  matrix  condition  by  a  scalar  condition 
is  examined.  It  is  shown  that  this  is  permissible  only  for  those  eases  where  a 
Mermitian  matrix  constructed  from  the  Mueller  matrix  is  positive  semidefinite. 

1.  Introduction 

Several  methods  have  been  used  in  the  description  of  the  polarization  state  of  a 
waveHeld.  While  the  Jones  method  [  1 1  and  the  Poincare  sphere  method  1 2)  are  useful 
for  the  description  of  fully  polarized  states,  the  coherency  matrix  method  (3]  and  the 
Mueller  Stokes  method  (4,  5]  can  handle  both  partially  and  fully  polarized  light.  It 
should  he  noted  that  all  these  methods  assume  the  radiation  Held  under  consideration 
to  he  an  ensemble  of  plane  waves  all  having  the  same  wave-vector.  It  is  only  relatively 
recently  that  a  systematic  procedure  for  handling  polarization  in  a  beam  Held  has 
been  developed  (6|.  In  the  following  we  assume,  however,  the  radiation  Held  is  of  the 
former  type. 

The  coherency  matrix  and  the  Stokes  vector  are  equivalent,  and  carry  exactly  the 
same  amount  of  information.  However,  when  the  passage  of  the  beam  through  an 
optical  system  is  encountered,  the  situation  becomes  quite  different:  the  usual 
transformation  law  of  the  coherency  matrix  via  the  Jones  matrix  of  the  optical 
element  corresponds  to  deterministic  (non-depolarizing)  systems;  while  the  trans¬ 
formation  of  the  Stokes  vector  through  the  Mueller  matrix  corresponds  to  more 
general  systems  including  non-deterministic  (depolarizing)  systems.  In  the  deter¬ 
ministic  case  the  Mueller  matrix  can  he  derived  from  the  Jones  matrix  of  the  system. 
A  non-deterministic  system,  on  the  other  hand,  has  a  well-defined  Mueller  matrix; 
hut  there  does  not  exist  a  Jones  matrix  from  which  it  can  he  derived.  This  is  to  he 
expected,  for  the  Jones  matrices  form  a  seven-parameter  family  (the  absolute  phase 
of  the  Jones  matrix  should  he  suppressed  in  am  comparison  with  the  Mueller  matrix 
since  it  does  not  allect  the  transformation  ol  the  coherency  matrix,  this  transform¬ 
ation  being  quadratic  in  the  Jones  matrix),  whereas  the  Mueller  matrices  form  a 
sixteen-parameter  family. 

In  \  iew  ol  this  situation  the  following  question  is  ul  much  practical  interest .  I  low 
can  one  determine  whether  an  experimentally  measured  Mueller  matrix  corresponds 
to  a  deterministic  or  a  non-deterministic  system?  This  question  was  first  poseil  and 
examined  by  Harakat  |7|  A  complete  answer  to  this  question  in  the  form  of  a 
neeessary  and  sufficient  matrix  condition  w  as  subsequently  presented  by  the  present 
author  |X|.  (ill  anti  Bernabeti  J *> J  have  recently  made  the  interesting  claim  that  this 
matrix  condition  can  he  replaced  In  a  scalar  condition.  In  the  present  paper  we 
analyse  this  claim  and  show  that  it  is  not  valid  for  all  situations. 
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In  $  2  wo  briefly  rooall  tho  relationship  between  the  descriptions  in  terms  of  tin- 
coherency  matrix  and  in  term.s  of  the  Stokes  vector  leading  to  the  matrix  condition 
IS],  Then  we  analyse  the  scalar  condition  of  Gil  and  Bernabeu  and  show  that  it  is  not 
equivalent  to  the  matrix  condition  in  general.  In  fact  we  show  that  it  is  only  in  those 
situations  where  a  particular  I  lermitian  matrix  constructed  from  the  Mueller  matrix 
is  positive  semidertnite  that  the  scalar  condition  is  equivalent  to  the  matrix  condition. 
In  §  3  we  present  a  simple  example  which  illustrates  these  results.  Section  4  contains 
some  concluding  remarks. 


2.  Jones  matrix,  Mueller  matrix  and  the  depolarization  criterion 

The  coherency  matrix  ip  describing  a  polarization  state  is  a  2x2  complex 
Hermitian  positive  semidetinite  matrix: 


•P  II  <P  12 


</  =  </>. 


(p  ^0. 


l'/>2l  </>22j 

Its  transformation  by  a  deterministic  (non-image-forming) 
Jones  matrix  J  is  given  by 


opt  ical 


(I) 

system  with 


(2) 


For  the  purpose  of  comparison  with  the  Mueller  Stokes  formalism  it  is  convenient 
to  associate  with  every  coherency  matrix,  </>,  a  four-element  column,  <t>,  in  the 
following  one-to-one  manner: 


>11 ' 

■o>„- 

‘Pi  2 

<P  2  1 

<t>. 

_‘P  2  2  _ 

L(l>( 

(3) 


The  Stokes  vector  .S’  describing  the  same  state  is  related  to  <1>  through  a  simple 
numerical  matrix  A.  We  have 


S 


(4) 


where 

'  i  o  o  r 

I  II  0  I 

I  (5) 

II  1  1  0 

(I  i  i  u 

It  can  he  easilv  checked  that  I  is  unitary,  except  lor  a  multiplicative  factor,  and  we 
have 


I  1  =  i.  IV  mi 

Since  .1  is  non-singular,  it  follows  from  (4)  that  the  Stokes  \  e>  tor  and  the  <  oherenev 
matrix  arc  in  one-to-one  correspondence,  and  hence  lontain  identical  information 
about  the  state  of  the  field.  Since  <p  is  Hermitian  .S'  is  real,  and  tin  positive 
semidefinitem-ss  of  </>  implies 

•S'f,  -S'T  -S'j  .V„  -n  1 7 ) 
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l  'tuler  the  action  ot  an  optical  element,  the  change  in  the  polarization  state  is 
ilesi Tilled  through  a  linear  transformation  M  on  .S': 


.S’ 

-* ,S"  = 

MS, 

-v/„, 

A/„, 

A/,,, 

A/,„ 

A/„ 

Ul2 

A/,., 

A/.„ 

a/2, 

A/,, 

A/,., 

A/  m 

a/„ 

A/.,, 

A/,, 

The  4  x  4  real  matrix  M  is  called  the  Mueller  matrix  of  the  optical  element.  In  the 
special  case  where  the  optical  element  under  consideration  is  deterministic  it  can  he 
described  either  through  a  Jones  matrix  .1  or  a  Mueller  matrix  M,  and  the  two  are 
related  through  |  10,  X] 

A/--.l(./(x)./*).J  ',  CM 

where  *  indicates  complex  conjugation  and  (x)  denotes  the  Kronecker  matrix 
product. 

In  [8]  we  defined  a  matrix  ,Y  through  the  elements  of  M.  This  is  shown  as 
equation  (1(1)  on  the  following  pane.  This  relationship  between  M  and  N  is  clearly 
one-to-one.  The  matrix  .V  is  manifestly  llermitian.  Its  trace  is  simply  related  to  M: 

Tr(,V)  =  2iW00  (11) 

In  the  following  we  will  need  to  use  another  relationship  between  M  and  .Y: 

Tr(,Ya)  'Vr(MM'),  (12) 

where  A/  denotes  the  matrix  transpose  of  A/.  The  relationship  (12)  can  easily  be 
verified  trom  the  explicit  form  of  Y  given  in  (10)  by  noting  that  the  left-hand  side  of 
(12)  is  the  sum  of  the  modulus  square  ol  all  the  I (i  elements  of  .Y,  In  v  irtue  of  the 
I  termitian  property  ol  .V;  v\  tide  the  right-hand  siile  is  the  sum  of  the  squares  of  the 
elements  of  the  real  matrix  A/. 

In  the  spirit  ol  (5)  we  write  the  2x2  complex  Jones  matrix  ./  in  (2)  as  a  tour- 
element  column  \  cctor. 


I.veti  though  we  use  the  same  symbol  ,/  for  the  2x2  matrix  as  for  the  four-element 
column,  no  confusion  is  expected  to  arise.  For  deterministic  systems  whose  Mueller 
matrix  is  related  to  ./  through  CM  the  matrix  ,Y  is  related  to  ./  in  a  simple  wax  |X|: 

■V,„  7.  /(  M.I.2.T  (14) 
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Now  assume  thill  the  optical  element  is  deterministic.  Then  it  lias  a  Jones  matrix 
J ,  and  the  \  matrix  ol  the  system  is  given  In  ( 1 4).  Squaring  the  matrix  equation  (14) 
we  have 


■■|Tr(iY)JAV  (15) 

That  is,  lor  deterministic  systems 

.Y:  =  [Tr(.Y)).Y.  (Ih) 

Conversely,  assume  that  the  optical  element  satisfies  (  Hi).  That  is,  from  the  Mueller 
matrix  ol  the  given  optical  element,  we  construct  the  .Y  matrix  according  to  the 
prescription  (10),  and  this  matrix  satisfies  ( 1(>).  Then  ( 1  (>)  implies  that  |Tr  (,Y)|  1  ,Y 
is  a  projection  operator,  and  hence  .Y  can  he  written  in  the  form  (14)  for  some  J.  In 
other  words,  the  system  can  he  described  through  a  Jones  matrix  and  lienee  is 
deterministic.  Thus  we  have  the  following  theorem  |X]:  The  necessary  and  sufficient 
condition  for  an  opticid  system  nilh  a  nii  en  Mueller  matrix  to  hr  deterministic  is  that  its 
\  matrix  formed  throuf’h  (1(1)  should  satisfy  the  matrix  condition  (10). 

1  laving  established  ( 1  <>)  we  are  ready  now  to  analyse  the  results  of  other  authors 
in  the  light  of  this  result.  The  matrix  condition  of  liarakat  will  not  he  analysed  here 
(see  |S|).  Assume  that  we  hav  e  a  detv  rniinistic  system.  Then  (  . .  )  is  satisfied.  Taking 
tlie  trace  of  (  Ui)  and  using  (II)  we  obtain  for  such  systems 

Tr(.Y2)  =  4.l/2(l,  (17) 

ami  hence  from  (12), 

Tr  (. \/.UT  )  =  4. Uf, „.  (IS) 

I  hus  ( I  H)  is  a  necessary  condition  tor  a  sv stem  to  he  deterministic.  I  lence  the  result 
ol  T  ry  anil  Kattaw ar  |  I  I  j  is  consistent  with  out  matrix  condition.  ( ill  and  Bernabeu 
have  claimed  that  it  is  also  the  sufficiency  condition.  To  see  if  this  is  so  we  have  to 
examine  vv liether 

Tr(.YJ)  -  |Tr(.Y)|-  (I'D 

is  equivalent  to  (l(>).  Clearly,  there  arc  two  cases  to  In-  distinguished: 

Case  I.  \  is  post  tire  semidefinite 

In  this  ease  it  can  he  seen  that  (  I’J)  is  indeed  eqmv  alcnt  to  (  1 (>)  Tills  is  most  easilv 
established  by  recalling  that  A  is  llermilian,  and  working  in  its  diagonal 
rep  resentation. 

(  use  ..’  Y  is  not  positive  semidelimte 

In  tin-  v  use  l  1 1 1 1  implies  (  I ') ).  vv  livrcas  (  I'M  d<  >es  not  implv  (  I  fi  I  This  too  i-  easilv 
seen  in  the  diagonal  representation  ol  A. 

'll  Ills  till*  (|IU  >t|nl)  is  l  edlU  I'd  to  olH1  of  wllethei  tlii'li  t  \|sf  \liiclli  I  > 

"  hosi  \  matru  t  s  w  ill  ha\ e  at  least  one  tu^alix  e  e  mem  a  hie.  Such  Mueller  mat ru  es 

110  indeed  exist,  .mil  it  is  preetselx  for  these  th.it  the  i  hum  of  (  ill  and  Mem. thru  breaks 

Mow  ti  \\  e  i»i\  e  examples  nt  such  matt  u  es  in  the  I  ol  low  mu  turn  Mut  lute  \\e  note 
that  the  scalar  eotuhtion  (INK  <ir  eipm  .dent  )\  does  not  replaee  the  matnx 

i  omlit  ion  <  I  f\)  iii  \  tew  of  tlu-  tai  t  th.it  the  \  tn.it  t  tx  1^  not  ret  pi  tied  1 1  *  he  posit  i  \  e 
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3.  Example 

As  a  simple  example  to  illustrate  our  results  in  the  last  section,  consider  the 
matrix 


M 


0  0  0 

I  0  0 

I)  I  0 

I)  0  -  I 


(2D) 


This  is  a  valid  Mueller  matrix.  The  conditions  (7)  mean  that  the  Stokes  vector  should 
he  a  ‘time-like-  vector  with  positive  time-component.  Hence  any  4x4  real  matrix 
w  hit’ll  maps  every  'time- like’  vector  w  ith  positive  time-component  into  a  r  ector  with 
these  properties  is  an  acceptable  Mueller  matrix,  .1/  in  (2(1)  clearly  meets  this 
rei|uirement.  1'nrm.illv,  these  conditions  are  ideniital  to  those  imposed  on  proper 
I.orcntz  t  runs  I  urination.  I  nil  now  there  is  no  restriction  o!  invariance  oi  the  ‘norm’  of 
the  vector. 


The  A  matrix  corresponding  to  this  V/  is 


.V 


(I  (ID  1 
D  I  D  D 

D  D  I  (1 

I  D  D  D 


(2!) 


It  is  easily  checked  that  (2D)  satisfies  (IS)  and,  c(|ili  valent  l\  ,  (2  1  )  satisfies  ( I D).  Yet ,  it 
does  not  sal isl y  the  matrix  condition  (!(>)  and  hence  does  not  represent  a 
deterministic  system  In  tact  there  exists  no  .1  matrix  Iron)  which  the  M  can  hr 
derivevl  in  the  lorm  (')).  This  simple-lookiim  Mueller  matrix  which  changes  neither 
(hr-  intensitv  nor  the  deitrec  ot  polarization  of  am  input  state  is  non-determimstic, 
for  its  \  matrix  is  not  positive  semulefinite;  if  has  eigenvalues  (I,  1,  I,  I). 

As  vet  another  simple  example  we  rile  the  matrix 


l  / 


I  D  I)  (I 

D  I  (I  D 

D  II  |  II 


[  II  (1  (I  ! 


(22) 


The  rvi  id.  r  can  verilv  that  it  icpies.nts  a  uon-detetministic  svstcni  even  thoutth  it 
satisfies  the  scalar  condition 


4.  Concluding  remarks 

We  have  altalvsed  the  tniidilions  t . » r  a  M’leller  maliix  corresponds  to  a 
dctermmixtn  system  The  lie.  cssarv  and  siitlicielit  coiidittoo  tor  this  is  tliv  en  h\  the 
matrix  ei|  nation  I  I  ID  In  situations  m  w  Inch  the  cik>eu\  allies  of  the  .Ymaliix  I  or  Died 
tiom  the  unco  Mueller  maliix  an  all  noii-netialiv  >  .  and  only  under  these 
i  in  uni  slant  is,  is  (he  mati  ix  t  i  mdition  i  I  (,  t  et  pin  a  lent  to  the  scalar  com!  it  ion  ( IN  1 . 

Ilius  the  M  uc  I  let  mall  lies  divide  naturaliv  into  I  vv  o  dis|oint  classes:  one  with 
positive  si  niidetiiiile  \  matrix  and  the  other  with  an  \  matrix  that  has  at  least  one 
i  tee  at "  e  ei  eel  i  v  aim.  1 1  is  t  it  Intel  est  to  note  that  .1  von.l  it  loll  similar  to  (  I  fi  I  is  alrea.lv 
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known  m  tlu  context  ot  dynamical  mapping  of  the  density  matrix  of  a  quantum- 
meehanieal  system  |12|. 

\  related  h-'iir  ol  interest  is  the  possibility  ,  or  otlu  rw  ise.  of  realizing  the  Mueller 
m.itux  ot  .1  noh-delerimmstu  system  as  an  ensemble  ol  Mueller  matrices  of 
deter immsl u  systems.  I  Ins  problem  has  been  examined  in  a  reeent  paper  of  Kim, 
\  landel  and  \\  < ill  |  I 5  |. 
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The  effect  of  a  linear  random  medium  on  the  state  of  polarization  of  the  transmitted  light  is  investigated,  and  the 
connection  between  the  Stokes  vector  formalism  and  the  coherence  or  polarization  matrix  formalism  is  discussed. 
It  is  shown  that  an  ensemble  of  Jones  matrices  corresponds  to  the  Mueller  matrix  in  general. 


1.  INTRODUCTION 

When  light  propagates  through  a  linear  medium,  its  polar¬ 
ization  propert  ies  are  usually  described  either  by  the  Stokes 
vector  formalisin'  or  by  the  coherence  matrix  (also  known  as 
the  polarization  matrix)  that  was  introduced  by  Wiener'2  and 
by  Wolf.’1  The  effect  of  many  non-image-forming  optical 
devices  on  the  light  is  then  to  transform  both  the  Stokes 
vector  and  the  coherence  or  polarization  matrix,  so  that  the 
device  can  lie  represented  by  a  transformation  matrix.  This 
transformation  is  usually  known  as  the  Mueller  matrix4 
when  it  acts  on  the  four-dimensional  Stokes  vector  or  as  the 
Jones  matrix’’  when  it  acts  on  the  2X2  polarization  ma¬ 
trix.8-7 

Even  though  there  exists  a  one-to-one  correspondence 
between  a  polarization  matrix  and  a  Stokes  vector,  the  de¬ 
scription  of  optical  systems  in  terms  of  Mueller  matrices 
appears  to  be  applicable  to  more  general  situations  than 
does  the  description  in  terms  of  Jones  matrices.  This  was 
already  pointed  out  by  Azzam  and  Bashara,8  and  Howell9 
has  shown  that  some  optical  devices  can  be  described  by 
Mueller  matrices  but  not  by  Jones  matrices. 

In  several  recent  publications  the  constraints  that  must  be 
satisfied  for  a  Mueller  matrix  to  correspond  to  a  Jones  ma¬ 
trix  were  investigated.'"  14  Simon'2  and  Barakat"  found 
nine  constraints  that  are  necessary  if  a  Mueller  matrix  is  to 
be  derivable  from  a  single  Jones  matrix,  and  more  recently 
Gil  and  Bernabeu14  found  a  single  condition  on  the  trace  of 
the  square  of  the  Mueller  matrix.  These  results  apply  to 
propagation  through  a  deterministic  optical  device. 

On  the  other  hand,  some  optical  systems  are  nonde'  ermin- 
istic,  and  they  can  be  represented  by  an  ensemble.  In  what 
follows  we  show  that  when  an  ensemble  of  transformations  is 
introduced  to  describe  certain  stochastic  non-image-form- 
ing  optical  systems,  the  two  descriptions  can  be  completely 
reconciled,  and  both  are  equally  general.  In  the  special  case 
when  the  ensemble  reduces  to  a  single  realization,  the  trace 
condition  of  Gil  and  Bernabeu  follows  naturally. 

In  Section  2  we  review  the  properties  of  polarization  ma¬ 
trices  and  of  Stokes  vector,  and  in  Section  3  we  describe  the 
mathematical  transformations  that  characterize  transmis¬ 
sion  through  a  deterministic  device.  In  Section  4  we  intro 
duee  an  ensemble  of  transformations  to  represent  a  random 
linear  device,  and  we  examine  the  corresponding  relation 
between  the  Jones  and  Mueller  matrices. 


2.  THE  COHERENCE  OR  POLARIZATION 
MATRIX  AND  THE  STOKES  VECTOR 

We  consider  an  optical  field  in  the  form  of  a  quasi-mono- 
chromatic  plane  wave  propagating  in  some  direction  charac¬ 
terized  by  the  unit  wave  vector  x,  say,  the  z  direction.  Let  E 
be  the  complex  analytic  signal  representing  the  transverse 
vector  field.  The  field  can  always  be  resolved  into  two 
orthogonal  components,  1  and  2; 

E  =  E,tj  +  E.jt2,  (1) 

where  tlt  <2  are  orthogonal  unit  vectors.  «i,  t2  could  be  real 
unit  vectors  in  the  x,  >•  directions,  corresponding  to  orthogo¬ 
nal  linear  polarizations.  However,  sometimes  it  is  more 
convenient  to  resolve  the  field  into  more  general  orthogonal 
states  of  elliptic  polarization,  in  which  case  1 1 ,  <2  are  complex. 
In  any  case  the  transversality  of  E  is  expressed  by  the  condi¬ 
tion 

*.t,  =  0  (i  =  l,2)  (2) 

and  the  orthonormality  of  <■,  t2  by 

\ ,  (i.y-1.2),  (3) 

If  the  field  is  fluctuating,  then  E 1,  E2  in  Eq.  (1)  are  random 
variables  described  by  an  ensemble,  which  we  shall  assume 
to  be  stationary.  The  2  X  2  polarization  matrix  J  is  the 
covariance  matrix  of  the  two  variates  E,t  E>  and  is  given  by 

■/„  =  <£,£,*>•  <4> 

where  ( >  denotes  the  ensemble  average.  By  definition,  J  is 
Hermitian  and  nonnegative  definite,  and  its  trace  is  a  mea¬ 
sure  of  the  mean  light  intensity  (E*  -E).  The  effect  on  J  of 
changing  from  one  set  to  another  set  of  base  vectors  <i  and  t  j 
is  describable  by  a  unitary  transformation  on  ,J.  It  follows 
that  there  always  exists  a  basis  <1,  (■  in  which  J  is  diagonal, 
because  every  Hermitian  matrix  can  be  diagonalized  by  a 
unitary  transformation. 

The  degree  of  polarization  P  of  the  light  can  be  expressed 
either  in  terms  of  eigenvalues  of  J  or  in  terms  of  the  unitary 
invariants  of  J  in  the  form"’ 

/’=  1 1  -  4  del  J/fTr  J)2]u~.  (5) 

It  follows  from  either  form  that  0  <  P  <  1  and  that  P  =  1 
when  det  J  =  0  or,  equivalently,  when  one  eigenvalue  of  J  is 
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zero,  so  that  only  one  kind  of  polarization  is  present.  On  the 
other  hand,  P  -  0  when  the  two  eigenvalues  are  equal,  and  J 
is  then  proportional  to  the  unit  matrix,  corresponding  to  an 
equal  mixture  of  both  polarization  components.  It  may  be 
shown16  that  any  polarization  matrix  J  can  be  uniquely  de¬ 
composed  into  a  fully  polarized  part  and  a  fully  unpolarized 
part. 

The  four  elements  of  J  can  also  be  used  to  construct  four 
real  parameters  known  as  the  Stokes  parameters,1  which  are 
given  by 

=  (£,£,*>  +  < E2E ,*), 

S,  =  <£,£,»>  -  {E2E.,*), 

S2  -  (ElE.,m'>  +  (E 2E  *\ 

S:,  =  /«£,£•,*>  -  (t'.ft.,*) |  (6) 

and  also  represent  the  state  of  polarization  of  the  field.  The 
four  parameters  are  often  considered  to  be  the  components 
of  a  four-vector  S,  known  as  the  Stokes  vector.  In  terms  of 
the  components  of  S  the  degree  of  polarization  is  then  given 
by 

P  =  (S,'  +  S.?  +  .S/’)1  "’AS,,.  (7) 

Another  connection  between  the  Stokes  parameters  S„  (g 
=  0,  1 ,  2,  3)  and  the  polarization  matrix  J  becomes  apparent 
if  we  express  J  as  a  linear  combination  of  the  four  linearly 
independent  2X2  Pauli  spin  matrices 


winch  form  a  complete  set  for  the  representation  of  any  2X2 
matrix  Then  we  find  that 

./ .  *  '  .S>, l,i  =  ti,  1. 2  :tl,  (9) 

where  summation  on  repeuted  indices  is  understood.  It 
follows  that  the  Stokes  parameters  are  simply  twice  the 
coefficients  in  the  expansion  of  the  polarization  matrix  J  in 
terms  of  Pauli  matrices. 

Thts  immediately  leads  to  another  expression  for  ,S'U.  Let 
us  multiply  both  sides  of  Kq.  (9)  on  the  right  by  another 
Pauli  spin  matrix  a'  1  (r  =■  0.  1,  2.  3)  and  take  the  trace  on 
tioth  sides.  Then  we  obtain 

Tr|./rt'  '|  -  ]  Trln'-VM. 

Hecalling  that  the  product  of  two  different  Pauli  matrices 
.  n-lds  one  of  t he  t  hree  Pauli  matrices  V • a1  which  is 
fan-less,  we  see  that  the  only  none  mushing  contribution 
,,i.  .;rs  when  ,i  =  r,  in  which  case  l'  '.race  i  quills  2.  Hnd  we 
obtain  finally 
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■s;  =  'Yr\Ja'“'\  =  (,,=0.  1,2,3),  (10) 

when  the  summation  convention  is  applied, 

3.  TRANSMISSION  THROUGH  A  LINEAR 
OPTICAL  SYSTEM 

When  the  light  beam  is  passed  through  some  non-image- 
forming  optical  device  such  that  it  enters  and  emerges  as  a 
plane  wave,  the  new  field  vector  E\  given  by 

E' =  £,'<, -t- <2-  HD 

has  components  £•/  that  are  often  linearly  related  to  the 
old  components  E\,  E2-  For  brevity  we  shall  henceforth 
refer  to  this  device  as  a  filter.  We  may  then  represent  the 
filter  by  the  2  X  2  transformation  matrix  T.  usually  known  as 
the  Jones  matrix,  such  that 

E;=T„E,  (1.7  =  1.2),  (12) 

where  summation  on  rep.  ted  indices  is  again  understood. 
For  the  moment  we  take  the  ei,i  nts  of  T  to  have  definite 
values,  i.e.,  they  are  not  random,  implicit  forms  of  T  for 
certain  common  filters,  such  as  a  compensator  or  phase 
plate,  a  differential  absorber,  an  optical  rotator,  and  a  polar¬ 
izer,  have  been  given. I,i-1’: 

Let  us  now  examine  how  the  polarization  matrix  J  and  the 
Stokes  vector  S  are  affected  under  this  transformation.  We 
find  from  the  definitions  |Eqs.  (4)  and  (12)]  that  the  new 
polarization  matrix  </'  is  given  by 

J,/  =  <£,'£/*> 

=  T,JEmEn*)Tin’ 

=  TimJmnTn /  (13) 

or,  in  matrix  form, 

J’  =  Tjr,  (14) 

where  T*  is  the  Hermitian  adjoint  of  T.  Hence  J'  is  related 
to  ■)  by  a  similarity  transformation  involving  the  same  ma¬ 
trix  T  that  transforms  E  to  E'.  However,  as  we  show  below, 
there  exist  linear  filters  whose  effects  are  desoribable  not  by 
transformation  ( 12)  or  (14)  hut  only  by  an  ensemble  of  such 
transformations  " 

Let  us  now  examine  the  corresponding  transformation 
rule  for  the  Stokes  vector  S.  Under  any  linear  transforma¬ 
tion  the  new  Stokes  vector  S'  is  related  to  the  old  one  by 

S‘-M  S’,  (g,.=0,  1,2, 3).  (IS) 

The  t  x  1  transformation  matrix  M„,  is  known  as  the 
Mueller  matrix.'1  We  may  readily  obtain  the  form  of  A/„, 
when  the  field  vector  obeys  the  transformation  Eq.  (12)  by 
making  use  of  Eqs  ( 101  and  ( 14).  We  then  find  that 

■s;  =  TrlJ'ii'"') 

=  Tr|7V7Vl‘‘]  (16a) 

or.  in  component  form. 

v;  =  . <  itih) 

We  now  substitute  for  J  in  Eq.  ( IS)  from  Eq  (9>  and  obtain 


r 
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i 


S„'  =  |  TrlT«rt',7V',,lS,.  (17) 

Comparison  with  Eq.  (15  shows  that  in  this  case  the  Mueller 
matrix  is  given  by 

Af„„  =  ’  Tr|7VMr ff1"’!  =  -  Tr[o{“'T<>MP\ 

=  *  U8) 

and  it  is  evidently  related  to  the  2X2  Jones  matrix. 

It  is  apparent  from  Eq.  (18)  that  to  every  Jones  matrix  T 
there  corresponds  a  Mueller  matrix  M,  but  the  converse  is 
not  necessarily  true.  As  we  show  below,  there  are  physically 
realizable  but  nondeterministic  linear  filters  whose  effect  on 
the  polarization  matrix  is  not  expressible  in  the  form  of  Eq. 
(14).  although  the  Stokes  vector  transforms  as  in  Eq.  (15). 
In  particular,  under  the  similarity  transformation  (14),  an 
initially  polarized  light  beam  always  remains  fully  polarized, 
although  the  degree  of  polarization  of  a  partially  polarized 
beam  can  increase  or  decrease  on  transmission  through  a 
linear  filter. 

In  order  to  show  this  let  us  choose  the  polarization  basis  in 
which  the  original  polarization  matrix  J  is  diagonal.  If  the 
light  is  fully  polarized,  only  one  eigenvalue,  say  1\,  is  nonzero, 
and  J  must  be  of  the  form 

Jmn  =  M-Af  <«> 


By  using  Eq.  (18)  we  may  readily  derive  the  simple  trace 
condition  on  the  Mueller  matrix  that  has  been  shown  to 
apply  to  any  nondepolarizing  optical  system.14  If  MT  de¬ 
notes  the  transpose  of  M ,  then  from  Eq.  (18) 

Tr(M'M)  =»  MU,MU„ 

=  1  Tr(7V',7V',,rrr|7V,'7V’,'|.  (22) 

4 

where  we  have  made  use  of  the  fact  that  the  trace  is  invariant 
under  cyclic  permutation  of  factors.  With  the  help  of  the 
general  matrix  rules 

Tr|A)Tr|B)  =  Tr(A  ®  B]  (23) 

and 

(A  ®  B){C  ®  D)  =  AC  ®  BD,  (24) 

where  ®  denotes  the  direct  product,  we  can  reexpress  Eq. 
(22)  in  the  form 

=  1  TrW'T'o'*'  ®  TV-’TV"'] 

4 

=  1  Tr[(7V‘''  ®  T, r",')(7V,‘'  ®  7V“')1 

4 

=  -  Tr|(T  ®  Die1"1  ®  aM){P  ®  ®  (t"‘')|. 

4 

(25) 


Needless  to  say,  under  these  conditions  det  J  =  0,  and  from 
Eq.  (5)  it  follows  that  the  degree  of  polarization  P  =  1.  Let 
us  now  calculate  the  degree  of  polarization  P  of  the  light 
beam  emerging  from  the  linear  filter.  With  the  help  of  Eq. 
(19)  we  have,  from  Eq.  (13), 


From  the  explicit  form  [Eqs.  (8)]  of  the  Pauli  matrices  we 
find  that 


1 

2 


1  ®  = 


10  0  0 
0  0  10 
0  10  0 
0  0  0  I_ 


(26) 


so  that 

det  J'  =  /,<T1I7V7’,1T21*  -  r2ITu*rnT21*)  =  0.  (20) 

Hence  P  =  1,  which  implies  that  polarized  light  remains 
polarized  after  passing  through  any  linear  filter  whose  effect 
is  described  by  Eq.  (12)  or  (14).  Evidently,  a  depolarizing 
filter  is  excluded  from  this  category.  However,  the  action  of 
a  depolarizing  filter  on  the  Stokes  vector  S  is  still  describable 
by  a  transformation  of  the  form  shown  in  Eq.  (15),  although 
the  actual  transformation  matrix  M  is  then  no  longer  given 
by  Eq.  (18).  For  example,  for  the  fully  polarized  light  de¬ 
scribed  by  Eq.  ( 19)  the  Stokes  vector  S  has  components  (/ lt 
0.  0.  / 1 ),  whereas  the  Stokes  vector  for  unpolarized  light  is 
always  of  the  form  (/,  0,  0,  0).  The  4X4  transformation 
matrix  of  the  form 

K  °i  b>  K 
0  «■’  h.  o 

‘  -  (21) 
0  a{  b{  0 

a4  b4 

converts  (/ 1,  0,  0, 1 1 )  into  (Kl\,  0.  0,  0)  and  therefore  repre¬ 
sents  a  fully  depolarizing  filter  device.  The  extra  degrees  of 
f  reedom  available  in  M  permit  this  possibility,  whereas  there 
is  no  2  X  2  transformation  matrix  T  to  represent  this  filter. 
In  this  sense  the  Mueller  transformation  matrix  appears  to 
be  of  more  general  applicability  than  the  Jones  matrix. 


and  this  matrix  commutes  with  any  4X4  matrix  of  the  form 
T  ®  T.  Moreover,  its  square  is  the  unit  4X4  matrix.  It 
follows  that 

Tr(Af7M)  -  Tr|(T®  79(7*  ®  P)] 
and,  with  the  help  of  Eqs.  (24)  and  (23), 

Tr(MrM)  =  Tt(TP  ®  TT") 

=  [Tr(7T*)]2.  (27) 

But  from  Eq.  (18)  we  have 

Mm  =  l/2Tr(7’7n)  (28) 

so  that  finally 

TrtAfM1)  =  4A fw2.  (29* 

This  is  the  necessary  and  sufficient  condition  found  by  Gil 
and  Bernabeu14  for  a  Mueller  matrix  to  represent  a  nonde¬ 
polarizing  optical  system.  We  see  that  it  holds  for  every 
deterministic  optical  system  described  by  Eq.  (12). 

4.  REPRESENTATION  OF  A  FILTER  BY  AN 
ENSEMBLE 

So  far  we  have  considered  only  deterministic  optical  sys¬ 
tems.  But  in  some  situations,  for  example  when  light  is 
passed  through  the  atmosphere,  the  system  is  no  longer 
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deterministic  and  must  be  described  by  an  ensemble  of  fil¬ 
ters.8  We  shall  represent  a  typical  element  of  the  ensemble 
by  a  2  X  2  Jones  transformation  matrix  7V*  and  assume  that 
it  occurs  with  probability  p,..  The  action  of  the  ensemble  of 
filters  on  the  incident  light  E  is  to  generate  another  ensem¬ 
ble  of  vector  fields  E'(*'1  generally  with  different  polarization 
states,  such  that 

(30) 

We  emphasize  that  the  new  ensemble  associated  with  the 
filter  is  in  addition  to  the  ensemble  formed  by  the  various 
realizations  of  the  incident  field.  In  constructing  the  ele¬ 
ments  J of  the  polarization  matrix  of  the  light  that  has 
passed  through  the  optical  system,  we  then  need  to  average 
over  the  e  ensemble  also.  Thus  we  obtain  from  Eqs.  (13)  and 

(30) 


=  (31) 


M  3  1  Tnp,‘ma  ’“'cr  Ul.  (37) 

When  Eq.  (18)  is  rewritten  with  'T„pqm  in  place  of  the  prod¬ 
uct,  T„pTqm'  because  the  ensemble  has  only  a  single  mem¬ 
ber,  then  Eqs.  (18)  and  (37)  become  formally  identical. 
Similarly,  if  we  make  the  same  substitutions  in  Eqs.  (13)  and 

(31 ),  the  two  equations  become  indistinguishable.  We  have 
therefore  demonstrated  a  one-to-one  correspondence  be¬ 
tween  the  transformation  laws  for  the  Jones  and  the  Mueller 
matrices  and  have  shown  that  ‘Tnpqm  completely  determines 
both  transformations.  Any  single  realization  of  the  Jones 
matrix  T{e)  is  clearly  inadequate  to  describe  the  optical  sys¬ 
tem.  The  ensemble-average  product  is  needed  for  the  calcu¬ 
lation  of  quantit  ies  such  as  J,/  that  are  of  the  second  order  in 
the  field. 

Finally,  we  consider  the  problem  of  inverting  Eq.  (37),  or 
deriving  'Tnpqm  from  the  Mueller  matrix  A#„,.  For  this  pur¬ 
pose  we  use  Eq.  (37)  to  construct  the  following  sum  over  the 
indices  p,  >•: 


M... 


=  T 


OI_  (.1 
q  aht 


(38) 


with  summation  over  repeated  indices  again  understood. 
But,  from  definition  (8), 


or 


J'  =  V  (p(.T,,'’jr",”t  =  yTr'JV'"),.,  (32) 

where  < is  a  shorthand  notation  for  the  average  over  the  e 
ensemble.  This  result  should  be  compared  with  Eq.  (14). 
In  this  case  J‘  can  no  longer  be  related  to  J  by  a  similarity 
transformation,  as  (ora  single  realization  of  the  ensemble. 

Moreover,  because  of  the  e  ensemble,  it  is  no  longer  true 
that  polarized  light  passing  through  the  optical  system  re¬ 
mains  fully  polarized.  Thus,  if  J  is  given  by  Eq.  (19)  as 
before,  when  we  calculate  the  determinant  of  J"  from  Eq. 

(32),  we  find  in  place  of  Eq.  (20) 

det.7'  =  /,|<T1/‘''Tnl<'*>„<T2ll'l7',1'l'*>, 

-  (r,1"-lT11",|*)r(Tll,,-’n1,"l*)„  i,  (33) 


where  \  >,.  again  denotes  the  average  over  the  filter  ensemble. 
I  n  general,  det  J’  is  not  zero,  so  that  P  *  1  and  the  emerging 
light  is  no  longer  fully  polarized.  It  follows  that  our  ensem¬ 
ble  representation  of  the  filter  can  now  accommodate  a  de¬ 
polarizing  filter  also. 

Next  we  calculate  the  effect  of  our  optical  system  on  the 
Stokes  vector  S.  By  going  through  the  same  procedure  as  in 
the  derivation  leading  to  Eqs.  (17)  and  (18),  we  now  find  that 
the  Mueller  transformation  matrix  is  given  by 


p,V“lr'V"p'" 


or  more  explicitly  in  component  form  by 


(34) 


*  ,  7'  u-l'p  1 1- It , 
2 '  '!/•  ivm 


|„|  in 

mu  pq 


(35) 


If  we  denote  the  ensemble  average  of  the  product  of  two 
Jones  matrices  by 


T 


npqm 


then  Eq.  (35)  becomes 


<  T  '"T  '•"> 


(36) 


=  2  K,Kr  <39) 

When  this  result  is  used  twice  in  Eq.  (3$),  we  arrive  at 

2  =  Vie  <40> 

which  is  the  inverse  of  Eq.  (37)  and  shows  explicitly  that 
‘Tnpqm  =  (Tnp(rlTqnie,,)e  is  completely  determined  by  the 
Mueller  matrix  Afu,.  However,  there  is  no  unique  procedure 
for  constructing  the  ensemble  of  Jones  matrices  Tnpie),  ex¬ 
cept  in  the  degenerate  case,  when  Eq.  (29)  holds  and  the 
ensemble  redures  to  a  single  realization. 
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We  showed  in  a  recent  report'  (see  also  refs  2-4)  that  the  normal¬ 
ized  spectrum  of  light  will,  in  general,  change  on  propagation  in 
free  space.  We  also  showed  that  the  normalized  spectrum  of  light 
emitted  by  a  source  of  a  well-defined  class  will,  however,  be  the 
same  throughout  the  far  zone  if  the  degree  of  spectral  coherence 
of  the  source  satisfies  a  certain  scaling  law.  The  usual  thermal 
sources  appear  to  he  of  this  kind.  These  theoretical  predictions 
were  subsequently  verified  by  experiments5.  Here,  we  demonstrate 
that  under  certain  circumstances  the  modification  of  the  normal¬ 
ized  spectrum  of  the  emitted  light  caused  by  the  correlations 
between  the  source  fluctuations  within  the  source  region  can  pro¬ 
duce  redshifts  of  spectral  lines  in  the  emitted  light.  Our  results 
suggest  a  possible  explanation  of  various  puzzling  features  of  the 
spectra  of  some  stellar  objects,  particularly  quasars. 

To  explain  why  source  correlations  influence  the  spectrum  of 
the  emitted  light  consider  a  very  simple  example.  Suppose  that 
two  point  sources  P,  and  l\  have  identical  spectra  SQito\  and 
that  measurements  on  the  emitted  field  are  made  at  some  point 
l\  Ehe  sources  are  assumed  at  test  relative  to  an  observer  at  P. 
Assuming  that  the  source  fluctuations  can  he  descrihed  by  a 
stationary  ensemble,  the  field  at  P  may  he  characterized  by  an 
ensemble  {  V(  P.  ui  1}  of  frequency-dependent  realizations'1,  each 
of  the  form 


e  ’  e  ; 

I't  P,  u>  |  -  Q(  P(1  u)  |  •  <?(  P  .  w  I  (I) 

Rj  R, 

where  ( Qt  P„  to )),  (/  -  1,2),  characterize  the  strengths  of  the 
two  fluctuating  point  sources,  Rt  and  R.  are  the  distances  from 
P,  to  P  and  from  P,  to  P  respectively  (see  Pig.  I )  and  k  -  w/c, 
c  being  the  speed  of  light  in  vacuo.  For  simplicity  polarization 
eltects  are  ignored  and  hence  V'  and  Q  are  taken  to  be  scalars 
The  spectrum  of  the  light  at  P  is  then  given  by 


■S’l  I  P,  ui  I  —  ( V5 1  P,  cl  i  V(  P,  ui  I)  ( 2 ) 

where  the  asterisk  denotes  the  complex  conjugate  and  the 
angular  brackets  denote  the  ensemble  average.  On  substituting 
from  equation  I  1)  into  equation  (2)  and  using  the  (act  that 

(y‘(P„.UlC>(P„ft/l)  -  <y‘lP,  wlOtP,,  <»)>  .Sy(m)  ft) 


the  following  expression  is  obtained  lor  the  spectrum  of  the 
emitted  light  at  P: 


,,liK  #<  ,  t 

i  IV  ,i  p, ,  P,  a/ 1  i  v  c 

R,R 


1 4  l 


their  strengths 

A  generalization  of  the  elementary  formula  (4)  for  radiation 
from  three-dimensional  steady-state  (that  is,  statistically  station¬ 
ary)  sources  of  any  state  of  coherence  is  known  Of  special 
interest  in  the  present  context  is  the  form  that  the  formula  lakes 
when  the  source  has  the  same  normalized  spectrum  xtJ( ui  I, 

(  f*  ipl  w )  dru  =  1 1  at  each  point  in  the  source  region  and  has  a 
degree  of  spectral  coherence'  (appropriately  normalized  cross- 
spectral  density)  pi.^lr,,  r  ,  w )  that  depends  on  the  position  vec¬ 
tors  r,  and  r:.  of  any  source  points  P,  and  P.  only  through  their 
difference  r,  -r,.  If.  in  addition,  lor  each  frequence  that  sig¬ 
nificantly  contributes  to  the  source  spectrum,  the  spectral  corre¬ 
lation  length  (the  eUeclive  spatial  width  ' Ar'|  of  ( r',  uz  U ]  is 
small  compared  to  the  linear  dimensions  of  the  source,  (he 
normalized  spectrum  x'v'  '(u,  ml  of  the  emitted  light  in  the  far 
zone,  in  a  direction  specified  by  a  unit  vector  u,  becomes  isee 
equation  tj.l  I )  of  ref.  .X> 

v„( m  lu,,(  ku,  ui  I 

.  "  -  -  (61 
V(j(  UI  )gl,_,l  KU,  m  I  dui 

where  Py(K,  ml  is  the  three-dimensional  spatial  Fourier  trans¬ 
form  of  the  degree  of  spectral  coherence  gylf.u/l  ; 
Pi/tU'  ~  *i.  to). 

Let  us  now  choose  as  the  normalized  source  spectrum  v,_,(w 1 
a  spectral  line  with  a  gaussian  profile 

'v(o;l  exp  [  -  (m  -  ■  ai„l'72f>:]  Ifi"  to„l  l7l 

fis  2ir 

and  suppose  that  at  each  effective  frequency  m.  the  source 
correlation  decreases  with  the  separation  jr'l  ir,.  r,|ofany  two 
source  points  in  a  gaussian  manner,  lhat  is 

gu(r',  a/l  cxp|  r';/2ir'lw)|  (X) 

On  taking  the  Fourier  transform  of  equation  (Ki  and  substituting 
the  resulting  expression  into  equation  (6)  we  obtain  the  follow¬ 
ing  expression  for  the  normalized  spectrum  of  the  emitted  light 
in  the  far  zone  Isee  equation  I.V21)  of  ref.  XI 

\>( w  u'h  I w  )  exp  ( -}|  kirti  (<«  1 |'| 
i  k  { a)  )  -  r , -■  '  ■  - —  -  *  ,  (9  > 

j„  r(>(  ui  lir^ ( a/ )  exp  { -s[kirH(  m  1 1'(  dm 

Here,  ,r\'  '(in)  is  written  in  place  of  x'v'lu,  w),  because  the 
spectrum  of  the  far  field  is  now  independent  of  u,  as  a  con¬ 
sequence  of  the  assumed  isotropy  of  (see  equation  (8)1. 

The  formula  I1))  shows  that  the  spectrum  of  the  emitted  light 
in  the  far  zone  depends  both  on  the  spectrum  of  the  source 
fluctuations  and  on  the  manner  in  which  the  elfective  source 
correlation  length  <zu(m)  depends  on  the  frequency  m. 

Let  u.s  consider  two  particular  cases.  Ill  Suppose  tost  lhat 
ir„(u>)  is  independent  of  to.  Letting  y  denote  the  (now  constantl 
value  of  and  with  \,,( m I  given  by  equation  (7),  one  can 
readily  evaluate  the  integral  in  the  denominator  on  the  right  of 
equation  (9|  and  one  then  finds  that 

'(ml  ^  ^  c\p  ^  «i  ”'j  lllb 

where 


lu,  m) -= 


Here 


tv,,  I  P  ,  m  I  nr.  I\'  If  )l  />,  I, 


(X| 


C) 


I  I  I  .1  I 


is  the  so-called  ohm  spectral  density  of  the  source  fluctuations 
and  v  c  denotes  the  complex  conjugate 

the  formula  (4i  shows  that  the  spectrum  \  I  /’.  m )  is.  in 
general,  not  lust  proportional  lo  S,,lm'  but  is  modified  by  the 
uirrel.il ion.  chuiacten/vd  h\  tt ,  ,i  i,  between  the  fluctu¬ 

ations  ol  the  two  soilin'  strengths  (Ji  l\.  ml  ami  (Jl  V  .  ml.  Only 
in  some  vers  special  vases,  lor  example,  when  the  source  llticlu 
atoms  arc  unvorielated  |  Vt,,l  l\,  P.,  I  l)|  will  ,V,  (  P,  ml  he 
proportional  lo  .S,ymi  Hcikc,  in  general,  the  spectrum  of  the 
light  generated  bv  two  point  sources  depends  not  only  on  their 
spe-ura  but  also  mi  the  v orrclaloui  between  the  fluctuations  of 


and 


( I  Ibl 


When  the  source  is  ctlcctixely  spatially  incoherent,  v  ■*<'.  Lhen 
according  to  equation  ( 1 1 1  A  ♦  x  and  o  •  1  and  it  follows  (tom 
equations  t  lot  and  (’I  that  in  (his  ease 


I  m  I  •  x,  ,l  /,/ 1 


1 1 2  > 


Hence,  in  (he  limiting  case  of  a  completely  incoherent  source 
of  the  class  that  is  considered  here,  the  normalized  spectrum  of 
the  emitted  light  in  the  far  zone  is  identical  with  the  normalized 


spectrum  of  the  source  fluctuations. 

However,  when  the  source  fluctuations  are  correlated  over  an 
effective  distance  f  >0,  equation  (10)  shows  that  the  spectrum 
s'v ’(<»),  although  it  is  also  a  line  with  a  gaussian  profile,  is 
centred  at  a  lower  frequency  w'a=  toj a} <  oi0.  Hence  the  source 
correlations  give  rise  to  a  spectral  line  sVc,(<*>)  that  is  redshifted 
with  respect  to  the  spectral  line  produced  by  the  completely 
spatially  incoherent  source  with  the  source  spectrum  r0(<u).  The 
shifted  line  is  narrower,  having  root-mean-square  width  S'  = 
S/n<~  S  and  has  a- times  greater  height.  Examples  of  spectra 
of  light  in  the  far  zone,  produced  by  several  sources  which  emit 
the  same  spectral  line  but  which  have  different  correlation 
lengths  are  shown  in  Fig.  2.  From  the  formula  (10)  one  can 
readily  deduce  that  the  relative  shift  of  the  line,  namely, 

(13) 

Ao 

( A0  =  2nc/wu*  Aj,- 27tc/o>o)  is  given  by 


which  shows  that  in  this  case  the  redshift  increases  quadratically 
with  the  spectral  source-correlation  length  (.  (2)  Next  consider 
the  situation  when  (<o)  ■■=  a/oi  where  a  is  a  positive  constant. 
The  expression  (9)  for  the  normalized  spectrum  of  the  emitted 
light  in  the  far  zone  now  reduces  to 

I  15) 

Jii  l *</«>)/ at  1  dm 

Note  that  this  expression  is  independent  of  the  value  of  the 
constant  a. 

When  Sy(m)  is  a  line  with  a  gaussian  profile,  given  by  equation 
(7),  the  spectrum  s'vr,(<o),  given  by  equation  (15)  is  no  longer 
strictly  gaussian  but  it  can  be  closely  approximated  by  a  gaussian 
and  can  be  shown  to  be  redshifted  with  respect  to  sQ(u>)  by  the 
relative  amount 


An  example  of  this  situation  is  illustrated  in  Fig.  3. 

This  case  [<r„(wl  -  a /to  |  is  of  special  interest  because,  accord¬ 
ing  to  equation  (8),  the  degree  of  spectral  coherence  is  now 
given  by 

M<z< r  ,  w )  =  exp  [  - <  kr')’/ 2(  a/ <•)*’],  (17) 


Fig.  2  Spectra  s\' '( <a )  of  the  far  field  from  sources  with  spectrum 
Vj(o>)  -  (fiv^rr)  'exp|  (ui-  <«*1,)i/2fi'' )  and  degree  of  spectral 
coherence  uy(r',  u>|  -- expl  -  r'J/2{’),  with  u„  =  3.887  x  I01'  s  1 
(A„  =  4,861  A)  and  S  —  9. 57 x  I0”s  ’,  for  several  selected  values 
of  the  effective  source-correl.iiion  length  (.  The  solid  curve  ( (  -»0( 
also  represents  the  source  spectrum  s0(«u). 

for  the  spectrum  of  the  light  emitted  by  a  planar  secondary 
source  of  a  well-defined  class  to  have  certain  invariance 
properties  on  propagation.  It  will  be  shown  in  another  publica¬ 
tion  (J.  T.  Foley  and  E.  Wolf,  in  preparation)  that  for  three- 
dimensional  primary  sources  of  an  analogous  class,  whose 
degree  of  coherence  satisfies  this  law,  the  spectrum  of  the  emitted 
light  has  similar  invariance  properties.  We  conjecture  that  the 
usual  thermal  sources  obey  such  a  scaling  law. 

Now  briefly  consider  the  question  of  a  physical  mechanism 
for  producing  source  correlations.  Such  correlations  must  clearly 
be  manifestations  of  some  cooperative  phenomena.  At  the 
atomic  level  possible  candidates  may  perhaps  be  superradiance 
and  superfluorescence''.  An  effect  of  this  kind  was  first  predicted 
by  Dicke  in  1954  when  he  showed1"  that  under  certain  circum¬ 
stances  energy  from  excited  atoms  may  be  released  cooperatively 
in  a  much  shorter  time  than  the  natural  lifetime  of  the  excited 
states  of  the  atoms  and  with  much  larger  emission  intensity  than 
would  be  obtained  were  the  atoms  radiating  independently. 
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utdO'V) 


that  is,  it  has  the  functional  form 


ml  /lAr'l  (A 


I  bus  the  degree  ol  speeli.il  coin-mu  c  ol  the  source  dislribution 
now  satisfies  the  three  dimensional  analogue  of  a  requirement 
i  called  the  scaling  law  I  derived  in  ref.  I ,  as  a  sufficient  condition 


Hg.  .1  The  spectrum  s\'  ’(cm)  of  the  far  held  from  a  source  with 
source  spectrum  (An2tt)  1  exp  1  U»  w„»7?A')  and 

degree  of  spevlial  coherence  ^(r, ».»)  c*p|  i  kr  )' /  ?ia  / 1  )  | 
(<j  im  mhitrnry  constant),  with  •*»„  IKH7*l<l,,s  '  (A,, 

4, KM  A)  and  A  ¥.57x10' 1  I  he  source  spcctium  is 

shown  for  comparison  Note  that  w)  now  obeys  the  scaling 

law 
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However  not  enough  is  known  .it  present  about  the  coherence 
properties  of  large  three-dimensional  systems  of  this  kind  to 
make  it  possible  to  determine  whether  superradiance  and  super- 
fluorescence  might  involve  correlations  that  could  give  rise  to 
spectral  line  shifts 

I  here  is,  however,  c|uite  a  different  mechanism,  which  can  be 
described  at  the  macroscopic  level,  arid  which  can  imitate  eliecls 
of  source  correlations,  namely  rtlec  s  of  correlations  between 
the  refractive  index  at  pairs  ol  points  m  a  spatially  random  but 
statistically  homogeneous,  time  invariant  medium  If  a  wave 
illuminates  Mich  a  medium,  say  a  dilute  gas,  then,  as  is  well 
known,  the  medium  acts  as  a  secondary  source,  namely  as  a  set 
of  oscillating  charges  set  in  motion  by  the  incident  wave.  The 
secondary  waves  produced  by  the  oscillating  charges  then  com¬ 
bine  with  each  other  and  with  the  incident  wave  and  generate 
the  scattered  held  It  the  gas  is  not  too  dilute  the  collective 
response  of  the  microscopic  charges  to  the  incident  lield  can  be 
described  bv  macroscopic  paiameters  such  as  the  dielectric 
susceptibility  oi  the  refractive  index  Now  within  the  accuracy 
of  the  lirst  Born  approximation  the  basic  equation  for  scattering 
i.  of  the  same  form  as  the  basic  equation  for  radiation  from 
primary  sources,  the  ‘equivalent  source'  for  scattering  being  the 
product  of  the  scattering  potential  (which  is  a  simple  function 
of  the  refractive  index  l  and  of  the  incident  wave  This  correspon¬ 
dence  elearlv  implies  that  our  results  regarding  the  effects  of 
source  correlations  on  the  spectrum  of  the  emitted  light  must 
have  analogues  regarding  the  effects  of  a  spatially  random 
medium  with  correlated  refractive  index  distribution  on  the 
spectrum  ol  the  light  that  is  scattered  by  it.  This  topic  will  be 
discussed  elsewhere. 

I  et  us  now  consider  some  implications  of  this  analysis  Using 
equation  (  14),  (he  spectral  line  in  f  ig.  2,  produced  by  the  source 
whose  correlation  length  £  is  readily  found  to  have  a  redshilt 
given  by  :  0.0241  with  respect  to  the  source  spectrum.  It  is  of 

interest  to  note  that  if  an  observer  detected  such  a  redshilt 
unaware  of  its  true  origin  and  interpreted  it  on  the  basis  of  the 
Doppler  shill  formula  v/r  2k\/A„  :  he  would  incorrectly 

conclude  that  the  source  was  weeding  from  him  with  a  speed 
i  00241  <  7,2.10  kins  V 
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It  seems  worthwhile  to  note  that  there  is  a  maximum  line  shift 
that  can  be  produced  by  source  correlations  This  can  he  seen 
from  the  basic  formula  (6)  which  indicates  that  vV‘(u,  <al  0 
when  vy(<a  I  •  0,  implying  that  the  spectrum  of  the  far  field  can 
only  contain  those  frequencies  that  are  already  present  in  the 
source  spectrum.  Consequently  the  maximum  attainable 
frequency  shift  of  the  line  cannot  exceed  its  effective  frequency 
range.  However,  any  frequency  contribution  front  the  source 
spectrum  to  the  normalized  spectrum  of  the  far  field  can  be 
greatly  magnified  or  greatly  reduced,  as  is  evident  from  equation 
((>)  and  front  lag  2. 

We  have  mainly  considered  eliecls  of  source  correlations 
under  circumstances  when  the  source  spectrum  consists  of  a 
single  line  and  when  the  degree  of  spectral  coherence  /a ,,  that 
chaiacterizes  the  source  correlations  depends  on  a  single  para 
meter  Preliminary  calculations  show  that  with  a  suitably  chosen 
fi0  which  depends  on  a  larger  number  of  parameters,  redshifts 
of  several  lines  may  he  produced,  all  of  which  will  have  approxi¬ 
mately  the  same  e-values 

In  this  article  we  have  considered  redshifts  of  spectral  lines. 
However,  it  is  not  difficult  to  specify  source  correlations  which 
will  produce  blueshifts.  Examples  of  this  kind  are  given  in  a 
forthcoming  publication". 

It  seems  plausible  that  the  mechanism  discussed  in  this  article 
may  he  responsible  for  some  of  the  so  far  unexplained  features 
of  quasar  spectra,  including  line  asymmetries  and  small  dilfer- 
ertces  in  the  observed  redshifts  of  different  lines.  In  this  connec¬ 
tion  it  is  of  interest  to  recall  that  the  role  of  coherence  in  the 
emission  of  radiation  from  quasars  was  stressed  by  Hoyle, 
Burbidge  and  Sargent  in  a  well-known  article1'. 

I  thank  Mr  A  Gamlicl  and  Mr  K.  Kim  for  carrying  out 
compulations  relating  to  the  analysis  presented  in  this  article 
The  fact  that  scattering  can  also  produce  shifts  of  spectral  lines 
was  noted  independently  by  Professor  franco  (iori,  who  in¬ 
formed  me  of  this  result  when  commenting  on  an  early  version 
of  the  manuscript  ol  this  article.  1  his  investigation  was 
supported  by  the  NSf  and  by  the  US  Air  force  Geophysical 
Laboratory 
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li  has  recently  been  shown  theoretically  that  correlations  between  fluctuations  of  the  source  distribu¬ 
tion  at  different  source  points  can  produce  red  shifts  or  blue  shifts  of  emitted  spectral  lines.  To  facilitate 
experimental  demonstration  of  this  effect  a  simple  example  is  analyzed.  It  involves  only  two  small  ap¬ 
propriately  correlated  sources  and  brings  out  the  essential  physical  features  of  this  new  phenomenon. 


The  main  features  of  this  theoretical  prediction  have 
been  confirmed  by  Bocko.  Douglass,  and  Knox,  using 
acoustical  rather  than  optical  sources.  An  account  of 
their  experiments  is  given  in  the  accompanying  Letter.5 

Let  us  consider  two  small  fluctuating  sources  located 
at  points  P\  and  Pi.  I  assume  that  the  fluctuations  are 
statistically  stationary.  Let  [Q(P\h<o)\  and  IpfPz.tu)} 
be  the  ensembles  that  represent  the  source  fluctuations6 
at  frequency  to.  Furthermore,  let  {U(P,<u)\  be  the  en¬ 
semble  that  represents  the  field  at  point  P  generated  by 
the  two  sources  (Fig.  1).  Each  realization  C/(P,<o )  may 
then  be  expressed  in  the  form7 

(kRx  *kR2 

U(P,(0)-QiPu<o)L~  +  Q(Pi.a>)~-.  (I) 

A  1  A2 

where  R  i  and  Ri  are  the  distances  from  P\  to  P  and 
from  Pi  to  P,  respectively,  and  k  — to/c ,  c  being  the  speed 
of  light  in  free  space.  The  spectrum  of  the  field  at  the 
point  P  is  given  by 

S(,(P,w)-<(/*(P,w)(/(P,«)>,  (2) 

where  the  angular  brackets  denote  ensemble  average. 

J  On  substitution  from  Eq.  (I)  into  Eq.  (2),  we  find  that 

S(  (P,m)-(l/Pf+l/P^)S<?(w)  +  [JFe(P1,P2,w)e'*<''!',f')/RiR2  +  c.c.J.  (3) 

Here 

S(,(a>)-(Q'(Pu(o)Q(Pt,(o))-(Q*(Pi,w)Q{P2,a>))  (4) 

is  the  spectrum  (assumed  to  be  the  same)  of  each  of  the  two  source  distributions, 

Wq{P \,Pi,m)  r“(Q* (P\,w)Q(Pi,to))  (5) 

is  the  cross-spectral  density  of  the  source  fluctuations  (first  paper  of  Ref.  6,  Eqs.  (3.3)  and  (5.9)1,  and  c.c.  denotes  the 
complex  conjugate. 

The  degree  of  spectral  coherence  at  frequency  to,  which  is  a  measure  of  correlation  that  may  exist  between  the  two 
fluctuating  sources,  is  given  by  the  formula8 

ptf(P1,P2,cu)-W/(;(P|.P2,ft))/Se((u).  (6) 

The  normalization  in  Eq.  (6)  ensures  that  0<  \pq(P\,Pi,oj) |  —  1.  The  extreme  value  \pq\  “1  characterizes  com¬ 
plete  correlation  (complete  spatial  coherence)  at  frequency  to.  The  other  extreme  value,  p  "0,  characterizes  complete 
absence  of  correlations  (complete  spatial  incoherence). 

On  substituting  for  Wq  from  Eq.  (6)  into  Eq.  (3),  we  find  that 

,Sj  (P,u>)-SQ(w)l\/Rf  +  l//?f  +  [p(?(w)e‘*(R)"fi',//?,P2  +  c.c.l},  (7) 


PACS  numbers  42  68  Hf,  07.65  -b,  42  10  Mg 

I  showed  not  long  ago  that  the  spectrum  of  light  pro¬ 
duced  by  a  fluctuating  source  depends  not  only  on  the 
source  spectrum  but  also  on  the  correlation  that  may  ex¬ 
ist  between  the  source  fluctuations  at  different  points 
within  the  domain  occupied  by  the  source. 1  This  result 
was  recently  confirmed  experimentally,2  (  also  showed 
that  under  certain  circumstances  source  correlations  may 
produce  red  shifts  or  blue  shifts  of  spectral  lines  in  the 
emitted  radiation.3,4  This  prediction  has  obviously  im¬ 
portant  implications,  particularly  for  astronomy,  and  it  is 
therefore  desirable  to  verify  it  also  by  experiment. 

In  this  Letter  I  analyze  theoretically  one  of  the  sim¬ 
plest  systems  that  will  generate  spectral  shifts  by  this 
mechanism;  namely,  two  small  correlated  sources,  with 
identical  spectra  consisting  of  a  single  line  of  Gaussian 
profile.  I  show  that  with  an  appropriate  choice  of  the 
correlation,  the  spectrum  of  the  emitted  radiation  will 
also  consist  of  a  single  line  with  a  Gaussian  profile;  how¬ 
ever,  this  emitted  line  will  be  red  shifted  or  blue  shifted 
with  respect  to  the  spectral  line  that  would  be  produced 
if  the  sources  were  uncorrelated,  the  nature  of  the  shift 
depending  on  the  choice  of  one  of  the  parameters  that 
specifies  the  exact  form  of  the  correlation  coefficient. 
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FIG.  1.  Geometry  and  notation  relating  to  the  determina¬ 
tion  of  the  spectrum  Su(P.a>)  of  the  field  at  P  produced  by  two 
small  sources  with  identical  spectra  Sq(o>)  located  at  P\  and 

Pi 

where  I  have  omitted  the  arguments  P\  and  Pi  in  p.  For 
the  sake  of  simplicity,  let  us  choose  the  field  point  P  to 
lie  on  the  perpendicular  bisector  of  the  line  joining  Pi 
and  Pi.  Then  R\—Ri  (-R,  say)  and  formula  (7) 
reduces  to 

Su(P,<o)-(2/R2)SQ((o)ll  +  Rc/ip(<u)l,  (8) 

where  Re  denotes  the  real  part. 

We  note  in  passing  that  when  either  jip(a>)=0  (mutu¬ 
ally  completely  uncorrclatcd  sources)  or  when  /jp(<u) Si 
(mutually  completely  correlated  sources),  the  spectrum 
SviP,m)  of  the  field  at  the  point  P  will  be  proportional 
to  the  spectrum  Sp(oi)  of  the  source  fluctuations.  How¬ 
ever,  in  general  this  will  not  be  the  case.  In  fact,  it  is 
clear  from  formula  (8)  that  the  field  spectrum  may  differ 
drastically  from  the  source  spectrum,  the  difference  de¬ 
pending  on  the  behavior  of  the  correlation  coefficient 
Hq(io)  as  a  function  of  frequency. 

Suppose  now  that  the  spectrum  of  each  of  the  two 
sources  consists  of  a  single  line  of  the  same  Gaussian 
profile, 


FIG  2.  Red  shifts  and  blue  shifts  of  spectral  lines  as  pre¬ 
dicted  by  formula  (12).  The  spectrum  Sq( a>)  of  each  of  the 
two  source  distributions  is  a  line  with  a  Gaussian  profile  given 
by  Eq.  (9)  with  A- 1,  n*,“4  32201  x  I0‘5  sec'1  (Hg  line 
X  —4358.33  A),  *“5x10*  sec'1,  (a)  The  field  spectrum 
Su(P,a>)  at  P  when  the  two  sources  are  uncorrelated  (pp=0). 
(b),(c)  The  field  spectra  at  P  when  the  two  sources  are  corre¬ 
lated  in  accordance  with  Eq.  (10),  with  a  “1.8,  61  “7.5x10* 
sec'1,  and  (b)  u>\  “«do~2*  (red-shifted  line),  (c)  <ui“a)o 
+2*  (blue-shifted  line). 


By  straightforward  calculation  one  can  show  that  this 
expression  may  be  rewritten  in  the  form 

SL/fP,a>)-A'e~{''-mi‘)'ni?, 

where 

A’  'm(2Aa/R7)e 

a>6“  (5fa>0+^wi  )/(*?  +  &}), 
and 


(12) 

(13) 

(14) 


Sp(w)  ■ A  e 


(9) 


where  A ,  too,  and  *  («w0)  are  positive  constants.  Sup¬ 
pose  further  that  the  correlation  between  the  two  sources 
is  characterized  by  the  degree  of  spectral  coherence 


pp(co)-fle-<*’-*’’,'/24f-l.  (10) 

where  a,  w,  and  S  («cO|)  are  also  positive  constants.  In 
order  that  expression  (10)  is  a  degree  of  spectral  coher¬ 
ence,  I  must  also  demand  that  a  <  2.  On  substituting 
from  Eqs.  (9)  and  (10)  into  Eq.  (8).  I  obtain  the  follow¬ 
ing  expression  for  the  spectrum  of  the  field  at  the  point 
P 


SL(P,a>) 


(II) 


\/tf-\/Si+\/Sf.  (15) 

On  the  other  hand,  were  the  two  sources  uncorrelated, 
the  correlation  coefficient  pp  would  have  zero  value  and 
we  would  then  have,  according  to  Eqs.  (8)  and  (9), 

i.Vt;(/’.a»)lUBcorr-(2/l//?2)e  (16) 

Comparison  of  Eq.  (12)  with  Eq.  (16)  shows  that  al¬ 
though  both  the  spectral  lines  have  Gaussian  profiles, 
they  differ  from  each  other.  Since  according  to  Eq.  (15) 
So<So ,  the  spectral  line  from  the  correlated  sources  is 
narrower  than  the  spectral  line  from  the  uncorrelated 
sources.  Further,  we  can  readily  deduce  from  Eq.  (14) 
that 


too  §  too 
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according  as 

(tl  |  $  too- 

(knee  if  to i  <  wo  the  spectral  line  (12)  produced  by  the 
correlated  sources  is  centered  on  a  lower  frequency  than 
the  spectral  line  (16)  from  two  uncorrelated  sources,  i.e., 
it  is  red  shifted  with  respect  to  it;  and  if  <o\  >  a>o  the 
spectral  line  (12)  is  blue  shifted  with  respect  to  the  spec¬ 
tral  line  (16)  Figure  2  illustrates  these  results  by  simple 
examples. 

The  preceding  considerations  show  clearly  the  possibil¬ 
ity  of  generating,  by  means  of  correlations  between 
source  fluctuations,  cither  red  shifts  or  blue  shifts  of  lines 
in  the  spectrum  of  radiation  emitted  by  sources  that  are 
stationary  with  respect  to  an  observer. 

1  am  obliged  to  Mr  A  Gamliel  for  carrying  out  the 
computations  relating  to  Fig.  2.  This  research  was  sup¬ 
ported  by  the  U  S.  National  Science  Foundation  and  by 
the  U  S.  Air  Force  Geophysics  Laboratory  under  Air 


Force  Office  of  Scientific  Research  Task  No.  2310G1. 
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Hilbert  space  operators  are  introduced  into  classical  wave  theory,  which  make  il  possible  to  associate  a  unique  operator  with 
the  cross-spectral  density  By  linearly  mapping  this  operator  onto  an  associated  phase-space  one  obtains  a  wide  class  of  generalired 
radiance  functions,  including  two  well-known  ones  that  were  introduced  by  Walther  in  a  different  manner.  When  the  source  is 
quasi-homogeneous  and  the  wavelength  is  short  enough  all  these  functions  become  identical,  and  this  unique  limit  is  found  to 
have  all  the  properties  of  the  traditional  radiance,  at  least  in  the  source  plane. 
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Hilbert  space  operators  arc  introduced  into  classical  wave  theory,  which  make  it  possible  to  associate  a  unique  operator  with 
the  cross-spectral  density.  By  linearly  mapping  this  operator  onto  an  associated  phase-space  one  obtains  a  w  ide  class  of  generalized 
radiance  functions,  including  two  well-known  ones  that  were  introduced  by  Walthcr  in  a  different  manner  When  the  source  is 
quasi-homogcneous  and  the  wavelength  is  shori  enough  all  these  functions  become  identical,  and  lliis  unique  limn  is  found  to 
have  all  the  properties  of  the  traditional  radiance,  at  least  in  the  source  plane 


I.  Introduction 

In  order  to  clarify  the  foundations  of  radiometry  a 
number  of  authors  proposed  various  expressions  for 
the  radiance,  in  terms  of  the  cross-spectral  density  of 
the  light  distribution  across  the  source  [  1-5].  Unfor¬ 
tunately  none  of  them  satisfies  all  the  postulates  of 
radiometry  for  any  state  of  coherence  of  the  light  and 
it  is  now  know  that  none  in  fact  exists,  if  the  radiance 
is  to  be  linearly  related  to  the  cross-spectral  density 
|b).  Very  recently  it  was  shown,  however,  that  if  the 
source  is  quast-homogeneous,  the  expressions  for  the 
radiance  proposed  by  Walthcr  [4,1  J  acquire,  in  the 
limit  of  short  wavelengths,  all  the  properties  that  one 
postulates  for  it  in  traditional  radiometry  ( 7,8  ] . 

It  was  noted  |‘>]  that  in  its  mathematical  structure 
radiometry  has  much  in  common  with  the  phase- 
space  representation  of  quantum  mechanics.  In  par¬ 
ticular  the  phase-space  representation  of  quantum 
mechanics  deals  with  functions  which  arc  c-number 
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representatives  of  pairs  of  conjugate  operators  and 
the  radiance  is  a  function  of  a  pair  of  variables  that 
are  conjugate  in  the  sense  of  Fourier  theory  .  In  ihe 
present  paper  we  investigate  this  similarity  further 
and  we  show  that  n  leads  to  a  clarification  of  the  true 
significance  of  the  radiance. 

We  show  in  sec.  2  that  even  within  the  framework 
of  classical  wave  theory  one  may  introduce  non¬ 
commuting  operators  :i  p  and  s  which  arc  associ¬ 
ated  with  position  {p)  and  direction  (s  )  respec¬ 
tively.  In  sec.  3  we  associate  a  unique  Hilbert  space 
operator  (i  =  (/(/>.  f,  )  with  the  cross-spectral  den¬ 
sity.  Using  lliis  fact  we  then  briefly  indicate  how  a 
whole  class  of  generalized  radiance  functions  »?,,(/>. 
s)  may  be  introduced  by  linearly  mapping  <i{p.  v  ) 
onto  an  associated  p,s  -phase  space.  In  sec.  4  we  give 
explicit  expression  for  such  mappings  and  we  illus¬ 
trate  the  results  by  showing  that  the  two  expressions 
for  radiance  proposed  by  Walthcr  ate  (apart  from  a 
trivial  factor  I  just  the  phase  space  representatives  of 
)  obtained  according  to  the  so-called  Weyl 
rule  and  the  antistandard  rule  of  mapping  (operator 
ordering).  In  the  concluding  section  (see.  5)  we  show 
that  in  the  short  wavelength  limit  all  the  generalized 
radiance  functions  become  identical.  Finally  we  show 
that  when  the  source  is  quasi-honiogcnoous.  this 
unique  limit  is  a  function  that  has  all  the  properties 

We  denote  operators  with  a  caret 
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postulated  for  the  radiance  in  traditional  radionte- 
try.  at  least  in  the  source  plane. 


2.  Position  operator,  its  conjugate  operator  and 
commutation  relations  for  classical  wave  theory 

Consider  a  steady-state  (i.e.  statistically  station¬ 
ary).  planar,  secondary  source,  occupy  ing  a  portion 
n  of  the  plane  r  _  0  and  radiating  into  the  half-space 
r  »  0.  The  radiant  intensity  Jjs).  i  e  the  rate  at  which 
the  source  radiates  energy  at  frequency  r  per  unit 
solid  angle  around  a  direction  specified  In  a  real  unit 
vectors.  is  known  lobe  given  by  the  expression  (10) 

./,.(*)  (  2;rkv  t  II  (  As  .  ks  ,i').  (2.1) 

Mere 

lf'(/i ./'.  i' )  (  2ir)  J  |  |  H  (p, .  p.v ) 


•  exp(  »(/,  p\  +fp-)\  d  >,  d  >  (2.2) 

is  (he  four-diniensionai  spatial  Fourier  transform  of 
the  cross-spectral  density  function  II  (/>,.  p  ,  v). 
Further  p,  and  p<  are  position  vectors  of  any  tsvo 
source  points  S,  and  S-.  s  *(»,.  v,.  0)  is  the  com¬ 
ponent  (considered  as  a  two-dimensional  vector)  of 
j  parallel  to  the  source  plane  and  i-  is  the  component 
of  j  along  the  normal  to  the  source  plane  (  see  fig.  I  ). 
In  the  domain  of  physical  ortics.  the  radiant  inten¬ 


sity  J,.(s)  is  the  mam  measurable  quantity  relating  lo 
radiation  generated  by  the  source.  In  radiometry  the 
chief  quantity  is  the  radiance  li,(p.s )  which  is 
regarded  as  representing  the  rate  at  which  energy  at 
frequency  r  is  radiated  from  a  source  element  of  unit 
area  at  p  into  a  unit  solid  angle  around  the  s-direc- 
tion.  A  basic  radiometric  formula,  which  is  always 
introduced  from  intuitive  geometrical  considera¬ 
tions.  expresses  the  radiant  intensity  in  terms  of  the 
radiance  as 

J,,(s)=s,  [  «,.(/>.  s )  d>  .  (2.3) 

•1 

In  attempts  to  clarify  the  relation  between  radi¬ 
ometry  and  classical  wave  theory  various  expres¬ 
sions  have  been  proposed  for  the  radiance  in  terms 
of  the  cross-spectral  density  ( (  1 .5 1).  For  reasons  that 
will  become  apparent  short;,  we  will  refer  to  the 
quantities  introduced  in  ibis  way  as  wuvahu’d radi¬ 
ance  functions  and  denote  them  by  the  symbol  rf,  ( p 
s)  When  appiopnatc  we  will  attach  a  superscript  to 
this  sy  mbol  to  distinguish  between  the  different  def¬ 
initions.  As  already  mentioned  in  the  introduction, 
the  various  expressions  proposed  for  the  radiance 
have  very  similar  mathematical  structure  as  some  of 
the  phase-space  representatives  of  quantum- 
mechanical  operators.  This  similarity  suggests  that 
the  generalized  radiance  functions  may  be  just  dif¬ 
ferent  phase-space  representatives  of  one  and  the 
same  Hilbert  space  operator  We  will  show  later  that 
this  indeed  is  the  vase  llelorc  doing  so  we  will  intro¬ 
duce  a  set  ot  non-commuting  operalois  into  classical 
wave  them  y 

We  consider  the  Hilbert  space  of  square-integrable 
functions  of  p  We  associate  with  each  cartesian 
coordinate  l  and  i  ol  p  operators  \  and  f  whose 
eigenstates  |  v  >  and  i  r  ■  are  defined  in  the  usual  way: 

x ;  x  x  i  x  .  lie  ill'  ( 2.4 ) 

Next  we  introduce  the  vector  operator/!  -  (  v.  f)  and 
the  tensor  product  eigenstates  \p  |  i  |  r  .  It  fol¬ 
lows  from  eqs.  (2.4)  lhat/M/z  P\P 

For  our  purposes  the  variables  conjugate  to  \  and 
rand  the  associated  operators  may  he  most  naturally 
introduced  in  the  following  way.  We  consider  mon¬ 
ochromatic  vvavefields  f'(r.  <'  I  exp(  2zr  ii'Z )  in  the 
half-space  r  -0.  which  behave  as  outgoing  spherical 
waves  at  infinity  in  (hat  half-space.  I  (  r  r  I  then  sat- 
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isfies.  throughout  that  half-space,  the  Helmholtz 
equation 


( 1  ■  4-  A ' )  ('(r.  !»)=() . 

(2.5) 

where 

k~  Inv/c . 

(2.6) 

c  being  the  speed  of  light  in  vacuo.  It  is  well  known 
that  under  very  general  conditions  f'(r,  v )  may  be 
expressed  in  the  form  of  an  angular  spectrum  of  plane 
waves,  v  1/  .  |  I  I  | 

L’(r,  v)=  f  u(s  ,  v)  exp(  \ks-r)  d:  s  (2.7) 

Here  s  =  (.•>,.  .v,.  ,v.)s=  (*  .  ,v  )  is  again  a  unit  vector 
but  it  may  now  take  on  complex  values.  More  pre¬ 
cisely,  v,  and  a,  with  0  <o ,  <  x .  0«r  v,  <  -/  are  real  and 

•v  --  (1  -  v ;  —  v ; )  ‘  when  a; -t-.v;  s;  1  ,  (2.8a) 

--  i(  v;  +.vf  -  I )'  '  when  s  '  t-  v ;’  >  I  .  (  2.8b) 

Let  us  now  specialize  eq.  (  2.7  )  to  the  representa¬ 
tion  of  the  field  at  points  p-  (  v.  r.  0)  in  the  source 
plane  r-  0  and  let  us  operate  on  it  with  the  trans¬ 
verse  laplacian.  I  =?  ( ii/Hx.  H/Hy).  We  then  find  that 

—  » X  C(p.  t>)  j  s  .  i’)exp(iA-j  p )  d's ,  . 

(2.9) 

where 

A  -f.iln  --  I /A  - c/2ni'  (2.10) 

is  the  “reduced"  wavelength.  Lq.  ( 2.9 )  suggests  that 
we  associate  with  a  an  operators  by  the  formula 

s  --  - 1 X »  .  (2.11) 

Its  cartesian  components 

i v,  iMrin  (2.1?) 

together  with  the  position  operators  v  and  i.  ma\ 
readily  be  shown  to  obey  the  commutation  relations 

Ix.'.l  IX.  [1.  v',1  IX  .  (213) 

which  are  analogous  to  the  quantum  mechanical 
communitation  relations  for  position  and  momen¬ 
tum  |  12). 


:  v  \|.o!  cm  1 

3.  A  class  of  generalized  radiance  functions 

In  an  important  paper  dealing  w  ith  radiometry  and 
coherence.  Walthcr  [  1 )  introduced  the  generalized 
radiance  function 

tf,'"  ‘(p.  s)  --  ( M2n): s  l  \""[p.s  ).  (3.1) 

where 

(p.s  )  |  u’(/m  \p  .p  ip  .i’) 

>  exp(  i A.*  •/>  )  d  p  .  (  3.2  ) 

In  these  formulas  s  denotes  a  real  unit  vector  and 
H  (p,.p:.  r)  represents  the  cross-spectral  density  in 
the  source  plane. 

In  the  phase-space  representation  of  quantum 
mechanics  a  function  /•  '"  1  of  the  form  giv  en  by  the 
integral  (  3.2  )  is  the  so-called  Wigner  representative 
of  an  operator  U  that  depends  on  a  pair  of  non-com¬ 
muting  variables  (13).  In  analogy  with  the  relation 
between  the  Wigner  distribution  function  '  and 
the  operator  (>  which  it  represents  in  phase  space,  we 
will  rewrite  eq  (  3  7)  in  the  Idrnt  ’ 

A!"  (p. .»  I  j  P  1  \p  \  hip  s  )i p  > 

xexp(  iAj  -p  )  d  /)  .  (3  3) 

where  ' 

(Pi  |  (Up.  J,  )  Ip:  >  -  H  (p,  .  p  . .  IM  .  (34) 

The  function  ’  (p.  s  ).  defined  by  these  two 
equations,  may  be  said  to  be  the  \\  igner  representa¬ 
tive  of  the  Hilbert  space  operator  (i  <i{p.s  )  fq 
( 3.4 )  shows  that  the  matrix  elements  of  this  operator 
are  just  the  appropriate  values  of  the  cross-special 
density  H  ( p , .  p..  e).  It  seems  worthwhile  to  stress 
that  in  spite  of  its  close  resemblance  to  the  phase 
space  representation  of  quantum  mechanics  the 
above  representation  is  based  cnltivlv  on  classical 
theory . 

fheopeiator  h  hip  v  )  that  we  nisi  inliodueed 

S  iiv.v  1  Mr  pi  t'vnl  I  tic  oj  \  i>  hjsrd '  hi  i  l.ixMv .(!  r  atbri  i/i.in  iju.r' 
lum  lhoor> .  A  I  \  miIki  I  ft  « fi  h  2r, .  h  bony  I’I.hk k  s  iciv 
sum  )  «i  p  pears  in  a]  MU 

\  In-  operator  <»(/>  \  »  and  also  tin’  opcratoi  Up  /» 
s  \  »  and  tin*  lum.  tion  iJ<  u  r)  introduced  W*low  depend  on 
f.  but  we  do  not  disj>!a>  this  dependence 
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v  ia  eq.  ( 3.4 )  may  be  expressed  in  many  different  ways 
by  the  use  of  the  commutation  relations  (2.13).  More 
specifically,  if  O'  is  expressed  in  the  form  of  a  power 
scries  in  the  cartesian  components  of  p  and  i , .  each 
term  involving  a  product  of  these  elementary  opera¬ 
tors  may  be  arranged  according  to  some  chosen  rule 
of  ordering  (cf.  [14]).  One  can  associate  with  each 
such  ordered  series  a  c-number  representative  ( phase- 
space  representative)  of  O'.  The  function 
defined  by  eq.  ( 3.3).  is  among  the  best  known  rep¬ 
resentations  of  the  operator.  It  may  be  shown  to  be 
associated  with  O'  via  the  so-called  Weyl  rule  of 
ordering. 

It  is  evident  that  other  generalized  radiance  func¬ 
tions  can  be  introduced  via  formulas  of  the  form 
(3.1)  and  (  3.2  ).  with  /•  [2s  1  (/>.  s  )  replaced  by  other 
phase-space  representatives  of  the  operator  O'.  If  we 
label  the  different  representatives  of  O'  by  superscript 
i>  we  will  then  have  in  place  of  eq.  (3.1)  the  formula 

A.’ip.s )  (X/2,7)'.v  l'',.‘J‘(p.  i  I  .  (3.5) 

We  will  consider  phase-space  representatives 
F\r'{p.  s  )  produced  by  mappings  of  the  class 
investigated  in  ref.  [14).  Each  {p.  s  )  is  then 
related  linearly  lo  (Up.  ,f  )  and  it  follows  from  eqs. 

I  3.4 )  and  (3.5)  that  the  associated  generalized  radi¬ 
ance  function 

n\  -!  (p.  j)  - :  if  ■</»,  p:.  i>) ; .  <  3.6) 

where  /.'  ''  denotes  a  linear  transformation. 

We  will  impose  on  .A1,.1'1  the  constraint  that 

-/,-<  s  -.v  |  *'“'{p.  s)  d'p  .  (3.7) 

where  ./,(.*)  is  the  radiant  intensity  given  by  the 
expression  (2.1)  of  physical  optics.  In  eq  (3.7)  the 
integral  extends  over  the  whole  source  plane  r-0.  It 
will,  ot  course,  reduce  to  the  radiometric  expression 
(23)  when  (/>.  v)  vanishes  for  all p-veciors  that 
speuty  points  in  the  source  plane  outside  the  region 
n  occupied  by  the  source. 


4.  Expressions  for  (he  generalized  radiance  functions 


operators  [  1 4 [  that  every  linear  mapping  (of  a  broad 
well-defined  class)  of  operator  functions  (Up.  s  ) 
onto  c-number  functions  I-].14' ip.  s  )  is  character¬ 
ized  by  a  filter  function  Q(u.  /•).  with  the  following 
properties: 

(a)  S2(u.  r)  is  an  entire  analytic  function  of  the 
four  complex  variables  u  -( u, .  ti. ).  r-(i , .  r, ), 

(b)  Q(u.i')  has  no  zeros  on  the  real  u,-.  i\-  and 

ty-axes, 

(c)  *2(0.0)  =  !. 

The  explicit  expression  for  l\u‘  in  terms  of  O'  is  :j 

/•[/■'(p.  s  )  { 2zrX)-  J  (p,UHp.s  ) 

xJ'^’ip  P,  s  .  -/  )  \P\  >  d‘/r,  .  (4.1) 

where 

i“‘'(p  j>.\  -s.)-(2r)  1  |  \  Diu.r) 


xcxp(  -  i [«•  (p -p)  +  »*• (s  -s  )j ;  d:/<  d  ’(■ , 


and 

(4.2) 

f!(U.  1’)  =  |*2(  — 1|.  —!’))  1  ■ 

(4.3) 

On  substituting  from  eq.  ( 4.2 )  into  eq.  (4.1)  and 
then  substituting  the  resulting  expression  for 
intoeq.  (3.5)  and  making  use  of  eq.  (3.4)  we  obtain, 
after  some  calculation,  the  following  expression  for 

*y''(p.s)^(2n)  *s,  j| J  *5(«,r) 

xexp[  \iu-p  ( -vs  ,  +  J  X «t)]  U’(p,  .p,  -  Xr,  v ) 

xcxp(i#7>,  )d’wd  i  d  'p,  .  (4.4) 

We  see  that  .A)}"  is  indeed  a  linear  transform  of  the 
cross-spectral  density  H'Jseeeq.  (3.6)]. 

Now  A\y'  must  also  satisfy  the  relation  (  3.7  ),  with 
the  radiant  intensity  ./.. ( s )  given  by  eq.  (2.1  ).  i.e.  it 
must  satisfy  the  relation 

)  4\.l‘'ip,  j)d-'/> -(2jr)\v,  \\  (ks . .  ks ,  .  v)  . 

(4.5) 


It  is  know  n  form  the  general  theory  of  phase-space 
representations  of  functions  on  non-commuting 


1  Manx  of  the  formulas  pertaining  to  the  mapping  theors  devel¬ 
oped  in  ref  (14]  contain  the  “reduced”  Planck's  constant  h. 
hut  thev  also  appl>  to  the  present  case  if  one  replaces  h  b>  X 
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This  requirement  places  a  certain  constraint  on  the  > 
admissible  filter  functions  £2tu.  «•).  It  is  a  straightfor¬ 
ward  matter  to  show  that  the  constraint  is 

£<0.  *•)  =  1  for  all  e.  (4.6) 

Let  us  now  consider  some  examples.  For  the  II  or/ 
rule  of  ordering  ([14],  sec.  VII)  £2  =  £2' w  1  where 

’(u.  r)  ~  I  .  (4.7) 

If  we  use  this  fact  in  the  general  formula  (4.4)  we 
find,  after  changing  one  of  the  variables  of  integra¬ 
tion  from  e  top’  =  e/\  that  t/>,.  -  n)',."  ' .  where 

■'Jl'*"(p,  J)  ( k/2n)'s,  |  Hfpcjp.p-  Ip'.i>) 

Xexp(  -  iks  ,  p  )  d’p'  .  (4.8) 

The  expression  on  the  right  is  precisely  the  first 
expression  proposed  by  Walther  [  I  ]  for  the  radiance 
function,  which  we  have  already  encountered  [eqs. 
(3.1 )  and  (3.2)  above). 

For  the  so-called  antistandard  rule  )  [  1 4 1 .  sec.  VII) 
of  ordering,  £2~.  £2'  'N’  where 

""(«.  v)  -exp(  -  iXii-r/2 )  .  (4.9) 

If  we  use  this  expression  in  the  general  formula  ( 4.4 ) 
we  find  that 

jfJ,ASI(p,  s)  -  (k/2n)2s.exp(iks .  -p) 

X  j  W  (pi,p,  i>)  exp(  -  iks ,  p, )  d‘p,  .  (4.10) 

The  expression  on  the  right  is  the  complex  form  of 
the  second  expression  proposed  by  Walther  [4|  for 
the  radiance  function. 


5.  The  short  wavelength  limit  with  quas: 
homogeneous  sources 

Although  the  procedure  outlined  in  the  previous 
sections  leads  to  a  large  class  of  generalized  radiance 
functions,  it  is  clear  from  the  remarks  made  in  the 
introduction  that  none  of  them  will  satisfy  all  the 
postulates  of  traditional  radiometry  for  sources  of  any 
state  of  coherence.  However,  as  we  will  now  show, 
our  theory  leads  to  a  very  general  result  regarding  the 
foundations  of  radiometry. 


We  have  seen  that  the  different  phase-space  repre¬ 
sentatives  /  “‘  of  (/'  and  consequently  the  various 
generalized  radiance  functions  are  associated 
with  different  rules  of  ordering  of  products  involving 
the  noncommuting  operators p  and  s  .  However,  it 
is  seen  from  eq.  (  2. 1 3  )  that  in  the  limit  as  A  >0.  these 
operators  will  commute  and  the  distinction  between 
the  different  types  of  ordering  will  then  disappear. 
Consequently  all  the  phase-space  representatives 
of  the  operator  <7  and  hence  also  all  the  generalized 
radiance  turn  lions  .it',.1''  ml!  heroine  identical  in  the 
short-wat  elength  limn.  However,  in  view  of  Friherg's 
theorem  (6)  this  limiting  expression  cannot  be 
expected  to  have  all  the  properties  attributed  to  the 
radiance  function  in  traditional  radiometry  for 
sources  of  any  state  of  coherence. 

It  was  recently  shown  (  7)  that  when  the  source  is 
quasi-homogeneous  the  generalized  radiance  func¬ 
tion  (4.10)  is  in  the  limit  as  A  •()  (more  precisely  in 
the  asymptotic  limit  as  k  I /A  >())  given  by  the 
expression 

rl.ip.s)  -  k  -  v  / ' " 1  ( p.  i’ )  ,<'""(  ks  ,  t<) 
when  pt  a 

=  0  when  pi  a.  (5.1) 

where  f""(p.  e)  represents  the  intensity  distribution 
across  the  source  and 

£“”(/.  *')  —  < 2n)  j  x""{p  .  t'lcxpl -i/-p  )  d'/' 

(5.2) 

is  the  two-dimensional  Fourier  transform  of  the 
degree  of  spectral  coherence  of  the  light  distribution 
in  the  source  plane.  The  expression  (5.1)  was  shown 
to  have  all  the  properties  attributed  to  the  radiance 
function  in  traditional  radiometry.  Ii  billows  from 
this  result  and  the  result  established  earlier  in  this 
section  (italicized  above)  that  when  the  source  is 
{jiMsi-hornopcncous.  all  the  generalized  radiance 
functions  M'fi''  (p.  s)  of  the  class  that  we  hair  consid¬ 
ered  in  tho  flop, -i  hate  the  same  asvniptoln  limn, 
given  l<y  ci/.  (  5. 1  ),  as  k  •  /  ;  and  this  ( oinmon  hom¬ 
ing  form  of  all  the  generalized  radiance  tunc: ions  may 
he  identified  with  the  radiance  of  traditional  ntdtome- 
trv.  at  least  at  all  points  In  the  source  plane. 

It  should  be  evident  that  we  have  not  proven,  in 
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the  strict  mathematical  sense  that  radiometry.  even 
for  the  restricted  class  of  quasi-homogeneous  sources, 
is  the  asymptotic  limit  for  large  wave  numbers  of  sta¬ 
tistical  wave  theory.  In  this  connection  it  might  be 
worthwhile  to  point  out  that  the  somewhat  analo¬ 
gous  statement  frequently  made  that  classical 
mechanics  is  the  limit  of  quantum  mechanics  if 
Planck’s  constant  h— 0  has  not  been  rigorously  justi¬ 
fied  to  this  day,  and  that  even  the  restricted  class  of 
systems  for  which  this  statement  may  perhaps  be  true 
has  not  been  precisely  defined.  Nevertheless  we  are 
ot  the  opinion  that  the  results  derived  in  this  note 
provide  a  genuine  insight  into  the  true  meaning  of 
the  radiance. 
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An  exilel  law  is  derived  fur  the  propagation  in  Iren  span*  u!  tin*  lirsl  gcnerall/ea  radiance  function  introduced  bv 
Wallher  A  simplified  lorm  ol  this  law  is  obtained  for  the  case  when  the  source  is  ipiasi  homogeneous  It  is  also 
shown  that  when  the  source  is  quasi  homogeneous  and  the  wave  nuniher  is  large  enough  (the  wavelength  is 
siillicientlv  short  1  (he  first  generalized  radiance  acquires  all  I  he  properties  of  the  radiance  of  traditional  radioiaetrv. 


1.  INTRODUCTION 

In  two  well-known  papers'  -  dealing  with  the  foundation  of 
rudiometry,  Walther  introduced  certain  generalized  radi¬ 
ance  functions.  These  functions  have  some  of  but  not  all  the 
properties  that  are  attributed  to  the  radiance  in  traditional 
rudiometry. 1  Much  subsequent  work,  aimed  at  clarifying 
the  connection  between  radiometry  and  physical  optics, 
made  considerable  use  of  these  functions.  We  will  refer  to 
the  generalized  radiance  function  introduced  in  Kefs.  1  and  2 
as  the  lirst  and  the  second  generalized  radiance  functions 
(g.r.f.'s),  respectively. 

Approximate  transport  equations  for  the  propagation  of 
either  of  these  two  g.r.f.'s  were  obtained  by  Walther.4  .)ann 
son,  '  Friberg,"  Pedersen,'  and  Bastiaans.8  An  exact  law  for 
the  propagation  of  the  second  g.r.f.  in  free  space  was  recently 
obtained  by  Foley  and  Wolf,"  who  also  showed  that  when  the 
source  uf  the  optical  field  is  quasi-homogeneous  this  func¬ 
tion  acquires,  in  the  short -wavelength  limit,  all  the  proper¬ 
ties  that  are  postulated  for  the  radiance  in  traditional  radi¬ 
ometry.  They  also  obtained  an  explicit  expression  for  this 
limiting  form  of  the  second  g.r.f.  in  terms  of  the  distribution 
of  the  intensity  and  of  the  degree  of  spectral  coherence  of 
light  across  the  sourer. 

In  the  present  paper  we  derive,  to  begin  with,  an  exact  law 
lor  propagation  ol  the  first  g.r.f.  in  free  space.  We  then 
consider  the  lortn  that  this  propagation  law  takes  when  the 
source  is  quasi-homogeneous.  Finally  we  conuder  the  as 
ymptotic  limit  for  large  wave  number  (short  wavelength)  of 
the  lirst  g.r.l.  in  optical  fields  produced  by  a  quasi -homoge¬ 
neous  source,  and  we  find  that  it  is  identical  to  the  eorre 
sponding  limiting  form  obtained  for  the  second  g.r  I  in  Kef 
I).  I  Older  these  circumstances  the  propagation  law  is  found 
to  reduce  to  the  usual  radiometric  transport  equation. 
These  results,  together  with  those  derived  in  Kel.  SI.  go  a  long 
way  toward  clarifying  the  foundation  of  radiometry. 

2.  PROPAGATION  LAW  FOR  WALTHER'S 
FIRST  GENERALIZED  RADIANCE 

Let  us  consoler  a  secondary  source  nccupvmg  n  lintlv  do 
loam  -i  ill  the  plane  z  =  (I  and  radiating  into  tin-  hall  space  z 


>  0.  We  assume  that  the  field  fluctuations  in  the  source 
plane  are  characterized  l»v  a  stationary  statistical  ensemble. 

We  denote  by  Wlri,  r ..  w)  the  cross-spectral  density  of  the 
emitted  light  at  any  two  points  l’i  and  l’:l,  specified  by 


position  vectors  iq  and  r>,  in  the  half-space  z  >  0.  Let  us 
choose  the  two  points  to  tie  located  in  some  plane  z  =  con¬ 
stant  >  0.  which  we  will  denote  by  It,  and  let 

r  =  ( r ,  +  r.)/2,  p  =  r,  -  r,.  (2.1a) 

Then 

r,  =  r  +  p/2,  r ,  =  r  -  p/2.  (2.11)) 


(See  Fig.  1,  where  1’  denotes  the  point  with  position  vector 
r.)  An  expression  for  the  first  generalized  radiance  intro¬ 
duced  by  Walther,1  and  defined  by  him  at  points  r  in  the 
source  plane  z  =  (I,  can  readily  he  generalized  to  apply  to 
field  points  r  in  any  transverse  plane  11  in  the  half-space  z  > 
t).  It  takes  the  form 

ViJ r.  s)  =  ^  j  v  |  lV(r  +  p/2,  r  -  p/2)exp(i7rs  ,  ■  p)d>, 

(2.2) 

where 

/,'  =  -  “*  (2.,'t) 

c  .\ 

is  the  wave  number  associated  with  the  frequency  w  and 
wavelengt  h  i  is  the  speed  of  light  in  vacua,*  (s,,  s\,  s.. )  is 
a  real  unit  vet  tor.  and  s  (v„  s  .  (It  is  its  transverse  compo¬ 
nent  i<  oit-idered  a  two  dimensional  vector). 

\\  i  will  now  derive  an  expression  lor  /, .  in  terms  of  the 
value  t.  licit  t lie  generalized  radiance  | F.q.  l2.2)|  takes  in 
t  lie  source  plane  .  -  u  For  t  Ins  purpose  vve  will  make  use  of 
the  lollovving  result  established  not  long  ago.1'1  The  cross 
spectral  detisiiv  H  ir,.r  .  *•)  nav  lie  represented  in  terms  of 
an  appropriate  ensemble  |(  i  r.  w l|  of  monochromatic  wave 
lunetions  iwnb  lone  dependeme  e\pl  is/l  understood) 
Ilia!  propagate  Irmii  the  source  plane  -  ('  into  the  ball 
spai  e  .  »  O.  m  I  lie  lurin 

ll'ir,.  r  '  •  f  '*( r , .  o  K  <r  wl  (2  1) 
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Here  the  angle  brackets  denote  the  average  taken  over  this 
ensemble.  Now  the  value  of  l Hr,  ui)  at  any  point  r  in  the 
half-space  z  >  0  may  be  expressed  in  terms  of  its  boundary 
values  ( o>)  at  point  r„  in  the  source  plane  z  =  0  by  the 
use  of  Rayleigh’s  formula" 


where  we  used  the  fact  ii.  u  d  r.^d-r,,.  =  d-r„d  p,>.  Let  us 
now  substitute  front  Kq.  (2  In  otto  the  expression  (2.2). 
interchange  the  order  of  integrations,  and  make  use  of  the 
fact  that  the  generalized  radiance  /L""(r,,.  s)  at  points  in  the 
source  plane  is  given  by  |cf.  Kq.  (2.2)] 


l '( r.  **)!  =  j  f,’(r  -  r„.  w)r,‘"(r„,  w)d-’r„,  (2.5)  /j  ""(r,„8)  =  ^  V  »  j  VL"’’i 


(r„  +  pj 2,  r„  -  p„/2,  w) 


where  (V(r  -  r„,  *>)  is  the  Green's  function 

1  ,•)  (~ exp(/A-ir  —  r,,l)~I 

<;<r-r„.w)  =  -  ‘  ,  ,  •  (2.6) 

2tt  itz  ^  Ir  -  r„l  ] 

On  substituting  from  Kq.  (2.5)  into  Kq.  (2.4),  we  readily 
find  that 


X  exp  Ukai  •p0)dv>ii,  (2.12) 

where  s.  is  the  z  component  of  the  unit  vector  s.  We  then 
obtain  the  following  expression  for  the  generalized  radiance 
ft  Jr,  s)  at  any  point  r  in  the  half-space  z  >  0  in  terms  of  its 
boundary  values  /L"(r<,,  s)  at  points  r,>  in  the  source  plane: 


W'(rl.r,.u>)  =  (.'*(R1,u.’)<7(R,,  w) 

J  Jl.-oil 

X  W,m(r0„  r(),)dJr,)1d"r„,.  (2.7) 

where  r,,i,  r.(2  are  the  position  vectors  of  two  typical  points  Qi 
and  Q  .  in  the  source  plane, 

R  =r,  -  r,„  (7=1,2).  (2.8) 

and 

VV""(r,,,,  r„„  u:)  =  <( . ■(rl,).ui)(f""(r11,.  w>  •  (2.9) 

is  the  cross-spectral  density  of  the  light  in  the  source  plane 
The  integration  on  the  right  hand  side  of  Kq.  (2.7)  is  taken 
twice  independent Iv  over  the  source  plane. 

Let  us  now  change  the  variables  in  Kq.  (2.7)  according  to 
Ihe  transformation  (2.1b)  and  according  to  a  similar  trails 
tormation  involving  the  source  variables: 

I'm  =  r„  +  e„/2.  r„,  =  r,,  -  p„/2.  (2.10) 

i  r,,  represents  the  point  Qn  in  Fig.  1.)  The  formula  (2.7)  then 
takes  the  form 


/(u.(r,s)=  K(r  -  r„.  s;u;)-^J",(rH.  s)d'Jf„.  (2.12) 

Ai-«ii 

The  kernel  A( r  —  r«,  s:  u i)  in  this  integral  is  given  by  the 
formula 

Air  -  r„.  s;  <p)  =  fi'*(r  -  r,,  +  p'/ 2.  wK.’lr  -  r„  -  p' 12,  w) 

Jn 

X  expOTs  -p'ldV-  (2.14) 

In  deriving  this  expression  we  changed  the  variable  of  inte¬ 
gration  from  p  top'  =  p  -  p„. 

The  formula  12.12)  represents  an  exact  law  lor  the  propa 
gation  of  Wall  tier's  first  radiance  function  from  Ihe  source 
plane  z  =  0  into  the  half  space  z  >  0. 

We  will  now  specialize  the  general  formula  (2.12)  to  fields 
generated  by  quasi -homogeneous  sources  that  are  of  partic¬ 
ular  interest  in  connection  with  the  foundation  of  radiome- 
try.'1 


( ,  3.  PROPAGATION  OF  THE  GENERALIZED 

Wtr  +  p/2,  r  -  p/2,  u>)  =  G"|r  -  r„  +  (p  -  p„)/2,  w|  RADIANCE  WHEN  THE  SOURCE  IS 

J  QUASI-HOMOGENEOUS 


8  <i[r  -  r„  -  (p  -  p„)/2.  w| 

<  VV  r„  +  p,,/2.  r,-»  .  at  r,,d  p„. 


When  the  source  is  quasi -homogeneous,  the  cross  sped ral 
(2.1 1 )  densil  v  function  H'""(r|.  r .,  w>  has  the  fonn1- 
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. . .  r„  ad  =  r^r'  *  r\  -  r„  a),  (3.1) 

where  /""(r,  a)  is  the  intensity  distribution  and  g'm(r't  a)  is 
the  degree  of  spectral  coherence  of  the  light  in  the  source 
plane.  /""( r,  a)  is  assumed  to  be  a  slow  function  of  r.  whereas 
A‘""(r',  a)  is  assumed  to  he  a  fast  function  of  r'. 

On  substituting  from  Kq.  (3.1)  into  Kq.  (2.12),  we  obtain 
the  following  expression  for  the  generalized  radiance  of  a 
quasi  homogeneous  source  at  any  point  r0  in  the  source 
plane: 

s)  =  fe's7/""(r„,  wig^ika  . ,  a).  (3.2) 


g""(f,  a)  =  -  1  f  gl0l(r',  a)exp(-)f  ■  r')d-’r'  (3.3) 

(2rr)2  J 

is  the  Fourier  transform  of  A'""(r',  a).  To  determine  the 
generalized  radiance  of  the  field  at  any  point  r  in  the  half 
space  z  >  0  produced  by  the  quasi  homogeneous  source,  we 
substitute  from  Eq.  (3.2)  into  Eq.  (2.13)  and  find  that 

'A  Jr,  s)  =  k-sj)m{ks  : ,  a)Af(r,  s:  a).  (3.4) 


Mir,  s;  a)  =  Kir  -  r,„  s;  a)/""(r„,  a)d'r„ 


and  the  kernel  K  is  given  by  Eq.  (2.14). 

4.  SHORT-WAVELENGTH  LIMIT  OF  THE 
GENERALIZED  RADIANCE  OF  A  FIELD 
GENERATED  BY  A  QUASI-HOMOGENEOUS 
SOURCE 

Let  us  now  consider  the  behavior  of  A  Jr,  a)  of  a  field  gener¬ 
ated  by  a  quasi -homogeneous  source  in  the  short-wave¬ 
length  limit,  more  precisely,  in  the  asymptotic  limit  as  the 
wave  number  h  =  2ir/X  *  «>.  For  this  purpose  we  first 
express  the  (ireen’s  function  (2.6)  in  a  more  explicit  form. 
On  carrying  out  the  differentiation  we  readily  find  that 


f/(R,  a)  =  —  *  ( ik  —  \  ^  Z 

2ir  |_\  ft)  It 


which,  for  sufficiently  large  values  of  hit,  may  be  approxi¬ 
mated  by 


(i(R,  a)  ■ 


ih  /z\ 

2ir  \/?/  It 


Next  we  substitute  from  expression  (4.2)  into  the  expression 
12.  If)  for  the  propagation  kernel  At  r  -  r, ,,  s;  a)  and  find  that 
when  It i  »  1,  hit  »  I: 


Kir  -  r,„  s;  a) 


fh:Y  (  ,.'**■ 


2*/  In  H , 


expli 'ha  ,  •  p  Hi-fi', 


Let  us  next  determine  the  asymptot ic  approximation  as  h 
-  *  for  the  factor  A/(r,  s,  a),  defined  by  Kq.  (3.5).  which 
enters  the  ex()ression  (3.4 )  for  the  radiance  fund  ion  of  a  field 
generated  by  a  quasi-homogeneous  source.  For  this  purposc 
we  substitute  for  Kir  -  r„,  s;  a)  from  expression  (4.3)  into 
Kq.  (3.5)  and  introduce  the  new  variables 


One  then  readily  obtains  the  following  expression  lor  A/(r. 
s,  w): 


Mir,  a:  a) 


f  expt-Mlr  -  p,l  , 

X  dP|  i  r  ' 

I  lr  -  p , r 

X  exp(— ika  -P]). 


exp(i/fs  •  p.J 


(T‘) 


where  we  have  made  use  of  the  fact  that  dVnd-p  =  d-pid-p.  . 

We  may  express  Kq.  (4.6)  in  a  more  symmetric  form  by 
making  use  of  the  fact  that  because  (he  source  was  assumed 
to  lie  quasi  homogeneous,  /""(p.  u.  )  will  change  slowly  with  p 
for  each  el  feel  i\  e  I  requeues  w  ( Iml  eon  I  ribules  to  (he  source 
spectrum  Hence  we  mas  make  I  In  approximation 

*/•'.*)  '  |/""(pl..4|'  |r"(p,.w)|'  -  (4.7) 

on  the  right-hand  side  of  expression  (4.(4).  The  resulting 
expression  for  Af(r.  s;  w)  may  then  be  written  in  the  form 


A/(r.  s;  w) 


Fi  r.s.w)  =  (/"%,.  *.4| 


I  Fir.  s;  s’lf’lr.  s;  w). 


expU'Alr  -  p,l) 


lr  -  pj" 


X  expd/.s  ,  -  p, Id'Pe  (4.91 

The  asymptotic  approximation  to  the  integral  on  the 
right-hand  side  of  Kq.  14.9)  may  be  determined  by  the  use  of 
the  two-dimensional  form  of  the  principle  of  stationary 
phase.1,1  In  carrying  out  the  calculations  we  ignore  the  de 
pendente  of  the  source  intensity  /""(p,  w)  on  w,  for  reasons 
indicated  in  connection  with  the  approximation  (4.7).  The 
result  is 


Fir.  s;  u-) 


I/,,., |p  -  i;/s.)a  .  sell*  '  exp(t/rs  •  rl 


when  X,  .  » 


R,  =  r  -  r„  +  p7 2. 


'  i)  when  .X,  a  (4.  Id) 

as  k  -  '»■.  where  .V,  is  the  point  in  the  source  plane  4  =  0 
whose  position  vector  is  p  -  (c/s.)s  (see  Fig.  2).  On  stihsli 
tut  ing  from  expressii  m  (  ! .  ID)  into  expression  (IS),  we  obtain 
lor  A/(r,  »;  cl  (he  asymptotic  approximation 

A/(r,  s:  u  )  -  /""Ip  -  (.'/>■  is  .  J  when  S  ,  * 

,  (till 

0  when  \, •  • 


R .  =  r  -  r„  -  p/2. 


(4.4b) 
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Kin-  -  Illustrating  the  notation  relating  to  the  formula  (4.14).  So 
is  the  point  in  the  source  plane  at  which  the  line  through  the  field 
point  I’  in  the  direction  of  the  unit  vector  s  intersects  that  plane. 
The  vector  p  -  I'/s.ls  ,  which  appears  in  expressions  (4  10),  (4. 1 1 ), 
and  14.1.3).  is  the  position  vector  of  the  point  S„. 

On  substituting  from  expression  (4.1 1)  into  Kq.  (.‘1.41,  we 
see  at  once  that 

fijr,  s)  ~  Bjr,  s)  as  ft  »  *>,  (4.12) 

where 

/(Jr,  s)  = /r’.s./''"|p  -  (z/.Os  ,  .  w|  when  ,S(I  (  a 
=  (I  when  .S’,,  s  n 

(4.13) 

Since  p  -  Iz/.vjs  ,  is  the  position  vector  of  the  point  S(l  at 
which  the  line  through  the  field  point  P  in  the  s  direction 
intersects  the  source  plune.  Kq.  (1.13)  may  he  rewritten  in 
i  he  form 

/<</’.  s)  =  /r.s  ./""(.s’„,  ud£,,"(*s  ,  a)  when  s  ►  SI,. 

=  I)  when  s  0  !/,. 

(4.14) 

where  S’, •  is  the  solid  angle  generated  by  the  lines  from  all  the 
source  points  to  P  (see  Kig.  2). 

I'he  expression  (4.14)  is  identical  (except  for  a  slight 
i  h.inge  in  notation)  with  the  expression  derived  in  Kef.  9  for 
the  asymptotic  limit  as  k  »  <®  of  the  second  generalized 
radiance  function  of  Walther,  under  the  assumption,  made 
also  in  the  present  paper,  that  the  source  of  the  field  is  a 
quasi  homogeneous  source.  It  has  also  been  shown  in  Ref.  9 
that  the  expression  14.14)  satisfies  all  I  he  postulates  of  trndi- 
. . .  radiometry  in  free  space.  We  may,  therefore,  con¬ 
clude  by  saying  that  the  analysis  presented  in  this  paper 
supports  the  view  that  traditional  radiometry  may  be  re¬ 
garded  as  l  he  asymptotic  limit  lor  large  wave  numbers  I  short 
wavelengths)  of  statistical  wave  theorv  ol  fields  produced  by 
quasi  homogeneous  sources. 
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Not?  addl’d  in  proof  Since  this  paper  was  written  anoth¬ 
er  paper  dealing  with  the  foundations  of  radiometry  was 
published  [M.  Nieto-Vesperinas,  “Classical  radiometry  and 
radiative  transfer  theory:  a  short -wavelength  limit  of  a  gen¬ 
eral  mapping  of  cross-sped  ral  densities  in  second-order  co¬ 
herence  theorv,"  3.  Opt.  Soc.  Am.  A  3,  1354-13fi9  (1986)). 
Unfortunately  the  main  conclusions  of  that  paper  are  incor¬ 
rect  because  the  analysis  contains  several  errors.  Specifical¬ 
ly,  the  expression  (29)  of  that  paper  is  not  the  onlv  expres¬ 
sion  that  satisfies  Kq.  (27).  Moreover,  the  flux  equation  (16) 
does  not  imply  that  Hr,  s)  is  positive  definite,  as  is  stated 
below  Kq.  (32). 
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The  .S  matrix  is  first  introduced  within  the  framework  of  the  angular  spectrum  representation  of  wave  fields 
interacting  with  linear  dielectric  bodies  of  arbitrary  shape.  By  using  some  universal  properties  ol  the  .V  matrix,  a 
number  of  relations  involving  certain  generalized  reflection  and  transmission  coefficients  are  derived.  '1  hese 
relations  may  be  regarded  as  generalizations  of  two  well-known  classic  reciprocity  relations  due  to  <!  (1  Stokes. 
Two  reciprocity  relations  involving  the  reflection  and  the  transmission  coefficients  for  interaction  ol  a  plane 
electromagnetic  wave  with  a  stratified  dielectric  medium  are  obtained  as  special  cases. 


1.  INTRODUCTION 

In  a  classic  paper  published  in  1849,  Stokes1  derived  two 
well  known  reciprocity  relations  involving  reflection  and 
transmission  of  light.  More  specifically,  he  considered  a 
plane  monochromatic  wave  incident  upon  a  plane  boundary 
separating  two  semi-infinite,  homogeneous,  isotropic  dielec¬ 
tric  media.  Suppose  that H,  and  M,  are  the  angles  of  incidence 
and  refraction,  respectively,  when  the  wave  propagates  from 
the  first  into  the  second  medium,  and  that  r  and  t  are  the 
corresponding  reflection  and  transmission  coefficients. 
Next  consider  the  situation  when  the  wave  is  incident  at  an 
angle  »/  =  II,  from  the  second  into  the  first  medium,  and  let  /i 
and  f  he  the  corresponding  relied ion  and  transmission  coef¬ 
ficients.  The  relations  derived  by  Stokes  are 

rt  +  r*  =  I,  (1.1a) 

p  +  r  =  0.  (1.1b) 

The  relations  (1.1)  were  later  generalized  to  somewhat 
more  complicated  situations  involving  stratified  media,  and 
they  have  played  a  useful  role  in  optics  of  thin  homogeneous 
films.-  More  recently,  relations  of  this  kind  have  li  -come  of 
importance  in  some  investigations  concerning  the  cancella¬ 
tion  of  distortions  by  the  technique  of  phase  conjugation.*-4 
All  these  situations  have  one  feature  in  common.  They 
involve  a  homogeneous  or  a  succession  of  homogeneous  di 
elect rics  with  mutually  parallel  planar  boundaries,  and.  con¬ 
sequently.  when  a  plane  wave  is  incident  upon  such  a  config¬ 
uration  only  one  reflected  and  one  transmitted  wave  is  gen¬ 
erated.  It  seems  natural  to  inquire  whether  one  can 
generalize  the  Stokes  relations  further,  so  that  they  apply  to 
situations  such  as  rough-surface  scattering  and  scattering 
from  an  inhomogeneous  plane-parallel  dielectric  slab  or  to 
phase  conjugation  of  waves  that  are  scattered  from  a  dielec- 
t  ric  body  of  arbitrary  shape.  In  the  present  paper  we  obtain 
sin  h  a  generalization  within  the  framework  of  the  scalar 
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wave  theory.  Our  derivation  utilizes  in  a  basic  way  the 
concepts  of  the  angular  spectrum  representation  of  wave 
fields  and  of  the  S  matrix.  The  combined  use  of  these  two 
concepts  has  already  proved  rather  useful  in  treatments  of 
other  problems,  which  yielded  interesting  results  relating  to 
the  theory  of  antennas  and  to  distortion  correction  by  phase 
conjugation.11  The  generalization  of  the  Stokes  relations 
presented  in  this  paper  does  not,  however,  appear  to  have 
been  obtained  previously. 

2.  SOME  GENERAL  RELATIONS  INVOLVING 
THE  ANGULAR  SPECTRUM 
REPRESENTATION  OF  WAVE  FIELDS  AND 
THE  S  MATRIX 

Consider  a  monochromatic  field,  not  necessarily  a  planar 
one,  incident  upon  a  dielectric  scatterer.  We  denote  by 
(/"’(r)exp(- iwf)  and  ( -'""(r)exp(-iu;/ )  the  incident  and  the 
scattered  fields,  respectively,  with  r  denoting  the  position 
vector  of  a  typical  point  in  space,  t  the  time,  and  w  the 
frequency.  The  total  field  ( '(r)exp(-iuit)  is.  of  course,  the 
sum  of  these  two  fields. 

I,el  us  choose  a  Cartesian-coordinate  system  of  axis  so 
that  the  scatterer  is  situated  within  Ihe  strip  0  <  z  <  L,  and 
let  /?  and  ft '  lie  Ihe  two  half-spaces  on  the  two  sides  of  the 
scatterer  (see  Fig.  II  It  is  well  known  that  under  very 
general  conditions  the  total  field  in  each  ot  the  two  hall 
spaces  may  be  represented  in  the  form  of  an  angular  spec¬ 
trum  of  plane  waves,  both  homogeneous  and  evanescent 
ones  -  I  he  amplitudes  of  the  evanescent  waves  decay  expo¬ 
nentially  with  increasing  distance  from  the  scatterer.  Be¬ 
cause  we  will  be  interested  only  in  the  field  faraway  from  the 
scatterer.  we  will  omit  the  contributions  of  the  evanescent 
waves.  The  angular  spectrum  representation  of  the  time 
independent  part  of  the  total  field  then  takes  the  following 
form  tn  (he  two  halt  spaces 
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z=0  z'=  L 

Fig.  i.  Illustrating  the  notation. 


In  /r: 

(’( r)  =  -'*|  <  ’  '(n)c'*n'dU  +  ‘k  I  />''(n)e'*n  rdU, 

2rr  2  it  !„  • 


incoming  at  infinity,  whereas  the  integrals  containing  the 
spectral  amplitudes  /)' 1 1  represent  a  field  that  is  outgoing  at 
infinity.'1 

We  will  assume  that  the  scatterer  is  linear,  i.e..  that  the 
outgoing  field  depends  linearly  on  the  incoming  field.  Con¬ 
sequently  the  C  amplitudes  and  the  I)  amplitudes  will  he 
coupled  by  a  relation  of  the  form1" 

D(n)  =  -  f  Stn,  n'lCtnldU',  (2.4) 


where  C  and  D  are  the  column  vectors 


D(nl 


//  +  l(n) 
/)■-’(  n) 


(2.a) 


and  Sis.  /or  each  pair  of  arguments  n  and  n'.  a  2  X  2  matrix. 
Written  out  more  explicitly,  Kq.  (2.4)  gives 


In  /C: 


42.1a) 


'  ’( rl 


,k  [  C,  +  ,(n)e‘*nrdS2+  (  //*'(n)c'*“  rdl>. 

2ir  J„i  - 1  2 ir 

(2.1b) 


In  these  formulas  n  -  («,.  n,,  n,)  are  real  unit  vectors.  A  =  m/c 
is  the  wave  number  associated  with  the  frequency  u>  (c  being 
the  speed  of  light  in  eavuo).  dll  is  the  element  of  solid  angle 
generated  by  the  unit  vector  n,  and  o'"  and  it1"'  are  unit 
hemispheres  in  n  space  defined  as 

a'*1:  rr  =  1,  n,  >  0,  (2.2a) 

a''1:  rr  =  1 ,  n.  <  0.  (2.2b) 


1  he  constants  multiplying  the  integrals  in  Kqs.  (2.1 1  have 
been  chosen  so  as  to  simplify  subsequent  formulas. 

In  the  representation  (2.1)  the  factors  and  /)'"  have 
the  physical  significance  of  (generally  complex)  amplitudes 
of  homogeneous  plane  waves  that  propagate  in  different 
directions  either  toward  the  scatterer  [waves  with  ampli¬ 
tudes  f'  * '  and  P  "'I  or  away  from  it  |waves  with  amplitudes 
i  "  and  !)'" 1.  However,  they  also  have  another  physical 
significance,  which  becomes  evident  when  one  examines  the 
behavior  of  the  total  field  far  away  from  t  he  scatterer.  One 
'  ben  finds,  for  example  by  the  use  of  the  principle  of  station¬ 
ary  phase.’ "  that  as  the  distance  r  of  the  field  point  from  the 
fixed  origin  0  in  (he  source  region  increases  along  any  fixed 
direction  specified  by  any  real  unit  vector  u  (u,,  uv,  u„), 


(  (rul 


C‘'(~u) 


+  I)' 


as  hr  *•  «>, 


(2.2) 


whore  the  upper  or  the  lower  signs  are  taken  on  the  right- 
hand  side  according  to  whether  the  field  point  r  =  ru  is 
located  in  the  half -space  ~H '  or  H  .i.e..  according  to  whether 
n.  >  (lorn.  <0. 

The  formula  (2.2)  expresses  the  far  field  in  each  of  the  two 
(tall  spaces  on  either  side  ot  the  scatterer  as  a  sum  of  a 
i  'merging  and  a  diverging  spherical  wave,  with  complex 
impbindes  C"  and  /)'''  (see  Fig.  2)  This  result  implies 
lh.il  (he  integrals  in  Kqs.  (2  1)  that  contain  the  (generally 
"inplexl  spectral  amplitudes  :>  presenl  a  field  that  is 


l)'*‘(nl  =  -  '(n,  n'lf"'  '(n'ldlF 


>  '  tn,  n'l  r'^'tn'ldt!',  (2.6al 


l)'  ’(n)  = 


.S"  'in.  n')C'" '(n')dl! 


with 


I, 


S""  (n,  n')  (" 4  '(n  idi!',  (2.6b) 


„  „  r.S"*-'(n,  n')  .S' f  "(n,  n'l 

S(n, n  )  = 

L.SM--'(n,  n')  .S""-"(n.  n  ) 


(2.7) 


From  the  significance  of  the  quantities  C 11  and  D 1,1  as 
complex  amplitudes  of  waves  that  propagate  either  toward 
or  away  from  the  scatterer,  and  recalling  the  definitions  (2.21 
of  the  domains  of  integration  in  Kqs.  (2.6),  it  is  clear  that  the 
four  elements  of  the  2x2  matrix  (2.2)  are  defined  only  for 
the  following  ranges  of  the  z  components  of  the  unit  vectors 
n  and  n’: 


S"  '(n,  n'):  n.  >  0.  n.' >  (), 


(2.8a) 


.s"  ’  ‘ 1 1  n ,  n  ):  n .  >  0.  n .'  <  0, 


(2.8b) 


c:  i-u) 


(-u)- 


uJn>  u 


( 1  (u) 


(u' 


t? 


SCAT TERFR 


ft 


ft' 


Fiji-  2  The  far  fields  in  (he  half  spa<  <*>  t\"  and  fi  on  either  side  ol 
the  soitterrr 
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S'  ■'(n.nl:  n.  <  (I.  n/>  0,  (2.8c) 

n  ):  n.  <  0,  n.'<0.  <2.8d) 

It  will  he  useful  at  this  point  to  make  a  contact  with  the 
usual  .S'  matrix  of  potential  scattering.  In  the  theory  of 
potential  scattering  the  expression  for  the  far  field,  which 
corresponds  to  expression  (2.  .'11,  would  he  written  in  the  form 

,,  O' 

I  Iru)  -*•  7\(u)  +/■'.,( u)  •  as  hr  ••  <»,  (2.9) 

r  ‘  r 

and  one  expresses  the  relution  between  the  complex  ampli¬ 
tudes  of  the  outgoing  and  incoming  waves  in  the  following 
form,  which  corresponds  to  Kip  (2.4): 

A',(n)  =  -  j  A'(n,n')F4-n')dsr,  (2.10) 

where  the  integration  extends  over  the  whole  unit  sphere 
generated  by  the  unit  vector  n'  |cf.  Ref.  11,  Eqs.  (20)  and 
(21 1 1 . 

It  is  clear  that  A  (n.  n)  and  JS(n,  n')  are  essentially  the 
same  <piant it ies,  hoth  being  continuous  matrices  whose  ele¬ 
ments  are  labeled  hy  pairs  of  real  unit  vectors  n  and  n'. 
However,  in  contrast  with  -V(n,  n'l,  each  element  of  l$(n,  n') 
is  itself  a  matrix,  arising  from  the  partition  of  A  (n.  n')  into 
four  separate  contributions  |K<|s.  (2.7)|.  Such  partition  is 
advantageous  when  the  field  in  each  of  the  two  half-spaces 
K  and  H  *  is  represented  in  the  form  of  an  angidar  spec¬ 
trum  of  plane  waves. 

For  later  purposes  we  recall  some  generul  properties  of  the 
usual  scattering  matrix  t  tn,  n't  It  is  well  known  that  when 
the  scatterer  is  dielectric  (i.e.,  lossless),  A'  is  unitary,  i.e.,  it 
obeys  the  relations  |cf.  Ref.  1 1,  Kqs.  (24)  and  (29)) 

J  A  *tn.  n')A(n,  n")dil  =  A(n'  -  n").  (2.11a) 


A'(o',  n)A*(n",  n)d(!  =  A(n'  -  n"),  (2.11b) 


where  the  asterisk  denotes  the  complex  conjugate  and  the 
integrations  extend  over  the  unit  sphere  generated  hy  the 
unit  vector  n  Further.  A(n'  -  n")  is  the  "spherical”  Dirac 
delta  function,  defined  hy  the  formula 


Ain'  -  n  ") 


AIM'  -  rt" )A(sr'  -  *-"> 
Isin  M'| 


(2.121 


where  (M’,  y')  and  (M",  *-")  are  the  spherical  polar  oordinates 
of  the  unit  vectors  n'  and  n",  respectively,  and  A  is  the  usual 
one  dimensional  Dirac  della  function. 

The  A  matrix  also  obeys  the  reciprocity  relation  |Ref.  1 1, 

K<|  128 1  [ 


A(-n'.  -  n)  =  Ain.  n  ).  12.1.1) 


We  show  in  Appendix  H  that  when  the  incident  field  is  a 
plane  wave. 

f  (n„  -  1).  (2.14) 


the  (actors  F ,  and  F  in  the  asymptotic  approximation  (2.9) 
ot  the  total  field  (incident  t  scattered)  are  given  by 
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/•j(n)  =  A(n,  n„).  (2.15a) 

/•' .( n )  =  -  A(n  +  n„l.  (2.15h) 

tk 

3.  THE  GENERALIZED  TRANSMISSION  AND 
REFLECTION  COEFFICIENTS 

We  will  now  show  that  the  four  elements  of  the  partitioned  .S’ 
matrix  that  we  introduced  through  Kqs.  (2.61  have  a  simple 
interpretation. 

Suppose  that  a  monochromatic  plane  wave  of  unit  ampli¬ 
tude  and  direction  of  propagation  specified  by  a  real  unit 
vector  n„,  i.e., 

,"o(r)  =  r  (IU) 

(with  time  periodic  factor  expl-iud)  not  shown),  is  incident 
upon  the  scatterer.  It  then  follows  from  Kqs.  (2.15).  (2.9), 
and  <2.;i)  that 

<" * *(n)  =  -  A(n  -  n„).  (3.2a) 

ik 

I)1  "l n)  =  A’(n,  n„),  (3.2b) 

ik 

the  upper  or  lower  signs  being  taken  on  the  left-hand  sides  of 
these  formulas  according  to  whether  ri,  *;  l).  Now  the  second 
integral  on  the  right  hand  sides  of  Kqs.  (2.1 )  represents  the 
outgoing  field.  I  n„)  say,  i.e.,  the  field  that  behaves  as  a 
diverging  spherical  wave  at  infinity.  Hence  it  follows,  on 
making  use  of  Kq.  (3.2b).  that  when  the  plane  wave  given  hy 
Eq.  (3.1)  is  incident  upon  the  scattering  object 

H-"(r;n(>  =  j  A  (n,  n„)e,,‘n  rdl!  when  r  i  /(' ,  (3,3a) 

A’ln.  n,,)c'*n'rdll  when  r  t  /Y*.  <3. 3b) 

where  a'  '  and  it'"  are  the  hemispheres  defined  hy  Kqs. 

(2.2). 

The  formulas  (3.3)  represent  the  outgoing  Held  in  each  of 
the  two  half  spaces  H  and  /t"  in  the  form  of  alt  angular 
spectrum  ol  plane  waves,  with  (generally  complex)  ampli¬ 
tudes  A(n.  nid  that  propagate  away  from  the  scatterer  in 
direct  ions  specified  by  unit  vectors  n  When  Ihe  c  compo 
nent,  tin...  ol  the  unit  propagation  vector  n,.  of  the  incident 
wave  is  positive  A  (n.  no)  clearly  has  the  physical  significance 
id  a  genera li;ctl  (runs mission  coc//icicn( .  (in,  n,,)  say,  w  hen 
ii,.  >  l)  and  of  a  gnicni/icerf  rc//cc(i on  coi7/iiient,  rln,  n, , ) 
say.  when  ti,  <  I),  /or  inmltnir  /tom  the  ho//  s/tticc  H  (see 
Fig.  3).  Recalling  expressions  (2.8a)  and  (2.8|»).  we  see  that 
these  coefficients  are  precisely  two  ol  the  elements  of  (he 
partitioned  .S'  matrix  (2.7).  vir.. 


(In,  n„) 

'in.  n(1). 

n .  >  0.  n0.  >  0, 

(3.4a) 

r(n,  n„) 

S' 

'(n.  n„i. 

n  <■'  0,  n,,  >  tK 

(.3  Hit 
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tln,n0l  ■  S1"  (n.n0) 


pin.iyj  S'"'in,nu) 


Fiji  3.  Illustrating  the  significance  of  the  elements  of  the  parti¬ 
tioned  .s'  matrix  as  generalized  transmission  and  reflection  coeffi- 


Similarly  the  quantities 

r(n,  n„>  S'  *  '(n.  n„).  n.  <  0,  n„.  <  0,  (3.4c) 


formulate  certain  generalized  Stokes  relations.  For  this 
purpose  we  first  separate  in  the  integral  that  expresses  the 
unitary  condition  (2.1Ja)  of  A(n,  n')  the  contributions  from 
the  hemispheres  <j,_i  and  <r,  +  l.  defined  by  Rqs.  (2.2): 

|  A’*(n,  n  )A’(n,  n")dl!  +  j  A’*(n,  n'M'(n,  n')d!2 

=  A(n'  -  n").  (4.1) 

if  we  chose  the  unit  vectors  n'  and  n"  to  have  positive  z 
components,  i.e.,  n.'  >  0  and  >  0.  and  recall  the  physical 
significance  of  A’(n,  n')  and  A(n,  n'')  discussed  in  Section  3, 
we  obtain  at  once  from  Eq.  (4.1)  the  following  relation: 


r*(n,  n')r(n,  n")d!2  +  t *(n.  n')2(n,  n")d!2 


=  A(n'  -  n").  (4.2a) 

Next  let  us  choose  the  unit  vectors  n'  and  n"  with  n.  <  0  and 
n."  >  0.  The  formula  (4.1)  then  gives,  if  we  also  use  Eqs. 
(3.4). 


pin,  n„)  S'*  "(n.  n„),  n,  >  0.  ft„,  <  0  (3.4d)  |  r*(n,  n')r(n.  n")dS2  +  |  ^(n.  n')/(n.  n”)dS2  =  0. 


have  the  physical  significance  of  a  generalized  transmission 
loiilmrnt  and  a  uencrnliicd  reflection  coefficient,  respec 
lively,  fur  incidence  from  the  half  s/niee  H'  (see  Fig.  3). 

Il  follows  (hat  in  terms  of  these  generalized  transmission 
and  relied  ion  coefficients,  the  partitioned  N  matrix  (2.7) 
may  lie  expressed  in  the  form 

i  flin.  n  )  gin,  n')l 
S(n.  n)=  ,  I-  (3.5) 

r( n,  n  )  r(n.  n  ) J 

l(  is  to  he  noted  that  in  view  of  the  relation  (2.13),  the 
generalized  transmission  and  reflection  coefficients  obey  the 
ret  iprneity  relations 

2(-n',  -  n)  =  r(n.  n'),  (3.6a) 

.-( -n  .  -  n)  = /(n,  n').  (3.6b) 

g( -n  .  -  n)  =  g(n,  n').  (3.6c) 


In  a  similar  way  we  obtain  with  tin1  choice  >  0.  n{"  <  0 
|  r*(n,  n')r(n,  n")dS!  +  j  2*(n.  n')p(n.  n")dS!  =  0, 


and,  with  the  choice  n/  <  0,  rt."  <  0. 
j  r*(n,  n')r(n,  n")di2  +  j  p*(n,  n  )p(n,  n")dl! 

=  A(n' -  n").  (4. 2d) 

One  can  readily  verify  that  the  four  relations  (4.2)  are 
equivalent  to  the  following  matrix  equation,  which  expresses 
the  unitarity  condition  (2.1  la)  in  terms  of  our  partitioned  N 
matrix  in  a  familiar  form: 


r(-n  .  -  n)  =  r(n.  n'l.  (3.6d) 

li  M  imis  worthwhile  to  point  out  that  our  definition  of  the 
Moralized  transmission  and  relied  ion  coefficients  depends 
,  ii  (In  choice  of  the  z  axis  II  is  possible  for  a  transmission 
i  oi  l  i icicnt  defined  with  respect  to  a  part ieiilnr  z  direction  to 
become  a  reflection  coefficient,  and  vice  versa,  when  the  z 
direction  is  chosen  differently  However,  in  many  situa¬ 
tions  ol  practical  interest  a  particular  direction  is  distin¬ 
guished  from  all  other  directions,  and  it  is  (hen  natural  to 
i  lioi.se  (lie  z  axis  along  this  special  direction.  Examples 
unhide  stratified  media,  inhomogeneous  plane  parallel 
plates,  and  rough  planar  surfaces  Scattering  f  rom  bodies  of 
irbitrarx  shape  in  the  presence  of  planar  phase  conjugate 
mirrors  also  belongs  in  (Ills  category 

4.  GENERALIZED  STOKES  RELATIONS 

V\  i  I;  i  he  interpretation  <■(  the  elements  of  (he  partitioned  .S' 
eia'tix  dial  wi  iu~t  discussed,  we  nr>  m  w  in  a  position  to 


j  S'(n.  n)S(n.  n")d!2  =  /A(n'  -  n’ ).  (4.3) 

Here  S’  is  the  Hertnitian  adjoint  of  Sand  /  is  the  identity 
matrix. 

In  a  similar  manner  that  led  to  the  relations  (4.2)  one  can 
derive  from  the  second  unitarity  condition  (2. 1 1  b)  of  Ain.  n  ) 
the  following  four  relations: 

j  p(n  ,  n)/i*(n ",  n)dS2  +  j  2<n  .  nl/*(n  ",  ntd!2 

=  A(n  -  n"),  ( 4.4a ) 

J  r(n  .  nlp’ln".  n)d’.2  +  |  r(n  .  nlf  *(n".  nidi?  =  (I. 

14.4b) 

I  /p(n  .  nlr’ln".  nidi’  f  j  f(n  ,  ntr’ln' ,  nldd  =  (V 

1 4.4c- 1 
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|  r(n\  nli*(n",  n)dll  4  j  r(n',  n)r*(n",  n)dii 

=  A(n'  -  n">.  (4.4d) 

The  relations  (4.4)  may  readily  be  shown  to  be  equivalent 
to  the  following  matrix  equation,  which  expresses  the  second 
unitary  condition  (2.1b)  in  terms  of  our  partitioned  S  ma- 

S( n',  nlS'ln",  nldfl  =  MOT,  it*).  (4.5) 

The  formulas  (4.2)  and  (4.4)  may  be  regarded  as  general¬ 
izations  of  the  Stokes  reciprocity  relations  (1.1).  We  verify 
in  Section  5  that  they  reduce  to  Kqs.  (1.1)  in  the  special  case 
considered  by  Stokes. 

Of  the  eight  relations  (4.2)  and  (4.4)  only  two  are  actually 
independent  of  each  other.  To  see  this  let  us  first  apply  the 
reciprocity  relations  (3.6a)  and  (3.6c)  to  Eq.  (4.4a)  and  take 
the  complex  conjugate  of  the  resulting  equation.  This  gives 

I  p*(-n,  -  n')p(-n,  -  n")dtl 


+  J  r*(— n,  -  n')r(-n,  -  n")dS2  =  A(n'  -  n").  (4.6) 

If  we  now  change  the  variables  by  letting  n  -  -n,  n'  *  -n', 
and  n”  -  -n",  the  relation  (4.6)  becomes 

j  p*(n,  n')p(n,  n")dl)  +  j  r*(n,  n')r(n,  n")dSl 

=  Mb'  -  n"),  (4.7) 

which  is  the  relation  (4. 2d).  In  a  strictly  similar  manner  one 
can  show,  with  the  help  of  the  reciprocity  relations  (3.6),  that 
Kqs.  (4.4b),  (4.4c),  and  (4.4d)  are  equivalent  to  Kqs.  (4.2c), 
(4.2b),  and  (4.2a),  respectively.  Hence  if  we  take  the  reci¬ 
procity  relations  (3.6)  into  account,  the  set  of  the  four  equa¬ 
tions  (4.2)  contains  the  same  information  as  the  set  (4.4). 
We  may,  therfore,  confine  our  attention  from  now  on  to  the 
set  (4.2)  only. 

Since  the  two  half-spaces  i<~  and  /?*  play  the  same  role  in 
the  present  t  heory .  it  is  clear  that  from  any  of  the  generalized 
Stokes  relations  that  we  just  derived  one  will  obtain  a  valid 
relation  through  the  simultaneous  transformations 

al+,-*<r  (4.8a) 

t--T.  (4.8b) 

P**r.  (4.8c) 

Two  formulas  that  transform  into  each  other  in  this  way  may 
be  said  to  be  dual  of  each  other.  Clearly  Eqs.  (4.2a)  and 
(4. 2d)  form  a  dual  pair,  and  so  do  Eqs.  (4.2b)  and  (4.2c). 
Hence  there  are  essentially  only  two  independent  relations 
of  the  type  that  we  are  considering,  which  we  may  take  to  be 
Eqs.  (4.2a)  and  (4.2b).  The  other  six  relations  may  be  ob¬ 
tained  from  them  bv  the  use  of  reciprocity  and  duality. 

5.  AN  EXAMPLE:  STOKES  RELATIONS  FOR 
STRATIFIED  DIELECTRIC  MEDIA 
SURROUNDED  BY  FREE  SPACE 

We  will  illustrate  the  use  of  the  general  relations  (4.2a)  and 
14. 2b)  by  applying  them  to  the  interaction  of  a  plane  mono 


chromatic  electromagnetic  wave  with  a  stratified  dielectric 
medium. 

Consider  a  si  ratified  dielectric  medium  that  occupies  the 
strip  0  5  z  <  /,,  and  let  /V  =  Nlz)  be  the  (real)  refractive 
index  function  of  the  medium.  We  assume  that  the  strati 
fied  medium  is  surrounded  by  free  space;  hence  N{;)  =  1 
when  z  <  0  and  when  z  >  L.  We  assume  further  that  the 
incident  electromagnetic  wave  is  linearly  polarized,  with  its 
electric  field  either  in  the  plane  of  incidence  ('I'M  wave)  or 
perpendicular  to  it  (TE  wave).  As  is  well  known  |Kef.  12, 
Sec.  1.6.1]  the  state  of  polarization  of  either  of  these  two 
waves  (modes)  does  not  change  on  interaction  with  the  strat 
ified  medium;  and  an  incident  wave  of  any  state  of  polariza¬ 
tion  may  be  expressed  as  a  linear  combination  of  these  two 
modes,  which,  moreover,  are  independent  of  each  other 
when  they  interact  with  the  stratified  medium. 

A.  Consequences  of  Eq.  (4.2a) 

Suppose  first  that  the  wave  is  incident  upon  the  stratified 
medium  from  the  half-space  z  <  0,  in  a  direction  specified  by 
a  unit  vector  n'(n,'  >  01,  and  let  nr(n’)  and  a, In')  be  the  unit 
vectors  in  the  direction  of  propagation  of  the  reflected  and 
the  transmitted  waves,  respectively  (Fig.  4(a)}.  The  func¬ 
tional  dependences  of  n,  and  of  n,  on  n'  are  given  by  the  laws 
of  reflection  and  refraction,  respectively,  for  stratified  media 
(Ref.  12,  Secs.  1.6.1  and  1.6.3).  Since  there  is  only  one 
reflected  and  one  transmitted  plane  wave,  the  generalized 
reflection  and  transmission  coefficients  will  evidently  be  of 
the  form 

r(n,  n')  =  r(n')A|n  -  nr(n')|  (5.1a) 

and 

t(n,  n  )  =  l(n')A|n  -  n,(iT)],  (5.1b) 

where  r  and  f  are  the  usual  reflection  and  transmission  coef- 


\*u.V  '(!(’•  Ml  Im  iV 
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Fitf.  4  lllustrat  mg  the  nutation  relating  to  the  derivation  nt  tin 
Stokes  relations  lor  stratified  dielectric  media 
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tn  ii'iil.-.  respect ivelv.  tor  incidence  tor  the  half-spare  z  <  0. 
I..  !  us  consider  the  first  term  on  the  left  hand  side  oi  the 
generalized  Stokes  relation  (4.2a).  If  we  Use  Kq.  (fi.la)  we 
hav  e 

r‘m.  n  ir(n  .  n")dl! 


-  i  r*(nlr(n"lA|n  n,(n  l|A|n  -  u  in")|di:  to. 2) 

i  rt  oi.  el.  in',  ).  and  (a",  y "  I  he  the  sphero  id  polar  angles  ol 
• ‘.nil  veelors  n.  n  .  and  n",  with  the  polar  axis  being  taken 
g  Hie  positive  ;  direction.  Then  we  have,  according  to 
!  hi  law  ot  reflection. 

o a  n  i  '  s  —  e  ).  i i.fn " '  '  '  ir  -  y" t,  i.i.’U 

•  ml  Kq  i  'i.2l  becomes 

■  'in.  n  irin  ,  n " tdSf 


nirin")  A|n  -  n.tn')|A|n  -  nr(n")|dli.  to. 4) 


. i  integral  on  the  right  hand  side  may  be  evaluated  at  once 
ie  esc  ot  Kq.  it '.4 1  of  Appendix  t’  and  the  fact  t  hat  n.  <  0. 
In  -  gives 

A|n  n  in  i)A|n  -  n. i n ” >J«1 12  =  A|n,ln)  -  n,(n")|. 


ol  the  delinition  i .’ d  'l  o|  itn  sphere  a I"  delta 
,iiid  make  use  , ,f  Kq.  id  o  .nm  n,i-  fa,  t  that  a  -  U, 
it:  mi  Kq  Id  fit  t he  relat e a, 


A[n  -  n.in  i|A|n  -  n.in  ;jd'..' 

•'•[  --"'-in-  ""  I -  y'l 
sin  " 

-Ain  -  n  e  u'l  i.i 

di  - '  .1  at  mg  I  roll  i  Kq  I  d.ti)  into  Kq  '  d  1 1  we  tind  that 

•'ll  n  o.n  n"hU!  -  e*m  irin"iAin  ~n  I  id.,') 

,  -'rati,  siiiniar  manner  We  in  i\  tvalu.ile  the  seiolid 
1  u  n  tie  li  P  ham)  side  o|  Kq.  '  I.L’al.  We  then  hnd 

'  n  n  n  n  nl'.'  :'in  inn  '.Vn  n  i  id. hi 

t  -M.h.i  in loi.  o.aii  Kqs  i  ,  Vi  and  id. Hi  into  (he 

,  i  ■ .  0,  si , ,  i  eh*  1  a  a  i  i  1  .  a  I  and  integrating  I  ml  h  sides 

'  1  ■  i ;  i  -phi  re  ge  i  a  l  a  I  ed  1  w  I  he  \  eel  <  a  n  "  -  l ,  y  '  l.  we 

■  ,  ..  Is  id. 

•■ft  >nn  I  )  .'•»!»  ">/»»»  )  ~  I  lFl.it) 

■'  ■  i  i  ,  \  id,  n  i  A  e  >  presses  ,  ,  a  r ■  iifi.tr;  ,  •/  ,  firq'i 

'•  •'  'I,  iieipoiihi  ie,  qiroi  n\  relat  ion  i  >  Kb)  t  tie  torniii 
•  •  . .  •  •>  id,U  in  i  \  pri  -  ,  Mi  to  r  in  •  it  a  Stokes 

a  a  .• .  a i  I.,,  V  a ,  I,.,  t  ■  is  purpose,  and 

r  i.a  i  r  oi  •  .  ii.,-  in  i  w  n  n  a  iiere.  at  Ion  ol  allot  her 


Stokes  rel  it  a i;'  os  t ;  i 
wave  prop  i/a ' .  •  ,r  n 
spe,  it  fed  i  .  loi,'.  . 
analog-.'  Willi  l.q  i  '  ',  . 


.eeie  the  incident 
the  direction 
•  /  mild  have,  by 


>n  '  ii  -  n  t  n  i 


t d job1 


1 1  ere  ,  anil  .  no  1 1  .  ,  -  i  a  1 1  ,  e  <  o  a  al  I  raiisniission  coelf  i 


I  1**1  1 1  K  l»S|  I.  l  •  i  J  ,  i  .  i'll  ,  !)  It- 

■  ■Ikl  f»  I  n  ).  ill  >< !  rt  i.  t.'  i  I  I.I  ■  r.  I  . 

relict  led  lin. I  Ii  ,.o:  I..  I, .  d  .  , 

determine  it.  it  ■ .  i  n  , 

i  i'-'elv 

W  I-  till v  it,.  ,  ,  • 


'  he  bull  spill  e  ■  -  I 

,  I :  ••  I  rat  ions  ol  the 
•■■  q  *  ,ih  h  may  be 
.  i  1  .  ,r,.  respec- 


Xow  s.iu  r  tile  - 1  i  i  e,,  i  i.  i  i.i!  to  be  sqr 

rounded  on  tioi i.  side.,  nv  l rt  ,  <•  ,  ,  z,  fraction  for 

stratified  media  | lift .  1  He.  la,  ,.  |q  ui'|  gives  at  once 
h  ( -n)  =  -  n.  alio  l  In  loti,  ’  i  >  becomes 

,  i- n',  n i  -  i  - n l Al  -  :t  i  til  <’>. 1 2) 

li  also  lollows  <r.  ,e  i.  ,  t‘  it  ■  aid  the  relations  n,(n  i  =  n‘ 

and  n  i/i  'i  n  •  .'tin  t.  .q.  ,o.  i •  .| |,  ••  it  mice  from  the  law  of 

refract  ion  l  that 

i  n  ..  ,ii  ,  . '  ■  i  ii  i  i  d  l.ii 

On  sided  itiitoip  . .  .  K  |  i  I  a  mu  ,  d.  12)  into  the  reciprot 

it  v  relation  i.i  l.i.i  and  . 1. 1 i  al  ,i,g  ,n  ,-r  I  hi'  unit  sphere  gen 
crated  by  the  via  (•  it  ai  lii.if  (fiat 


i  ni  mil 


Tins  lormiila  old  ious,y  <  *j„  c  ■  a  t  • ,  [m  tninsmis 

S tori  hi  II  tliuh/tt  ,/  ■••<«  than, 

On  substiiutii',:  lion.  Kq  ,■  i  mo. i,  I  j  (ah)  we  obtain 
the  relation 


I  ills  ,o,  is  #  1  ...  •  ;or  stratified 

•fididri.  n,o/ii,  . . . lit,,  i. ,,  i.  11  known  in  the 

'.la  or .  ol  ;  I.I,  1,.  t„.  . , ,  id,,,  iff  ;  ,  I  m 

It  I  .01.1"  qta  ,,.  1  ail  L,|  (I  aj 

Next  Wi  .,ili  „|.  I  lie  ol.|  In  alii.,  mil,  e  r.tlized  Stokes 
relation  il  21.)  Hnppo.,  that  die  pi,, a,  t.ve  is  incident 

lip. a ,  tin  . . |  li.edlo...  O  .11  -la  ||  .1  pail  "•  /..  ill  a 

•  lire,  tiiai  -  j . .  1 . .  .  l.l  m  I  In  siib-li 

t ,  it  i  ■  ip  fiom  Kqs  ih  1 1 1 1  - 1  ,.o  a  i '  i  .i  ,ii.  I-  ..,st  iniegriil  on 
l  111  It'll  hand  s,,|,m  ,  a  (•  q  ..  i  |  '  id  a  l  .1  o  I  f  •  d 

I  i  *  i  n  li  nl  i,  n  I 


.‘ill  elf,  j  ql.  ,  I.  q  ,,n  tl  I  ft"  l|dii.  id. lti) 

I  i.c  mnvi  i)  i.  . .  .  i  ,,.l  .  .  ,  .  .  .  .x  •), Iv  he  evaluated 

),\  the  os,  al  |v ,  it'-.,.-:  i.  n x  '  in  I  the  f  ,.  t  that  n,  <  0 

II  we  al-o  'oak'  is  a  llu  l  i  iiiiit  i  i,i  n  ,  w ha  h  follows 
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r*(n,  n'Mn,  it")di> 


Ism  O'  I 


In  a  similar  way  one  may  evaluate  the  second  integral  on 
I  he  right  -hand  side  of  K<|.  (5.16),  also  making  use  of  the  fact 
that  with  n'  ((/',  y'),  n,.(n')  =  (ir  —  O’,  y').  One  then  readily 
finds  that 


e*(n,  n'lttn,  n")dSt 


=  j»*(n')t(n") 


Mir  +  ii"  -  „)Mf'  -  «-”) 
Isin  O'  I 


(5.18) 


On  substituting  from  Kqs.  (5.17)  and  (5.18)  into  the  general¬ 
ized  Stokes  relation  (4.2b),  integrating  both  sides  of  the 
equation  over  the  unit  sphere  generated  by  the  unit  vector  n' 
=  nr,  v-'),  and  making  use  of  Kq.  (5.2),  we  find  that 

r*lnr")r(n">  +  ,,*(n/')l(n")  =  0,  (5.19) 

where  n,"  nr(n")  is  the  unit  vector  along  the  direction  of 
the  reflected  wave  when  the  incident  wave  propagates  in  the 
direction  specified  by  the  unit  vector  n"  In/'  >  (1). 

If  in  K*|.  (5.19)  we  make  use  of  the  reciprocity  relation 
(5.14)  and  write  n'  in  place  of  n",  we  obtain  the  formula 

t*(-n/)r(n')  +  g*(n/)j(n')  =  0.  (5.20) 

There  tire  other  forms  in  which  the  relations  (5.19)  and 
(5.20)  can  he  expressed.  For  example,  since  the  unit  vectors 
n'  and  -nr'  make  the  same  angle  with  the  z  axis,  t(— n,')  = 
t(n').  For  the  same  reason  />(n/)  =  p(-n').  Making  use  of 
these  relations  in  Kq.  (5.20),  we  obtain  the  formula 

ftn')r(n')  +  Hn'bX-n')  =  0.  (5.21) 

This  formula  is  another  Stokes  relation  for  stratified  dielec¬ 
tric  media  and  is  of  a  form  well  known  in  the  theory  of 
dielectric  films  (Ref.  2,  p.  172). 

The  relations  (5.21)  and  (5.15)  recently  played  an  impor¬ 
tant  role  in  the  theory  of  correction  of  distortions  hv  the 
technique  of  phase  conjugation.1 


APPENDIX  A:  DERIVATION  OF  A  FORMAL 
ASYMPTOTIC  APPROXIMATION 

\4  e  begin  by  recalling  that  under  very  general  condil ions  any 
solution  Vtr),  valid  throughout  the  whole  si  ace,  of  the 
Helmholtz  equation 

v  T'(r)  f  h’Vtr)  =  I)  (Al) 

may  be  expressed  in  the  form  of  an  angular  spectrum  of 
homogeneous  plane  waves,  all  with  the  same  wave  number  h, 
that  propagate  in  all  possible  directions1': 

V  ( r)  =  I  ntnle'^MS!.  <A2) 

The  complex  spectral  amplitude  function  «(n)  can  be  de 
rived  from  the  knowledge  of  V '( r )  by  the  inversion  formula 
|Ke|  |:t,  Kq  (Hi  1)| 


o(n) 


k~  ,  lim  f  dK  V(r)e  ,Kn',d:V. 
(2xV*  •  •*>>  lk-> 


(Ad) 


The  asymptotic  behavior  of  Vtru)  as  hr  *  w  in  the  direc¬ 
tion  of  a  real  unit  vector  u  may  be  obtained  by  the  use  of  Kqs. 
(2.1)  and  formula  (2.2)  or,  somewhat  more  directly,  from  a 
mathematical  lemma  due  to  Jones.11  The  result  is 


Vtru) 


o(-u) 


as  hr  •■<». 


(A4) 


Let  us  now  apply  these  results  to  the  case  when  V(r)  is  a 
plane  wave  of  unit  amplitude  that  propagates  in  the  direc¬ 
tion  of  a  unit  vector  n„: 


Vtr)  =  , 


iAtn|l-r 


( A5) 


To  determine  the  angular  spectrum  amplitude  function  otn) 
of  this  field  we  first  note  that,  when  V(r)  is  given  by  Kq.  ( A5), 
the  Fourier  transform  that  appears  in  Kq.(AJ)  becomes 

I  Vtr )e  lKn  ,d  V  =  j  exp|i(Ann  -  Kn)  •  r|d  'r 

=  (2ir)  V  "(frn„  -  K'n),  (A6) 


where  A,:n  is,  of  course,  the  three-dimensional  Dirac  delta 
function.  On  substituting  from  Kq.  (A6)  into  Kq.  (A2)  we 
find  that  <i(n)  is  now  given  by  the  formula 

ft*. 

n(n)  =  k~  lim  .V  ”(/;n„  -  Kn)d K.  (A7) 

-  h.  , 

To  evaluate  the  integral  on  the  right-hand  side  of  Kq.  (A7) 
we  make  use  of  the  representation  of  the  three  dimensional 
Dirac  delta  function  in  spherical  polar  coordinates. ,r'  One 
then  finds  at  once  that 

fimlhn„  -  Kn)  =  (l//r)A(n  -  n„)<)(/.'  -  K),  ( A8) 

where  A  is  defined  by  Kq.  (2.12)  and  A  Ik  -  K)  is  the  one- 
dimensional  Dirac  delta  function.  On  substituting  from  Kq. 
(A8)  into  Kq.  (A7)  and  carrying  out  the  trivial  integration 
with  respect  to  K,  we  find  that  the  angular  spectrum  ampli¬ 
tude  function  of  the  plane-wave  field  is  simply 

otn)  =  A  In  -  n„).  (A9) 

Finally,  on  substituting  from  Kqs.  (A5)  and  (A9)  into  the 
asymptotic  formula  (A4)  we  obtain  the  formal  asymptotic 
approximation 


ilov.-r  2ir 

I’  ^ 

ih 


A(n  -  n,,) 


Atn  +  n„) 


as  hr  »<*■. 


(Aid) 


It  is  of  interest  to  note  the  form  of  the  angular  spectrum 
representation  of  a  plane  wave.  It  is  dear  on  comparing  the 
right  hand  sides  of  Kqs.  (Alt))  and  (2.2)  that  for  the  plane 
wave  defined  by  Kq.  (A5) 

("*'(n)  =  -  Atn  -  n„l,  (Alla) 

ih 


/)"'( n) 


2n 

,h 


Atn  -  n„). 


(Alibi 


On  taking  in  (Kq.  2.1 )  Mr)  =  exptiAn,,  ■  r)  and  on  suhsl  it  lit 
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ing  for  the  C  and  D  coefficients  the  expressions  (All),  we 
find  that 


F,(n)  =  ~  4’(n,  n0). 


(B6) 


In  H~: 

=  I  A(n  -  n„)e‘*"rdS2  +  j  A(n  -  n0)e'',ordfi, 

(A12a) 


In  VU: 


ln"'mj  A(n  •  nllli',l,"ill!  *  |  Atn  -n„)r'*"MW. 

(A  12b) 


One  can  verify  by  direct  evaluations  of  the  integrals  on  the 
right-hand  sides  of  these  equations  that  these  formulas  hold 


The  formulas  (B6)  and  (B4b)  are  Eqs.  (2.15)  of  the  text. 


APPENDIX  C:  DERIVATION  OF  THE 
FORMULA  /4*  A(n  -  n')A(n  -  n")dfi  =  A(n'  -  n") 

We  have,  according  to  the  definition  (2.1 2)  of  the  “spherical” 
Dirac  delta  function 


where  (»,  y)  and  (O',  y')  are  the  spherical  polar  angles  of  the 
unit  vectors  n  and  n’,  respectively.  We  have  a  similar  ex¬ 
pression  for  A(n  -  n").  Hence  it  follows  that 


I  A(n  -  n')A(n  -  n")dfi  =  .  -.- — : 

lain  tf'l 


6(0  —  0r)6(^  ~  *p')M0  —  0")6(ip  ■ 


(sin  6 1 


ic") 

- sin  ddtfdy. 


(C2) 


throughout  a  wider  domain  than  indicated  here;  in  fact,  each 
of  the  two  Eqs.  (A12)  is  a  valid  representation  of  the  plane 
wave  exp(-ifcn<i  •  r)  throughout  the  whole  space. 

APPENDIX  B:  DERIVATION  OF  FORMULAS 
(2.15) 

Suppose  that  the  field  incident  upon  the  scatterer  is  a  plane 
wave  of  unit  amplitude  that  propagates  in  the  direction  of  a 
unit  vector  n„: 

D‘"»r>  =  o‘ha"".  (Bl) 

The  total  field  (incident  +  scattered)  in  the  far  zone  is  given 
by  a  formula  of  the  form 

t’(rn)  ~  cl<",°r  +  A(n,  nn)  - — •  asfer-*®,  (B2) 

where  A(n,  n„)  is  the  scattering  amplitude.  If  we  substitute 
in  Kq.  (B2)  for  expU'fenn  •  r)  its  formal  asymptotic  approxi¬ 
mation  given  hv  formula  (A10),  the  expression  (B2)  for 
i  '(rn)  acquires  the  form  (2.9),  viz., 

-ikr  „-l  hr 

( '(rn)  ~  F,(n)  —  +  F,(n)  —  -  as  kr  -  ».  (B3) 
r  "  r 


where 

F.ln)  =  A(n  -  n0)  +  A(n,  n„),  (B4a) 

iff 

F,tn)  =  -  “/  A(n  +  n,,).  (B4b) 

ik 

The  expression  (B4a)  may  readily  he  expressed  in  terms  of 
the  -V  matrix.  To  do  so.  we  substitute  from  Eqs.  (B4)  into 
tin-  formula  1 2. 10)  that  may  he  regarded  as  a  definition  of  the 
>'  matrix.  We  then  find,  after  trivial  calculation,  that 

A(n.  n„)  =  |T(n.  n„)  -  A(n  -  n„)|.  (B5) 

i)r 

On  comparing  Eqs.  (Bo)  and  (B4a!  we  see  at  once  that 


Now  for  0  <  ft  <  ir.  isin  0\  =  sin  0,  and  Eq.  (C2)  therefore 
reduces  to 

I  A(n  -  n')A(n  -  n'')di!  =  I  6(0  -  0')6(0  -  0")d0 

J(4r>  J() 


•2* 

x  (Mv  -  y')6(v  -  v"ldv.  (C.3) 
Jo 


By  an  elementary  property  of  the  Dirac  delta  function  (Ref. 
12,  App.  IV,  Eq.  (12)]  the  first  integral  on  the  right  hand  side 
is  equal  to  MO'  -  0”)  and  the  second  to  Me'  -  v").  Using 
these  facts,  Eq.  (C2)  reduces  to 


f  A(n  -  n')A(n  -  n")dil  = 

;|4«| 


M0’-0")  «(*'-*"> 
Isinrl  1 


or,  recalling  again  the  definition  of  the  “spherical”  Dirac 
delta  function  [see  Eq.  (Cl)], 


A(n  -  n')A(n  -  n")dS!  =  A(n'  -  n"). 


(C5) 
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Abstract.  The  transmission  of  a  generalized  radiance  across  a  planar  boundary 
separating  two  homogeneous  media  is  considered.  It  is  assumed  that  the  optical 
field  remains  continuous  at  the  interface  and  reflection  is  neglected.  A  result  is 
obtained  which  may  he  regarded  as  a  generalization  of  the  conventional  radiance 
theorem  for  fields  of  any  state  of  coherence.  This  result  differs  from  the 
conventional  theorem  by  a  factor  that  depends,  in  general,  both  on  the  optical 
intensity  and  on  the  degree  of  coherence  of  the  field.  However,  oxer  a  wide  range 
of  circumstances  the  generalized  radiance  theorem  is  shown  to  he  in  good 
agreement  with  the  conventional  theorem. 


1.  Introduction 

One  of  the  basic  principles  of  conventional  radiometry  [1]  is  the  so-called 
radiance  (or  brightness)  theorem  that  pertains,  in  its  most  general  form,  to  the 
relationship  between  the  radiance  of  an  object  and  the  radiance  of  its  image  formed 
by  any  specular  optical  system. t  Within  the  framework  of  linear  theory,  an 
arbitrarily  complicated  specular  optical  system  may  be  considered  simply  as  being 
composed  of  a  sequence  of  uniform  media  separated  by  sharp  boundaries.  The 
conventional  radiance  theorem  then  follows  directly  from  the  phenomenological 
laws  that  govern  the  transmission  of  the  radiance  through  a  uniform  medium  and 
across  a  boundary  separating  two  uniform  media  with  different  indices  of  refraction. 

The  propagation  of  the  conventional  radiance  is  governed,  under  general 
circumstances,  by  the  equation  of  radiative  transfer  ([2],  chapter  I ,  equation  (47)).  1 1 
implies  that  in  a  uniform  medium  (that  does  not  contain  sources  or  absorbers)  the 
radiance  function  /J,„(r,S)  at  some  frequency  to,  measured  in  the  direction  specified 
by  the  unit  vector  s,  remains  invariant  on  the  line  in  the  direction  S  through  the  point 
represented  by  the  vector  r.  In  a  number  of  recent  publications  (see,  for  example, 
[3  XJ,  the  validity  of  the  equation  of  radiative  transler  has  been  investigated  in  a 
(statistically)  homogeneous  medium  with  scalar  fields  of  arbitrary  states  of 
coherence.  The  discussion  has  also  been  extended  into  the  domain  of  electro¬ 
magnetic  fields  both  w  ithin  the  framework  of  classical  |  9,  10]  and  quantized  [11,12] 
wave  theories. 


t  A  specular  optical  system  in  this  context  is  one  that  does  not  contain  diffusely 
transmitting  (or  reflecting)  surfaces. 
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The  transformation  of  the  radiance  across  a  sharp  boundary  has,  however, 
received  considerably  less  attention.  Although  the  reflection  and  refraction  of  wave 
fields  at  boundaries  have  been  extensively  studied,  no  work  is  known  to  the  author 
that  deals  specifically  with  the  transmission  of  the  radiance  associated  with  a 
fluctuating  optical  held  through  a  medium  discontinuity  surface. +  For  this  reason  we 
concentrate  in  this  paper  on  examining  the  conventional  radiance  theorem  and  its 
range  of  validity  with  partially  coherent  light  at  a  single  planar  boundary  between 
two  homogeneous  media. 

In  the  present  paper  we  will  adopt  a  relatively  simple  and  straightforward 
approach  that  is  based  on  the  scalar  theory  of  light.  The  fundamental  assumption,  as 
is  customary  in  physical  optics,  is  that  the  (monochromatic)  optical  field  remains 
continuous  across  the  interface.  Nonetheless,  the  method  employed  takes  into 
account  some  interesting  physical  phenomena  such  as  the  conversion  of  evanescent 
waves  into  propagating  plane  waves.  Moreover,  it  offers  several  valuable  clues  to  a 
future  improvement  of  the  analysis. 

2.  Radiance  theorem  of  conventional  radiometry 

We  begin  by  briefly  recalling  the  phenomenological  form  of  the  radiance  theorem 
as  it  is  traditionally  encountered  in  radiometry.  Conventional  radiometrv  deals  with 
the  problem  of  energy  transport  at  some  temporal  frequency  «>.  With  reference  to 
figure  I,  the  conventional  radiance  theorem  at  a  single  refracting  surface  may  be 
expressed  in  the  form  [14] 

«,( r.s,)  _  fl,(r.*2) 

7  ’  2  ~  *  '  *  * 

"I  ”2 

where  /{,  and  ll2  are  the  values  of  the  radiance  on  the  two  sides  of  the  interface  (the 
explicit  d)  dependence  is  omitted),  w,  and  n2  are  the  refractive  indices  of  the  two 
uniform  media,  and  r  denotes  the  position  vector  of  an  element  Jer  of  the  boundary. 
Further,  it  is  important  to  note  that  m  the  conventional  radiance  theorem  ( 1 )  S,  and 
S2  are  unit  vectors  that  specify  the  path  of  a  geometrical  light  ray  across  the  surface 
element.  The  effects  of  reflection  have  been  neglected  in  the  derivation  of  equation 
0). 


Figure  I  Illustration  of  the  notation  relating  to  the  radiance  theorem  of  conventional 

radiometry. 

4  1'he  propagation  ol  generalized  radiance  functions  (that  pertain  to  fields  of  any  state  of 
coherence)  in  lens  systems  has  hecn  studied  in  [4]  and  [X  j  on  the  basis  of  the  usual  Fourier 
optu  s  approximations  Also,  a  rather  general  (one -dimensional )  analvsis  of  the  propagat  ion  of 
a  generalized  radiance  through  lenses  was  presented  in  [13]  making  use  ol  asymptotic 
approximations  fgcomrtrn  al  optics  limit)  based  on  the  principle  ol  ihe  stationary  phase. 


liiiiluniti  t/u  iiii  in  fctth  partmlh  i  ti/irH  iit  huh/  1371 

The  radiance  theorem  (I)  of  conventional  radiometrv  expresses  merely  the 
conservation  of  energy  (at  frequency  to)  that  is  incident  on  the  element  of  do  from  the 
differential  solid  angle  dll,  around  the  direction  St  and  emerges  into  the  differential 
solid  angle  df}2  around  the  direction  S2  (see  figure  1).  Therefore,  it  is  clear  that  in 
connection  with  the  conventional  radiance  theorem,  the  size  and  direction  of  the 
element  dfl2  are  directly  determined  by  dQ,  through  Snell’s  law  of  refraction. 

The  losses  due  to  reflection  could  be  included  in  the  conventional  radiance 
theorem  ( 1 )  by  introducing  a  phenomenological  coefficient  of  reflection  that  depends 
on  both  position  and  direction.  Formally  this  would  involve  the  use  of  the  differential 
scattering  coefficient  ([2],  chapter  1 ,  §  3)  that  appears  in  the  conventional  equation  of 
radiative  transfer.  Approximate  values  for  the  reflection  coefficient  could  be 
obtained,  for  example,  from  measurements  or  from  the  customary  Fresnel  equations 
(see,  for  example,  [IS]). 


3.  Generalized  radiance  and  the  radiance  theorem  with  light  of  any  state 
of  coherence 

In  the  context  of  fluctuating  optical  fields,  geometrical  optics  cannot  be  used  to 
couple  the  energy  transport  on  the  two  sides  of  the  boundary.  For  this  reason  we  will 
make  the  assumptions,  common  in  physical  optics  [16],  that  the  optical  field  remains 
continuous  in  passing  across  the  boundary  surface  and  that  on  the  surface  it  is  given 
simply  by  the  incident  field.  These  assumptions  are  the  cornerstones  of  the 
customary  analysis  of  scalar-wave  propagation  in  optical  systems,  w  here  the  various 
elements  such  as  lenses  are  represented  by  complex-amplitude  transmission 
functions.  Clearly,  the  assumed  field  properties  then  also  imply  that  reflection  at  the 
interface  has  been  neglected.  +  Since  the  continuity  holds  for  each  realization  of  the 
statistical  ensemble  (assumed  to  be  stationary),  the  cross-spectral  density  function 
[18]  lF(r, ,  r2)  that  characterizes  the  spatial  coherence  properties  of  the  field  at 
frequency  (U,  will  also  be  continuous  across  the  boundary. 

For  the  sake  of  simplicity,  we  take  the  refracting  surface  separating  the  two 
homogeneous  media  to  be  a  plane  a  =  constant,  say  s  =  su,  and  consider  a  wavefield 
propagating  across  the  boundary  into  the  half-space  z  >  s„  (figure  2).  We  may 
then  associate  with  the  field  distribution  in  any  transverse  plane  s  =  constant  a 
generalized  radiance  function  defined  by  the  expression^  ([16],  equation  (21)) 


H{p,  s  )  —  {k  in)1  cos  0 


!('(/>  I  1/2 p,p  1  2/»')exp(  i/es  * />" ) d 2  /»’ . 


(2) 


where  IF  (p  +  1/2 p’,  p  -  ’  2 p  )  denotes  the  cross-spectral  density  (at  frequency  o))  of 
the  light  at  the  points  Pi’  =  p  +  1  Up  and  p2  =  p  —  I  jlp  in  that  plane,  anti 

l<  -  nk0  =  n(it);c),  (3) 

+  In  scalar  optics  one  sometimes  requires  that  froth  the  optical  held  and  its  normal 
derivativ  e  remain  continuous  across  a  sharp  boundary.  These  boundary  conditions  then  give 
use  also  to  a  reflet  ted  held  component.  In  particular,  \\  itli  a  planar  boundary  and  an  incident 
plane  wave,  the  resulting  coefficients  for  reflection  and  refraction  are  the  usual  Fresnel 
equations  for  the  case  when  the  electric  field  is  perpendicular  to  the  plane  of  incidence 
i  compare,  tor  example,  [17],  equations  ( 5)  and  (6),  and  [  1  5],  equations  (4.34)  anti  (4.35)). 

j  Since  the  (.  artesian  components  >,, and  .*.  ot  the  unit  \  ector  s  are  related  In  the  itletit  ii  \ 
,»»;  +  (-  l ,  only  two  of  the  three  components  are  independent.  We  will,  therefore,  regal tl 
die  generalized  radiance  as  being,  in  its  directional  dependence,  a  function  of  t lie  two- 
hmcnsional  transverse  vector  s  (*,,*, )  and  denote  the  radiance  by  Hip,  S.  ). 


’  *  *  »  *  «  */  *  »  *  I  * 
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with  n  being  the  refractive  index  of  the  medium  and  r  the  speed  of  I, gin  ,r.  c.ti  uurn. 
Further,  in  equation  (2),  0  denotes  the  angle  between  tin  unit  \e. ■;<  t  i  t.u  tie- 
positive  s  axis  and  S  is  the  projection  of  S  (consideied  ax  a  wwi-dui"  ,,  .  ti  t  r  . ,  , 

onto  the  plane  -  constant,  i.e  it  S  t _  t  tin  1 1  h  (>.  i  I  n 

Finally,  the  integration  m  equation  (2)  extends  thronglm.i  1 1  it  >  ■  i  - , . .  (  i 

z  =  constant. 


medium  1  '< 

j 
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Figure  2.  A  planar  boundary  :  —  z, ,  separating  two  hon.ugcncoip  media  with  lefri  live 

indices  «,  and  n2,  respectively. 


For  later  analysis  it  will  be  convenient  to  introduce  an  au\ui.,rv  quantity,  know n 
as  the  \Vigner  (distribution)  function  [20],  that  is  closely  related  m  tin-  •  nnnei  ru'd 
definition  of  the  generalized  radiance  given  in  equation  (2i  k\  e  q,  lira  On  \\  igm  i 
function  by  the  formula 

!* 

X(P,f)  =  J  U'(P  +  1/2 p'.p-  l/2p  ) exp (  -  if  •  p  id', i  l  [I 

where  f  is  a  (real)  two-dimensional  vector.  An  analogous  quantity,  pc  h  ..tin  lie 
identified  [21]  in  some  sense  with  the  ‘local  spatial-frequency  >|ic»  limn  o|  tl.,  |u  id, 
has  been  studied  extensively  by  Bastiaans  [22,23]  in  geometrical  and  ii|uteiu 
physical  optics. 

On  comparing  equations  (  2)  and  (4),  we  see  at  once  that 


/ftp.  st  (/,■  2tr t  “  I  »S  tl/l  p  I  S  I  '  I 

a  relation  which  is  valid  in  any  transverse  plane  a  -  con-  uiiti.  Wi  >■  m  ,,  „i  . 
particular,  two  planes  displaced  by  a  small  distance  t:  ft mn  t|.,  ,  |  inc  ,.ip, 

media  1  (refractive  index  w, )  anil  2  (index  n2),  and  denoi.  q.c  h,  ,.i  i  .q/i ,,  , 
functions  in  these  planes  by  subscripts  -  and  f  ,  respei  u  .  iv  i  u  i  iqoii  c  1;  AlaM.,; 
use  of  the  facts  that  there  arc  no  back  ward -propagating  |ie|,|.,  j , .  (|i ,  (joii  p,  m  gh  .  u  qj 
and  that  the  cross-spectral  density  function  If  Ip, ,  p2 )  rent. on  , .  out  unions  .a  toss  t  (■.- 
boundary,  we  then  obtain,  in  the  limit  as  i:  ~»(),  from  equation.  (  j)  p,  qu  t,  |a|j,.,,' 

►  ,  .  '  ;  I 

/<  .  (p.S, .  )  =  (w:/h,  )i.3/(p;  S,  t  ,Sj ,  )H  O, 

where 


3/(p;S, ,  ,S2j )  = 


cos  02  jpp.  k , S 

c‘>»  /(#*,  A  t«, , 


f 


<  I 
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In  this  expression,  =  nxk 0  and  k2  =  n2k 0  with  k0  =  oj/c,  as  before,  and  0t  and  02  are 
the  angles  that  the  directions  St  and  S2  make  with  the  positive  z  axis. 

Formula  (6)  may  now  be  regarded  as  a  generalization  (at  a  planar  boundary)  of  the 
conventional  radiance  theorem  (1 )  for  optical  wavehelds  of  any  state  of  coherence.  It 
is,  in  essence,  an  identity  that  follows  directly  from  the  basic  assumptions  of  physical 
optics.  Comparison  of  equations  (1)  and  (6)  reveals  that  the  new  relation  (6)  contains 
an  additional  factor  M{p\  S,  2 ,  S2l),  which  is  given  by  equation  (7).  This  factor, 
determined  primarily  by  the  state  of  coherence  of  the  optical  field  at  the  boundary,  is 
a  measure  of  the  extent  to  which  the  present  generalized  result  differs  from  the 
conventional  law  connecting  the  radiances  B,  and  B2 .  Through  the  cross-spectral 
density  appearing  in  the  definition  (4)  of  the  Wigner  function  y,  the  factor  M 
depends,  in  general,  both  on  the  optical  intensity  and  on  the  complex  degree  of 
spatial  coherence  of  the  light  at  the  interface.  In  fact,  we  see  from  equation  (7)  that, 
apart  from  a  purely  geometrical  part,  the  factor  M  is  simply  the  ratio  of  the  values 
that  the  associated  function  y(p,  f)  assumes  with  the  arguments  f2  =  fc2S2i  and 
f  i  =  ^1*1 1 • 

If  the  directional  vectors  and  s2  specify  a  geometrical  ray  path  across  the 
boundary,  then  k2S21  =  according  to  Snell’s  law  and  the  additional  factor  M 

reduces  in  this  case  to  the  ratio  cos#2/cos0, .  This  result  is  a  consequence  of  the 
requirements  that  there  is  no  reflected  wave  corresponding  to  a  wave  incident  from 
the  direction  S,  and  that  the  transmitted  wave  in  the  direction  S2  matches  the  values 
of  the  incident  wave  at  the  interface.  The  result  holds  separately  for  any  incident- 
wave  direction  and  also  implies  that  under  the  present  assumptions  energy  is  not 
strictly  conserved  in  passage  across  the  surface.  However,  if  the  appropriate 
reflection  and  transmission  coefficients  are  included,  the  energy  conservation  for 
plane  waves  is  restored.  Conversely,  straightforward  calculations  using  the  general 
results  ((>)  and  (7)  show  that  if  II  {p,  S, ;)  is  zero  except  for  some  value  s„  ,  then 
H  +  {p,  S2i)  will  also  differ  from  zero  only  when  n2 S21  =  >/,S0i,  in  accordance  w  ith 
geometrical  optics. 

We  will  emphasize,  furthermore,  that  unlike  in  the  conventional  radiance 
theorem  (1),  the  variables  S12  and  S2i  in  the  generalized  result  expressed  by 
equations  (6)  and  (7)  are  projections  of  quite  arbitrary  unit  vectors  that  point  towards 
the  half-space  z  >  0.  This  makes  it  possible  to  use  the  analytic  properties  of  the 
generalized  radiance  function  R(p,  Sx).  It  implies  also,  for  example,  that  in  the  case 
when  w,  <«2,  the  generalized  radiance  B+(p,  S2i)  may  be  non-zero  even  in  the 
domain  w,/n2  <  |S2i|  <  1 ,  corresponding  to  angles  02  larger  than  the  critical  angle  of 
total  internal  reflection.  Physically,  such  a  situation  represents  the  phenomenon 
where  evanescent  wav  s  are  turned  into  homogeneous  (propagating)  waves  by 
refraction  at  the  discontinuity. 


4.  Radiance  theorem  with  quasi-homogeneous  light 

Let  us  assume  now  that  the  optical  field  at  the  interface  is  quasi-homogeneous, 
i.e.  one  that  is  characterized  by  a  cross-spectral  density  function  of  the  form  [24] 

H’<Pi,p2)  =  /((Pi  +Pi)/2)g(pi  - p2 ),  (S) 

where  /(p),  the  optical  intensity,  is  a  'slow'  function  of  p  and  g(p'),  the  complex 
degree  of  spatial  coherence  [18],  is  a  'fast’  function  of  p'  (see  [24],  §  1 1).  The  Wigner 

( 
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function  associated  with  such  a  field  distribution  is  readily  found  from  equation  (4) 
to  be  given  by  the  expression 

/■(/».  f )  =  (2n)2l(p)g(i ),  (9) 

where 


c'ffi-d  2tt)  | nip' ) exp (  if  •  p)  d  ’ ft  (Mi 

is  the  two-dimensional  spatial  Fourier  transform  of  g(p).  On  substituting  from 
equation  (9)  into  the  general  formula  (7),  we  obtain  the  expression 


'/(,>;  S,  ,  ,  S2l 


cos  02g{k2 S2l  ) 
cos  0 1  g(/t  i  S  |  ^  ) 


(11) 


Fquation  (II)  show  s  that  for  a  quasi-homogeneous  field  the  factor  A/,  which  is  absent 
in  the  usual  formula  (1),  is  independent  of  the  optical  intensity  of  the  light 
distribution  at  the  interface.  Moreover,  since g(p')  is  a  ‘fast’  function  of  p',  its  Fourier 
transform  g(f )  is  a  ‘slow  ’  function  of  f.  C  onsequently  for  quasi -homogeneous  light 
\I{p,Sit,S2L)  ^  \  and  the  generalized  result  (6)  is  seen  to  approximate  the 
conventional  radiance  theorem  (I)  with  relatively  good  accuracy  over  a  range  of 
directions  S,  and  S2  such  that  |S ,  ,  |  5:  |S2 t 1.  We  note  briefly  also  that  for  statistically 
homogeneous  fields  equation  (11)  remains  valid  even  when  g(p')  is  not  a  sharply- 
peaked  function.  For  such  fields  the  generalized  result  (6)  is  seen  to  depend  on  the 
functional  form  of  the  complex  degree  of  spatial  coherence  g(p  ). 

Let  us  now  assume,  furthermore,  that  n{  <  tii  and  that  the  complex  degree  of 
spatial  coherence  of  the  light  at  the  interface  is  given  by  the  expression 


Jf(p')  = 


sin  fe,p' 
ktp' 


(12) 


where  p  -  )p'|.  This  expression  is  characteristic  of  a  Lambertian  radiator,  such  as 
hlackbody  radiation  source  [25].  Making  use  of  the  definition  (2),  the  generalized 
radiance  tt  (p,  S, ,  )  is  then  found  to  be  independent  of  the  directional  variable  S,  2  , 
and  we  will  denote  it  by  /f,,(p).  On  substituting  from  equation  (12)  into  equation  (11) 
and  making  use  of  formula  ((>),  we  obtain  the  following  result: 


/  7 

/ Y 

C‘*S?02  1 

\»j 

1  -  («2  )‘  sin*  02  J 

H  ,  (p,  S21  )  = 


-  0 

I  lere  the  angle  ()t  is  defined  In  the  relation 


if  0,  >  0... 


(13) 


sin  ()c  =  w, ,  /), . 


(14) 


The  formula  (15)  shows  that  it  the  field  at  the  interface  is  quasi-homogcncous  (or 
strutlv  homogeneous)  with  its  complex  degree  of  spatial  coherence  given  by 
equation  (12).  then  there  is  a  maximum  angle,  0C,  beyond  which  no  energy  is 
l ransmitted.t  According  to  equation  (14/,  this  angle  is  precisely  the  critical  angle  of 
total  internal  reflection 


f  This  result  is  a  consequence  of  the  fact  that  the  spatial  Fourier  transform  of  the 
ci >r relation  turn  non  (  I  2 )  is  identically  zero  outside  the  domain  |f  |  <  /t,  (see  [25],  tj  1 1 ).  I  fence 
the  field  iiii  idem  mi  the  boundary  contains  no  evanescent  waves  that  could  he  turned  In  the 
disi  ontimiitv  siirtaee  tnlo  homogeneous  waves  propagating  at  angles  larger  than  (), 
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Formula  (13)  also  implies  that  even  lor  angles  less  than  the  maximum  angle  (h  , 
there  is  strictly  speaking  an  angular  dependence  that  is  not  present  in  the 
conventional  form  (1)  of  the  radiance  theorem.  This  angular  dependence  is  a 
consequence  of  the  basic  assumptions  of  physical  optics  (compare  the  discussion 
following  equations  (6)  and  (7)),  and  it  is  illustrated  in  figure  3  for  two  values  of  the 
ratio  n 2  n,  .  These  curves,  calculated  according  to  equation  (13),  indicate  that  the 
ratio  R ,  (p,  s2 1 ):  R„(p)  remains  substantially  constant  over  a  relatively  wide  range  of 
angles  tl2,  in  agreement  with  the  conventional  radiance  theorem.  Such  a  behaviour 
becomes  even  more  dominant  as  the  ratio  «,  //,  is  decreased  In  the  limit  as  n,  nt 
approaches  unity,  there  is  no  retracting  surface  and  R ,  (p,S,  ,  )  becomes,  of  course, 
identical  to  R„{p) 


BJ?.sZL){B0l 9) 


I  'Cure  .V  I fependem e  ot  II ,  (/*,  s2 .  )  ltn(p)  on  the  angle  0 ,  (sin  ll2  -  |s, ,  |)  for  two  values  of 
the  ratio  n,  «, .  namely  I  I  and  15,  when  the  degree  ol  field  correlation  at  the  interface 
is  given  by  equation  (12).  The  angle  (Jc  denotes  the  maximum  angle  beyond  which  no 
energy  is  transmitted. 

5.  Summary  and  discussion 

I  n  this  paper  we  studied  the  radiance  theorem  in  the  context  of  partially  coherent 
waves  and  considered  only  retraction  at  a  planar  interface  separating  two  homo¬ 
geneous  media.  The  analysis  was  carried  out  within  the  framework  of  the  scalar 
theory  of  light.  It  was  based  on  the  assumptions  that  the  optical  field  remains 
continuous  across  the  boundary  and  that,  as  is  customary  in  physical  optics,  the 
effects  ot  reflection  can  be  neglected.  The  discontinuity  may  therefore  be  thought  of 
merely  as  a  limiting  case  ol  an  optical  element  represented  by  an  amplitude 
transmission  function  Up),  with  Up)  approaching  unity.  Since  the  transmission 
function  is  independent  ol  the  properties  of  tin-  incident  field,  such  as  its  direction  of 
propagation,  this  method  typically  leads  to  results  that  can  be  expected  to  hold  only 
in  the  paraxial  regime. 

Our  analysis  showed  that  the  radiance  theorem  with  light  of  any  state  of 
coherence  contains  an  additional  factor,  not  present  in  the  conventional  radiance 
theorem,  that  depends  in  general  both  on  the  optical  intensity  and  on  the  complex 
degree  ol  spatial  coherence  ol  the  light  at  the  interlace.  For  a  qiiasi-homogeneous 
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held,  tills  factor  bet nines  imit-pcn  lent  nt  tin-  nptiu;  ir.t  -ii'.i’ .  elisfrihu'  mi  mid 
HIUII,  apart  from  .1  |4<  oillcti  Hill  purl  ms!i;  in\  li,  ,  li'  ,;n  :  1  •  i  r  ■  e  tht  ■  .1  • 

transform  of  tin-  uiin|iK'\  dcijrcc  of  spatial  t  ohen  me  n  il  t  ■  1  •  ,u  \<  ■  •  r 

f,  /f  ,Sj  and  f,  =  /•-■ , S ,  .  rcspcclivcK  (see  ci)u.ltli,;  (  '.  I  •  i,  .  .•<••• 

dearie  that  tin -cm-ige  transmission  across  a  medium  u.  •  "'1  1 

till  till'  state'  ot  ciilici'i'in  e  1  >1  th;-  ua\  t  helti 

It  was  also  shown  that,  foi  bi  le  kbody  raeliation  lie  kb.  :),  .  .,r  •  > 

radiance  theorem  lhat  we  obtained  is  in  rcl.il  i\ civ  . .  1  ,i>  • 

aim  etitional  radians  e  theorem  <w  nr  a  wide  rat  up  o{  ,11 je~.  w  i  |  {1 

h^lit  m. i\  he  assumed  to  he-  .1  reasonable  approe  imut,  0  \u  u  u  ■  .  I  .  I 

theorem  with  partially  e  oliei  ent  Intht  at  l.um  1  angles  oj  ,  ,,  ,,|.  1  .  .  ,i>  ,  .  . 

use  ot  the-  lull  electromagnet n  theory,  with  propel  .  ■  i.  ..  1  ,1  , 

conditions  at  the  dist  out  mint  \  mu  hue 
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The  question  is  raised  as  to  whether  the  normalized  spectrum  of  light  remains  unchanged  on 
propagation  through  free  space  It  is  shown  that  for  sources  of  a  certain  class  that  includes  the  usu¬ 
al  thermal  sources,  the  normalized  spectrum  will,  in  general,  depend  on  the  location  of  the  observa 
non  pomi  unless  the  degree  of  spectral  coherence  of  the  light  across  the  source  obeys  a  certain  scal¬ 
ing  law  Possible  implications  of  the  analysts  for  astrophysics  are  mentioned 
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Measurements  of  the  spectrum  of  light  are  generally 
made  some  distance  away  from  its  sources  and  in 
mans  cases,  as  for  example  in  astronomy,  they  are 
made  exceedingly  far  away  It  is  taken  for  granted  that 
the  normalized  spectral  distribution  of  the  light  in¬ 
cident  on  a  detector  after  propagation  from  the  source 
through  free  space  is  the  same  as  that  of  the  light  in 
the  source  region  I  will  refer  to  this  assumption  as 
the  assumption  of  tnvananie  oj  the  spectrum  on  prop¬ 
agation  This  assumption,  which  is  implicit  in  all  of 
spectroscopy,  does  not  appear  to  have  been  previously 
questioned,  probably  because  with  light  from  tradition¬ 
al  sources  one  has  never  encountered  any  problems 
with  n  However,  with  the  gradual  development  of 
rather  unconventional  light  sources  and  with  the  rela¬ 
tively  frequent  discoveries  of  stellar  objects  of  an  un¬ 
familiar  kind,  it  is  obviously  desirable  to  understand 
whether  all  such  sources  generate  light  whose  Spec¬ 
trum  is  invariant  on  propagation,  and  if  so,  what  the 
reasons  for  n  are  Actually  u  is  not  difficult  to  con¬ 
ceive  of  sources  that  generate  light  whose  spectrum  is 
not  invariant  on  propagation  In  this  note  I  will  shew 
what  are  the  characteristics  of  a  certain  class  of  sources 
that  generate  light  whose  spectrum  is  invariant,  at  least 
in  the  far  zone 

From  the  standpoint  of  optical  coherence  theory,  in¬ 
variance  of  the  spectrum  of  light  on  propagation  from 
conventional  sources  is  a  rather  remarkable  fact,  as  can 
be  seen  from  the  following  simple  argument  Consid¬ 
er  an  optical  field  generated  by  a  stationary  source  in 
free  space  The  basic  field  variable,  say  the  electric 
field  strength  at  the  space-time  point  (r.r),  nay  be 
represented  by  its  complex  analytic  signal1  J  £(r,f). 
According  to  the  Wiener-Khintchine  theorem3  the 
spectral  density  of  the  light  at  the  point  r  is  then 
represented  by  the  Fourier  transform. 


S  (  r.  u,  )  < 


I  <r.  rle'-’dr. 


of  the  autocorrelation  function  (known  in  the  optical 
context  as  the  self-coherence  function)  of  the  field 
variable  It  is  defined  as 

Mr  r  >  =  <  ft  r.r  >f  ( r,r  r  )  > .  (2) 


where  the  angular  brackets  denote  the  ensemble  aver¬ 
age.  Now  the  spectral  density  and  the  self-coherence 
function  are  the  "diagonal  elements”  ( r2  —  r1  =  r )  of 
two  basic  optical  correlation  functions,  viz.,  the  cross- 
spectral  density 

W(r,,r2.<p)-  J7„^(^|.r2.T)e,•,,  dr.  (3) 

and  the  mutual  coherence  function 

r(r|t  r2.  r )  -  (E’(rt.t  )£(r2,/  +  t  )) .  (4) 

It  is  well  known  that  both  the  mutual  coherence  func 
tion  and  the  cross-spectral  density  obey  precise  prop 
agation  laws  For  example,  in  free  space4 

(V/  +  AJ)M  <r,.r2.ru)  =  0  (y  -  1.2).  (5) 


k  —  ui/c.  (6) 

with  c  being  the  speed  of  light  m  vacuo  and  V;  being 
the  Laplacian  operator  acting  with  respect  to  the  van 
able  Tj.  Consequently,  both  the  mutual  coherence 
function  and  the  cross-spectral  density  and.  in  fact, 
also  their  normalized  values  change  appreciably  on 
propagation  For  example,  for  a  spatially  incoherent 
planar  source  W ( rt . r2. w )  and  r(r|,r2,  t)  will  be 
essentially  6  correlated  with  respect  to  t,  and  r2  at  the 
source  plane  but  will  have  nonzero  values  for  widely 
separated  pairs  of  points  which  are  sufficiently  far  away 
from  the  source  This  is  the  essence  of  the  well 
known  van  Cittert-Zernike  theorem  (Ref  1,  Sect 
10.4.2).  In  physical  terms,  the  correlation  in  the  field 
generated  by  a  spatially  incoherent  source  may  be 
shown  to  have  its  origin  in  the  process  of  superposi¬ 
tion  We  thus  have  the  following  rather  strange  situa¬ 
tion:  The  correlations  of  the  light  may  change  drasti¬ 
cally  on  propagation;  yet,  under  commonly  occurring 
circumstances,  their  (suitably  normalized)  diagonal 
elements,  which  represent  the  spectrum  of  the  light  or 
its  Fourier  transform,  remain  unchanged 
To  obtain  some  insight  into  this  problem  we  con 
sider  light  generated  by  a  very  simple  mode!  source, 
namely,  a  planar  source  occupying  a  finite  domain  I)  of 
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a  plane  r  «=  0  and  radiating  into  the  half  space  r  >  0. 
which  has  the  same  spectral  distribution  S"”(cu)  at 
each  source  point  Pip )  and  whose  degree  of  spectral 
coherence'  p(0'(pi. P:  <"  >  is  statisticall>  homogene¬ 
ous.  i.e..  has  the  functional  form  #i<Cl,(p:-p,.a») 
The  cross-spectral  density  of  the  light  across  the 
source  plane  is  then  given  by 

^,0,(pi.p2.Oi) 

-<(p,  )<(pj)Sl®1(«i>)/j<0l(p2  —  pi.«o),  (7) 

where  «(p)=  1  or  0  according  to  whether  the  point 
Pip)  is  located  within  or  outside  the  source  area  D  in 
the  plane  r  —  0 

We  will  also  assume  that  at  each  effective  frequency 
<«  present  in  the  source  spectrum,  the  linear  dimen¬ 
sions  of  the  source  are  much  larger  than  the  spectral 
correlation  length  [the  effective  width  A  of  lpl0’(p 
iu)|)  Sources  of  this  kind  belong  to  the  class  of  so- 
called  quasihomogeneous  sources  *  which  have  been  ex¬ 
tensively  studied  in  coherence  theory  in  recent  years 
Most  of  the  usual  thermal  sources  are  of  this  kind 
The  radiant  intensity  i  e  .  the  rate  at  which 

energy  is  radiated  at  frequency  w  per  unit  solid  angle 
around  a  direction  specified  by  a  unit  vector  u,  is  given 
by  the  expression  [cf  Ref  6.  Eq  (48)1 

JJ  u>  «=  k  i  AS'0' (uj  )p'"'U  u.,o. )  cos*  ft  (8) 

In  this  formula,  A  is  the  area  of  the  source, 

u  t  f .  cu )  = - - — ,  (>■<„. «>.  f  p  d‘  u  (9) 

<2n)‘J 

is  the  two-dimensional  spatial  Fourier  transform  of  the 
degree  of  spectral  coherence,  u.  is  the  transverse  part 
of  ihe  unit  vector  u.  i  e  .  the  component  of  u  (con¬ 
sidered  as  a  two-dimensional  vector)  perpendicular  to 
the  r  axis,  and  ft  is  the  angle  between  the  u  and  the  : 
directions  (see  Fig  I)  Evidently  the  normalized  spec¬ 
tral  density  S'”'(u.a/>  at  a  point  in  the  far  zone,  in 
the  direction  specified  by  the  unit  vector  u,  is  given  by 

S'" '( u.wl.yju  )//•/„<  uWcu  (10) 


FKi  I  lllusiranon  of  die  notation 


On  substituting  Eq.  (8)  into  Eq  (10)  we  obtain  for  the 
normalized  spectrum  in  the  far  zone  the  expression 

k  o  ( Of  )  u  (Au.loi 

S  “’(u.w)-- - - - -  till 

j  S1®1  (<i>  )  p  ik\l.m)diu 

It  is  clear  from  Eq  (11)  that  the  normalized  spc\- 
trum  of  the  light  depends  on  the  direction  u.  i.e  .  it  is 
in  general  not  invariant  throughout  the  fur  zone 
However,  it  is  seen  at  once  from  Eq  (11)  that  it  will 
be  invariant  throughout  the  far  zone  if  the  Fourier 
transform  of  the  degree  of  spectral  coherence  of  the 
light  in  the  source  plane  is  the  product  of  a  function  ol 
frequency  and  a  function  of  direction,  i.e  .  it  is  of  the 
form 

^“"Uuj.wl-FUJ/'/tu.  )  tl.'i 

In  this  case  Eq  (11)  reduces  to 

— 'fLL  .  *J5,<"(  0,)FU> 

J  l  U .  U)  )  —jT-  - -  .  II  .'- 

j  •’  ' (w  )  F  (w  I  dw 

and  the  expression  on  the  right  is  independent  of  the 
dtrecuon  u 

I  will  now  show  that  the  condition  (12)  has  some  in¬ 
teresting  implications,  which  follow  from  the  lad  th 
pl0>  is  a  correlation  coefficient.  Before  doing  this  w-„ 
note  that  since  u  is  a  unit  vector.  iu_  ;  <  1  Howese* 
we  will  now  assume  that  the  factorization  cor.ditn  • 
(12)  holds  for  all  two-dimensional  sector-  u  ’’ 
*S  |ua  I  <  oo ).  This  assumption  will  be  trivially  s.-.t - 
Tied  if  the  degree  of  spectral  coherence  p"  (p  .<-  1  i- 
at  each  effective  temporal  frequency  w.  band  limits  c 
in  the  spatial  frequency  plane  to  a  circle  ol  rudiu-  - 
about  the  origin,  in  more  physical  terms  this  conditi- 
means  that  p<0,(p  .  <u  )  docs  not  vary  appres  ubls  n\ . 
distances  of  the  order  of  the  wavelength  a  -  2-  . 

With  this  being  understood  let  us  take  the  lou*  d 
transform  of  Eq.  (12).  We  then  find  at  once  m.i* 

pl0>(p'.<u ) 

~F(w)J'/y(ui)exp(/AuJ  p  )  d'ik  u  i  < !  c 

i.e., 

p'0,(p'.cu)  =  A*7u)/yup>, 

where  H  is,  of  course,  the  two-dimensional  Found 
transform  of  H  Since  pl0,(p'.w)  is  a  correlation  u 
efficient  it  has  the  value  unity  when  p  «=  0,  ic  . 

w )  “  1 ,  for  all  <u.  i)t 

and  hence  Eq  (15)  implies  that 

A 2  7  (  oj  )  —  1/7(0))-'  <  r  > 

Since  Ihe  left-hand  side  of  Eq  (I7i  depends  or.  uv. 
frequency  but  the  right-hand  side  is  independent  of  it. 
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each  side  must  be  a  constant  (o  say)  and  consequently 

ft  (18) 

7 wo  important  conclusions  follow  at  once  from 
inese  results  If  we  substitute  Eq  (18)  into  Eq  (13) 
hl  obtain  the  following  expression  for  the  normalized 
spectrum  of  light  in  the  far  zone 

S  ~  (u  u,  -y-S  ~(—  —  .  (19) 

J  S(0’(ii>)d«o 

7  his  formula  shows  that  not  only  is  the  normalized 
spectrum  of  the  light  now  the  same  throughout  the  far 
zone,  but  it  is  also  equal  to  the  normalized  spectrum  of 
;he  light  at  each  source  point 

Next  we  substitute  Eq  (I8i  into  Eq  (15)  and  set 
uH  «  h.  p  =■  p;  -  f,t  We  then  obtain  for  p101  the  ex¬ 
pression 

|i'  =  —  p ,)} 

U-Wc);  (20) 

i  e  the  complex  degree  of  spectral  coherence  is  a 
I  unction  of  the  variable  £  =  klp;-p,)  only  We  will 
reter  to  Eq  (20)  as  the  si almy  Ian  Obviously  for  a 
source  that  satisfies  this  law.  the  knowledge  of  the  de¬ 
gree  of  spectral  coherence  of  the  light  in  the  source 
plane  ai  ani  particular  frequency  cu  specifies  it  for  all 
Irequenues 

1  he  scaling  law  (20).  which  ensures  that  for  sources 
hi  the  class  that  we  are  considering  the  normalized 
spectrum  of  the  light  is  the  same  throughout  the  far 
zone  and  is  equal  to  the  normalized  spectrum  of  the 
iight  j(  ea^h  source  poini  [Eq  1 19)].  is  the  mam  result 
nl  inis  note 

li  is  natural  to  inquire  whether  sources  are  known 
thai  obec  this  scaling  law  The  answer  is  affirmative. 
Many  oi  tne  commonly  occurring  sources,  including 
NackhoJy  sources,  obey  Lambert's  radiation  law  [Ref 
1  Sect  4  8  1)  It  is  known  that  all  quasi-homo- 
geneous  Lambertian  sources  have  the  same  degree  of 
sp-cciral  coherence,  viz 

p"  *  P:  ~  Pi- «**  *  “  Sin <  k  \p2  ~  pi  I  )!k  \p2  ~  P 1 1 .  (21) 

which  is  seen  to  satisfy  the  scaling  law  (20).  Accord¬ 
ing  to  the  preceding  analysis  such  sources  will  generate 
light  whose  normalized  spectrum  is  the  same 
throughout  the  far  zone  and  is  equal  to  the  normalized 
spectrum  at  each  source  point  This  fact  is  undoubted¬ 


ly  largely  responsible  for  the  commonly  held,  but 
nevertheless  incorrect,  belief  that  spectral  invariance  is 
a  general  property  of  light. 

This  Letter  has  dealt  with  what  is  probably  the  sim¬ 
plest  problem  regarding  spectral  invariance  on  prop¬ 
agation  It  would  seem  that  some  significant  questions 
in  this  area  might  be  profitably  studied  Among  them 
are  the  elucidation  of  the  physical  origin  of  the  scaling 
law,  spectral  properties  of  light  from  a  broader  class  of 
sources  than  considered  here,  the  relation  between  the 
scaling  law  and  Mandel's  results  regarding  cross- 
spectrally  pure  light,8  9  and  relativistic  effects  Appli¬ 
cations  of  the  results  to  problems  of  astrophysics 
might  be  of  particular  interest;  at  this  stage  one  might 
only  speculate  whether  source  correlations  may 
perhaps  not  give  rise  to  differences  between  the  spec¬ 
trum  of  the  emitted  light  and  the  spectrum  of  the 
detected  light  that  originates  in  some  stellar  sources 
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It  is  shown  that  for  fields  produced  by  quasi-homogeneous  sources,  one  of  the  generalized  radiance  functions  introduced  by 
Walther  reduces,  in  the  asymptotic  limit  of  large  wave  number,  to  a  function  that  has  all  the  basic  properties  of  the  radiance  of 
traditional  radiometrv  An  explicit  expression  for  this  radiance  is  obtained  in  terms  of  the  intensity  distribution  across  the 
source  and  the  degree  of  spectral  coherence  of  the  source.  The  results  provide  a  rigorous  foundation  for  radiometrv  in  free 
space,  on  the  basis  of  statistical  wave  theory 
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It  n  th.it  lot  fields  piodmcd  hi  quasi  homogeneous  sources.  one  of  the  genet  .1 1 1/cit  taili.ittie  (uiictuuis  introiluieil  hi 

\X  alt  her  teihiees.  in  the  -isi  inpioln  limit  nl  large  wave  iiumhet .  h*  a  I  unction  that  It.is  all  the  ham  piupemes  of  the  radiance  al 
tiaihtiiiti.il  r ailu uiiet i \  Ait  esplisit  esptessnm  lut  litis  lathanee  is  obtained  in  teuus  nl  the  intensiti  ihsnihtitinn  airnss  the 
■.untie  .nisi  the  degree  ol  speelr.il  coherence  ul  the  source  1  he  lesults  pi  ovule  a  rigoious  foumJafioii  lor  i.ulnuiieiri  in  lice 
space,  nn  the  hasu  nl  statistical  nuie  theors 


V. 

>, 

V. 

V 

>. 


1 .  Introduction 

During  (tie  last  two  decades  several  attempts  have 
been  made  to  elucidate  the  Inundations  of  radiom¬ 
ens  In  particular  several  authors  [  1  5]  proposed 
expressions  for  the  basic  quantity  of  radiometry, 
namely  the  (spectral)  radiance,  in  terms  of  various 
second-order  correlation  functions  of  the  optical 
field.  Allhough  each  of  the  proposed  expressions  ex- 
hibiis  some  ol  the  well-known  properties  that  are  al- 
trihuicd  to  the  radiance  in  traditional  radiometry, 
none  of  them  possesses  all  of  them,  for  sources  and 
fields  of  arbitrary  slate  of  coherence.  In  particular, 
some  of  the  proposed  expressions  for  the  radiance 
can  take  on  negative  values,  a  result  that  contradicts 
the  physical  meaning  of  radiance.  More  recently  it 
was  shown  |6|  that  it  is  not  possible  to  define  a  radi¬ 
ance  for  a  planar  source  which  depends  linearly  on  a 
second-order  correlation  function  of  the  source  field 
and  which  satisfies  three  basic  postulates  of  radionr- 
etry  for  every  possible  state  of  coherence  of  the 
source  1 1 

Kose.irrh  Mipporicil  b\  the  National  Soeiuc  I  iniml.iiu'n 
under  <jf.if»r  s?  PHY-X  and  the  Ah  I  orcc  (jcophssuN 

l  ahoratorv  under  Al  OSR  I a\k 
1  On  le.i\c  from  the  Department  ol  PhvMw  Mississippi  State 
I  mveiMH.  Mississippi  State.  MS  W7ft2.  I  SA. 

Also  at  the  Institute  »>f  Optics.  \  mversitv  of  Rochester 


We  believe  that  the  difficulties  just  mentioned 
arose  because  the  previous  investigations  did  not  take 
into  account  the  fact  that  traditional  radiometry  deals 
with  sources  that  are  spatially  highly  incoherent 
(namely  thermal  sources)  and  that  they  generate  radia¬ 
tion  whose  effective  wavelengths  X  are  very  small 
compared  with  their  linear  dimensions  1 2 .  We  show 
in  this  note  that  when  these  facts  arc  taken  into  ac¬ 
count  a  consistent  formulation  of  radiometry  is  ob¬ 
tained  on  the  basis  of  second-order  cohcren.c  theory, 
at  least  for  sources  and  fields  in  free  space.  ..lore  spec¬ 
ifically,  we  show  that  traditional  radiometry  correctly 
describes  the  behavior  ot  fields  generated  by  planar 
quasi-hotnogeneous  sources  1 1 1]  in  free  space,  in  the 
asymptotic  limit  as  the  wave  number  A  -  2n/X  -*  ®. 


"  In  an  interesting  recent  paper  1 7  ]  a  definition  of  radiance 
was  proposed  which  depends  non-lincarly  on  a  second- 
order  correlation  I  unction  ul  the  source  and  which  satis¬ 
fies  the  three  postulates.  It  appears,  however,  that  this 
radiance  does  not  obey  the  radiometric  law  lor  (lie  propa- 
palion  ol  radiance  m  tree  space. 

12  Allusions  lo  the  possibility  ih.il  traditional  radiometry 
implies  such  restriilions  have  teen  made  from  time  to 
time  |4h|.  |K  I0|.  hul  tile  appropriate  mathematical 
lustificatmn  has  not  he  n  previously  provided. 
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2.  Generalized  radiance 

1  ci  us  consider  a  second. irv  so.iree  n,  occupyinL'  a 
finite  portion  o!  the  plane  :  il  and  radiaiinp  into  ilie 
hall-space  r  >0.  We  assume  that  the  source  fluctua¬ 
tions  are  statistically  stationary.  We  will  denote  hy  P 
the  two-dimensional  vector  specilvin.'  the  location  ol 
a  source  point  S  and  hy  r  the  three-dimensional  veetoi 
speeit  vine  the  location  of  a  field  point  !’  m  the  hall  - 
space  :  2-  0.  hotli  referred  to  a  lixed  orieiu  0  in  the 
s'  mice  ration  ( see  tip.  I  ). 

let  h’lt'i .  r,.  t  )  he  the  eross-speettal  density  ol  the 
Held  venerated  In  the  somce  at  two  |ioints  I’ j  and  I’,. 

1 1  is  know  n  l  ha  i  the  cioss-spei  1 1  a  I  densitv  may  he 
n'piesenled  in  I o i  ms  nl  an  ensemble  ol  monochrom.il- 
ic  wave  liehls  f  t>\r.  t<)  e\p(  )  '■ .  all  ol  the  same 

t  rec|uenev  r.  as  |  I  2 1 

h’h|.  /•,.  r)  -  <  ( <  *(/-|.  »■>  (■(/■>.  i’)>.  (-.11 

■vheie  the  anunlar  bracket  on  the  iieht-hand  side  ol 
ei|.  i  2  1 1  denotes  the  aveiape  over  this  ensemble.  I  he 
spate  dependent  pail  of  each  member  ol  this  ensemble 
.  ihevs.  throimhout  the  hall -space  c  . *  0,  the  Helmholtz 
equation 

(7-  ♦  kd  fir.  p)  =  0.  (2.2) 

w hem  k  -  .2 .7 /  o  and  i  is  the  speeu  ol  lielit  in  vacuo,  and 
il  behaves  as  an  oiiteoimr  sirhetie.il  wave  al  intimiv  in 
diis  half-space  \s  is  well  know n.  sneh  solnlions  can 
Ire  represented,  under  very  i\  neral  conditions,  in  the 
form  of, in  aneular  sped  mm  ol  plane  waves,  i.e.  in  the 
Hum  |  I  ''I 

I  (r .  r)  j  </|  v  i)  espt  i  ks  r )  d’.vj .  (  2.2  ) 


where 


*  "  (SX.S,..  )- 

i,  1  (v,  .  v,..  0). 

(2.4) 

f  4 

-r 

(i 

it  's  -1. 

-  +«sf  i»l/2 

if  .1,  !  -1. 

( 2  ' ) 

and  the  intemalion 

on  the  iivlil-haiid  side  ol  eq. 

(27 

extends  over  the  whole  .s  plane  (O'  v;  r  v|  ') 

In  the  above  notation  the  complex  veision  ol  the 
generalized  radiance  function  mtioduccd  In  Walthei 
in  refs.  [4|  may  he  written  in  the  Imm 

si  i')a(s, .  rp  exp(ik.v  r |  (2.(0 

I  oi  pm  :>o  e ,  o :  lalei  ili  si  nssion  we  will  ex  pi  ess 
H  v )  m  lei  lie.  "i  1 1-  ■  values  It  J"'l  p .  v ).  w  I  in  h  il 
takes  on  (lie -.oui,  e  |  .ate  c  o  I  his  ,an  leadtlv  he 
done  iry  making  use  ol  the  tael  that  the  outgoing  solu¬ 
tion  of  the  I  lelnihol:/  :  ■  n ,  may  he  exptessed  in 
teims  of  its  lioiin  !  n\  vain  .  '  p  /’)  in  the  plane 

:  =  ()  In  the  Rayleigh  tormula  |  I  a, 

//(r.  r)  -  i  )  l p. ;  ■)  d-p.  (27) 

n 

wheie  R  ~  r  p  and  (H  R.i>)  is  the  (.teen  s  'unci ;■  m 

(,'{R  i  )  ( 1 ,2n )(<'i/dc)|exp(t^A’ I  R\.  (271 

with  R  i R  (In  suirslitutmi'  loi  /  (r.  r)  lioiti  eq. 

( 27)  into  eq  ( 2.<>)  and  mteicliangmg  the  ordeix  ot 
integration  and  avciaging  one  readily  finds  that 

lf,.(r.  si  expli/t.vr)  f(!*(R.i) 

n 

X  HI'.",(  p.  .vlexp(  itj,  •  p)  d’/j.  (2.0) 

where  il  |l.",(p.s|  is  the  net ierali/cd  t.idtance  at  a 
)>•  lint  in  the  somee  plane,  viz .. 

|(  [n,(  p .  s)  “•  \ .  I 'l0|d  p,  i>)  atsi ,  r)i  exp(  iks,  p  ) 

(2.101 


P- 


I  ' 


I  i:  1 .  Kliisii.ilinr  lire  not. emu 


In  general  |1l(r..vlix  complex,  lienee  it  cannot 
have  the  phvxtcal  significance  ol  i.uliancc.  I  he  same 
is  1 1  tie  ol  tlic  teal  pat  t  ol  H  (r.  v  I.  hec.mse.  as  w  as 
show  n  el  sew  licit'  |  I  s  |  .  il  can  soi  net  lines  take  on  in  fa¬ 
me  \  a  hie  .  |  low  evei ,  as  show  n  in  tel  |-),i  |  it  s  ie  il 
p.u  l  V  iclds  l  he  i .  meet  value  hn  the  i  ad  laid  in  ten  -  it  \ 
via  one  ol  tin  sl.ind.nd  loimtilas  ol  tiadition.il  t .idiom 
etiv 
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3.  Generalized  radiance  of  a  field  generated  by  a 
qiiasi-homogcncuus  source 

He  ml!  now  specialize  llic  expressions  (2.  10)  and 
(2  '>)  loi  the  genciali/ed  radiance  to  the  case  when 
the  source  is  qiiasi-lioinogeneous  1 1  I J.  flic  degree  of 
spec t tal  coherence  ol  a  quusi-homogcncous  source 
ilepends  on  the  two  source  variables  p|.p->on!v 
through  the  ditlerencc  P->  Pj.  Consequently  its 
cross-spectral  density  has  the  form 

l»’"”(P,.  P2,p) 

=  |/,())l  P| .  P)|  *  r'2  [7,0>(  p  2,  p)]  l/2^tn»(p2  p,,p), 

(3.1) 

where  /,0|(  p.  p)  represents  the  optical  intensity  at  the 
source  point  S  and g,n’  (p,  p, ,  p)  is  the  degree  of 
spec 1 1 a)  coherence  of  the  light  at  two  source 
points  S,  and  S2(with  position  vectors  p|  and  p2 
respectively).  Moreover,  for  sources  of  this  class 

. . .  r)  changes  so  slowly  with  the  position  (p) 

across  the  si  mice  that  it  is  essentially  constant  over 
regions  whose  linear  dimensions  are  of  the  order  of 
the  effective  range  of.g,,),(  p2  p(,  p),  i.e.  of  the 
older  of  the  spectial  correlation  length,  lv  say,  of  the 
light  aemss  the  source.  It  is  also  assumed  that  the  lin¬ 
ear  dimensions  of  the  source  are  large  compared  both 
with  /,  and  with  the  wavelength  X  =  c/v. 

Sources  of  this  class  include  the  usual  thermal  (c.g. 
blackhody )  sources  for  which  lv  is  of  the  order  of  the 
wavelength,  hut  other  types  of  sources,  for  which  lv 
may  he  much  greater  than  the  wavelength,  also  belong 
to  this  category.  However,  all  quasi-honiogeneous 
sources  may  be  said  to  W  globally  incoherent,  since 
the  domain  which  they  occupy  is  very  much  larger 
than  their  (spectral)  coherence  area  (=*rr/2). 

To  determine  the  generalized  radiance  ol  the  field 
produced  by  a  quasi-honiogeneous  source  we  proceed 
as  follows.  We  lusl  set  2  =  0  in  eq.  (2.3)  and  (lien  lake 
the  t  ourier  inveise  ol  the  resulting  formula.  T  his  gives 
;m  expression  lor  t')  in  terms  of  the  boundary 
values.  f/,0,(  p  0f  (j{r  ,q  jM  ||)C  plane  r  =  0.  Next 
wc  substitute  this  expression  into  eq.  (2. 10)  and  ob¬ 
tain  the  following  expression  for  the  generalized  radi¬ 
ance  in  the  source  plane. 


238 


P.  s)  =  (k/ 2tt)2,v,  expliXs  ,  ■  p ) 

X  J"  IV(0*(  p,  p'.  i>)  exp(  iA.v,  p’)d2//.  (3.2) 

O 

where  IV(0>(  p.  p',  p)  is  the  cross-spectral  density  of 
the  field  in  the  source  plane.  In  deriving  eq.  (3.2),  eq. 

(2.1 )  was  used. 

For  a  quasi-honiogeneous  source  H,(<b  is  given  by 
eq.  (3.1 ).  and  if  wc  use  that  equation  the  formula 

(3.2)  becomes 

T3  J,0|( p.  .v)  =  (k/2n)~s:  exp( iXr.»1  •  p ) 

X  |/(0)(P.  p)]1/2  f  (/(0)(P\p)1  l,2g(0)(P'  P.p) 

O 

X  cxp(  \ksL  p')  d2p’.  (3.3) 

Since  for  a  quusi-hoinogeneous  source  the  optical  in¬ 
tensity  /(0)(P.  p)  (with  p  fixed)  remains  sensibly  con¬ 
stant  over  regions  whose  linear  dimensions  aie  of  the 
order  of  the  effective  lanpe  lr  ol  g,,,i<  wo  may  replace 
the  factor  p'.  p)|  hv  ( la"{  p,  p)(  ^2  in  eq. 

(3.3)  and  then  take  it  outside  the  integral  sign.  More¬ 
over.  since  the  linear  dimensions  of  a  quasi-ltomogc- 
ncous  source  ate  much  greater  than  /„,  the  integration 
over  o  may  be  taken  over  the  whole  p  '-plane  without 
introducing  an  appreciable  error.  Eq.  (3.3)  then  gives, 
with  very  high  degree  of  accuracy,  the  following  ex¬ 
pression  for  B  J,°\p. s): 

T0j,n)(P ,s)  =  k2szfl°\p,v)g(0>(ksl,v) .  (3.4 ) 

Merej'ttbf/.  js  the  two-dimensional  Fourier  trans¬ 
form  of  g®*(  p',  p),  i.e. 

|<0)(/.  p)  =  (27T)  2  J g,n)(  p'.p)exp(  i/p’)d2p'. 

(3.5) 

The  foimula  (3.4)  shows  that  the  behavior  of  the 
generalized  radiance  of  a  qitasi-homogciieoiis  source 
at  a  point  S  in  (lie  source  plane,  in  a  direction  speci¬ 
fied  by  the  unit  vector  s,  is  determined  by  the  value 
of  the  optical  intensity  at  that  point  and  hv  the  spa¬ 
tial  Fourier  component  labeled  by  the  spatial-frequen¬ 
cy  vector  ksL  of  the  degree  of  spectral  coherence  of 
the  light  in  the  source  plane.  This  result  was  obtained 
previously  by  a  slightly  different  argument  in  ref. 

1 1  l|,eq.(AI0) 
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An  expression  lor  the  general  i/ed  ladiance  of  the 
held  generated  hy  the  qtiasi-honiogcncous  source  at 
any  point  1’  the  hall-space  e  >  rt  is  obtained  at  once 
on  substituting  I  rout  eq  (4.4)  into  eq.  (  2  '))  anil  one 
finds  that 

H,  (r..«)  =  A  •1v..y,l,,(Aji1 ,  i'K'*(r,  s, )  e\p(  tA.v  r),  ( 4.<>) 
where 

( ,\r,s, )  f(>{R.  i  )  /<U)(  P,  i’i  expl  lA.v  p)  d-’p.  (4.7) 


4.  The  asymptotic  limit  A  •  >*  of  the  gcnerali/.erl 
radiance  for  a  field  generated  hy  a  quasi-homogeneous 
source 

l.et  us  now  consider  the  behavior  of  the  expres¬ 
sion  (4.A)  lot  very  short  wavelengths  A  or,  untie  pre¬ 
cisely.  determine  its  asymptotic  limit  as  the  wave 
number  A  2iri\  *  r*=.  I  or  litis  pm  pose  we  carry  out 
the  ihiteientiation  on  the  i  ighl-liand  side  of  eq.  (  2.X) 
and  siitistilule  the  lesullm g  expression  for  the  (iieen's 
I u in  lion  < ,  Mi  1  >  III  eq.  I  t  ’  t  V.  e  tin'll  I  Hid  lll.ll 

f',,u.si » -  r|'  *(/-. -v,  1 1  rJ’V.  v  I.  (4  i) 

w  lieie 


<>,.V  r  k::  />»tp.c)‘ d, 
~n'J.  R- 


^’Vt  I-  -;-r/np.t')"P|i^  p!!  d-p,  (4.4) 

77  *•'  rt  < 


'■i|  R ,  p  )  -  R  *■  x ,  p  ( 4 .4 1 

I  ac!i  ot  die  integrals  in  C(|s.  (4  2 )  and  (4  4)  de- 
ivnds  mi  A  in  two  w  a\  s  via  the  exponential  n  un 
e\p|iA> !>{R.  Pt|  and  via  the  A  dependence  implicit  in 
the  oplic.il  inlcKilv  f>nii  P  .  i>)  /,0,(  p.  Ac/.’rr)  As  A 

hecoiin l.tfct  a.id  t.uee/.  ihe  esponcnlial  Icon  will, 
in  '.’cii.'ial.  oscillate  nioie  and  mote  upidly  .is  the 
point  S  exploies  the  domain  of  in  ten  at  ion.  On  the 
olhei  hand,  loi  auv  fixed  value  ol  A.  (lie  optical  in- 
lensin  i  I  a  tjii.isi-liomoiieiicons  source  v.n ics  slowly 


with  p  :  hence  its  A -dependence  may  he  neglected  in 
the  asymptotic  evaluation  ol  Cj1 '  and  (’J"1.  as  is  clear 
Irom  the  piinciple  ol  stationary  phase  |  l(>|.  Moieovei 
it  is  evident  liom  couipai isoii  ol  the  expressions  mi 
the  i  iidi  l  -iiaiid  sides  ol  cqs.  I  )  2 )  .mil  14  > I  lll.ll  as 
A  *  1  is  ol  liiehei  older  in  I  A  lhauC].'1  Hence 

we  only  need  to  coni  me  om  a  1 1  cm  ion  to  the  asymp¬ 
totic  approximation  to  Cj1  1 

Stiaightloiwaid  application  ol  the  piinciple  ol 
slalionaiy  phase  slums  that,  in  gencial.  ihe  inlegtal 
in  eq.  (4.2)  has  eilhei  one  critical  point  of  the  liist 
kind  or  none  at  all.  I  el  us  set 

r  -  (.v.  r.  :).  r.  ~  (v.  y.  ())  (4.5  ) 

and  let  us  denote  h\  S(l  the  point  specified  try  the  po. 
sit  ion  vectoi 

P  o-r.  {zj\z  ).v, .  (4  (>) 

which  may  readily  he  shown  n  the  plane  :  (I. 

One  finds  that  it  Sn  lies  within  the  sum  domain  n. 
it  is  the  ci it  at  point  ot  the  lusi  kind,  and  that  if Sn 
lies  outside  n.  the  inlei'ial  does  not  have  a  ci  it  ical 
poult  ol  ihe  liisl  kind  We  will  sec  shoilh  lhal  lie 
S(,  lias  a  simple  j'eoineliical  sipml  Main  e 

4V lien  S,,  is  located  within  n.  ihe  asyinploiic  ap 
proximaiimi  to  (\.(r.  ,v)  is  loumJ  n>  he 

C\.(r. )  '  /"V,  (r/y.lv  .iqexpbAi  r) 

as  A  •  117.,) 

When  S()  is  located  outside  o.  the  asymptotic  appio\. 
mtation  comes  from  eontrihutions  ol  critical  points 
of  the  second  kind,  and  isofliighet  order  in  I  A  than 
(lie  expression  on  the  fight-hand  side  of  eq  |4  7;)> 
and  we  may  express  this  tact  ( taking  some  liberty 
w  ith  the  interpretation  of  the  asymptotic  symbol  I  hy 
writing 

f]ir.  s  )  -  ('  as  A  ( 4  "b ) 

On  snhsliiuling  1 1 1 >iii  eqs.  (I  / )  into  eq.  (  4.0)  and 
using  the  tact  that  /,n|  is  /eio  when  ils  ar.Mimem  lies 
outside  ol  o.  we  finally  obtain  the  follow  mi'  asvmp- 
lolie  appioxiinatiou  to  the  genei ah/ed  ladiance  tunc 
lion  ol  a  tield  renei.ited  hv  quasi  In  >■  1 1 >  >n  neons 

souice: 

fi,.(r.  si  -  //,.(/•.  v)  asA  *  "\  ( 4  H ) 

wlieie 
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li,,(r.s)  -  /c 2.v,/<,,,(r|  (z/s,  )s1.i')j?,0,(fcs),  u).  (4,‘>) 

(lie  asymptotic  approximation  (4.8)  (44J)  lo  (lie 
general i/.ei I  radiance  is  t Ire  main  restill  of  this  mile. 

We  will  show  that  il  has  a  niinilier  of  imporlanl  con¬ 
sequences. 

hirst  wc  note  that  according  to  traditional  radioni 
etry  the  rate  at  which  energy  crosses  an  area  element 
d/1  per  unit  solid  angle  around  a  direction  specified  by 
a  real  unit  vector  s  is  given  by  Bv(r,  s)s’n  dA  ,  where  n 
is  the  unit  normal  to  dA.  In  particular  it  follows  from 
this  formula  that  the  rate  at  which  energy  is  radiated 
into  the  far  /.one  per  unit  solid  angle  around  the  s- 
direction  (i.e.  the  radiant  intensity)1  3  across  any  plane 
z  =  z0  =  const.  >  0  is  given  by 

:fijs)=  s,  J  //1,<r,.v)ihf  dj.  (4.10) 

-  -  -  a 

On  sirlislilulini'.  Iron)  ei|.  (4  0)  into  er|.  (4.10)  we 
icadily  find  that 

/>,,(«)  =  ( 2nk  )2sj  /  ""(0,  p)#,0,(  ksl.  p),  (4. 1 1 ) 

where 

7,0)(0.  v)=  (2tr)  2f/,n>(fi.v)d2p.  (4.12) 

If  we  recall  that  ,v.  =  cos  0,  where  0  is  the  angle  that 
the  (real)s-direction  makes  with  the  normal  to  the 
source  plane,  the  right-hand  side  of  the  formula  (4.1 1 ) 
is  found  to  be  precisely  the  expression  for  radiant  in¬ 
tensity  from  a  quasi-homogeneous  source,  calculated 
by  physical  optics  (ref.  ( 1 1 1,  eq.  (4,8)|. 

Il  will  be  convenient  for  the  purpose  of  subsequent 
discussion  to  express  the  formula  (4.0)  in  two  alter¬ 
native  forms,  hirst  we  rewrite  it  as 

W,,('l\s>  r  A2.v./,m(S0.  r>)?,,,,(*s1 .  n)  il  S„  o 

~  0  il  S,, '/  o 

(4.14) 

further  it  follows  from  elementary  geometry  that  the 

'  3  Rigorous  justification  for  the  identification  of  the  expres¬ 
sion  (4  1(1)  with  tile  radiant  intensity  of  physical  optics 
requires  some  additional  considerations,  which  wc  plan  to 
present  in  another  paper 

240 


r 


l  ip.  2.  Illustrating  the  notation  relating  to  the  formulas  (4  14) 
So  is  the  point  in  the  source  plane  whose  position  vector  Po 
is  given  by  eq.  (4.6);  il  is  the  point  of  intersection  with  the 
source  plane  of  the  line  through  the  point  I’  in  the  direction 
of  the  real  unit  vector  j, 

point  Sp,  whose  position  vector  is  given  by  eq.  (4.0). 
is  precisely  the  point  at  which  the  line  through  (’.  in 
the  direction  specified  hy  the  unit  veclot  s  (again  as¬ 
sumed  lo  lie  leal),  iiitciseels  the  source  plane  :  0. 

Hence  eq.  (4. 1 4)  implies  that 

BJ I*,  s)  =  A2s./*0,(S0,  i')£,0,(Asj ,  v)  if  s  6 

=  0  ifj*np 

(4.14) 

where  fip  denotes  the  solid  angle  generated  by  the 
lines  pointing  from  the  source  points  to  P  (see  fig.  2). 

The  first  tluce  terms  on  the  right-hand  side  of  the 
first  line  of  eq.  (4.14)  are  evidently  non-negative.  So 
is  t lie  last  tenn  j?,0,(As, ,  n).  since  it  is  the  Fourier 
transform  of  a  non-negative  definite  function  1 17). 
lienee 

«„<P.s)  "  0.  (4.15) 

I  et  B\ I'1  denote  the  limiting  value  of  when  (lie  spa¬ 
tial  argument  (r  oi  P)  approaches  the  source  plane 
:  -  0.  Since  the  optical  intensity  is  zero  ai  any  point 
Pin  that  plane  which  is  located  outside  the  source 
aiea  o.  wc  have  from  eq.  (4.14) 

«;0,(P..v)M)  ifP£o.  (4.10) 

Finally  wc  see  at  once  from  eq.  (4.14)  that 

B„{V.s)=  H®\Sn.s)  (4.17) 

I  bis  formula  implies  that  BjP.s)  is  constant  along 
each  line  in  the  half-space  ;  >  0. 
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The  tact  that  eq.  (4. 10),  with  Bv  given  ny  eq.  (4.4), 
represents  the  radiant  intensity  as  calculated  from 
physical  optics,  as  well  as  the  results  expressed  by  eqs. 
(4.1 5),  (4.16)  and  (4.17),  show  that  Bv  has  all  the 
main  properties  attributed  to  radiance  We  may, 
therefore,  conclude  that  traditional  radiometry,  with 
the  radiance  given  by  eq.  (4.9),  correctly  describes 
the  behavior  of  fields  generated  by  quasi-homogeneous 
planar  sources  in  free  space,  in  the  asymptotic  limit  as 
k  =  djr/X  -► 

Finally  we  wish  to  remark  that  although  we  derived 
the  expression  (4.4)  lor  the  radiance  from  one  partic¬ 
ular  definition  of  a  generalized  radiance  function  (in¬ 
troduced  in  refs.  |4] ),  we  believe  that  the  same  ex¬ 
pression  will  follow,  in  the  asymptotic  limit  of  large 
wave  number,  from  some  of  the  other  (non-equivalent) 
definitions  of  generalized  radiance  functions,  when 
they  are  specialized  to  fields  generated  hy  quasi-homo¬ 
geneous  sources 

Ack  now  Jedgemen  t 

Wo  wish  to  acknowledge  stimulating  discussions 
with  Prof.  W.T  Weltord,  I  K.S.  and  with  Prof. 

R  Winston  about  the  subject  nuittci  of  this  note. 


1 2  |  V  .1  1  alarskii,  The  e  feel  of  the  turbulent  atmosphere 
i  a  fju  propagatin',  •!  .8.  I)ep  .  Commerce,  Na- 
ttoi.il  ie.|.i,i  al  | a i on  Service.  Spiingfield.  V A, 

I  97  1 1  ,-,t.  ni 

bee  also  j  In, nr.  f  1  i  iUdiofizika  7  fI964)  559 
( 3 1  <il  <  V  i  luruur:  an.  s  I  i  inrskii  Hadinphys.  C7uant. 

I  leclt'.it  I 3  1 19  72/  | lip 

j4  |  lal  V  Wallin  I  )  (i|  I  S-  •....  i  ,il  'i|  ifc22, 

I  bl  V  IVuiilit*  1  Oj  >  S-  \  >i.  i.  r  .  i ’J  lit  )  oil6 

[5  |  I  W  .  il  la,  s  II,"  1 1 1  t  1 1  ■>  in  i  i)f.  1 
K>|  A  i  1  rin.  .j  i  n,  i  \ .o  69  .  i  f  in:. 

j7j  k.M.nlm  'lie, mi*  ,,,*)  r  M  M* ; . -.  .1  i  ipi  Sue.  Am. 

\ I  ( ivm ,  •  so 

1 8  j  A  1  1  nbe.r  i  Ij.ii..,  \. ;  n  ( 1  48 i  i  26 1 
(9|  r  Jannson,  J.  i,|,,  S  *■  V,i.  .  7<  ■  I ‘<Ki<>  1 S44. 

|l(i|  H.M.  Pedersen,  0,.,t. .,  v  t,.  J  m  r"i  ii  877. 

1 1  I  |  M  il  Carter  and  1  Wui.ji,,  i  ■  A.n  67  (1977 ) 

785. 

1121  I-  Wolf.)  OptS-H  ill,  Tj  ,  |9h2)  74  3 
[  I J  |  (  ]  houwkamp,  Rep.  Progi  I’hy*  ilondon:  Ilir  Phv  st¬ 
eal  Sueieiy)  17  (I954|  35 

1 14  |  lord  Kay  leigh.  life  iln.iiy  sound  'reprinted  bv  Dover, 
New  Wk  Vol  It  14s  ,  See  278  (with  a  m<'UilKa1i"n 
aPlnopiMie  1,,  me  ii,ne  ,1*.  pendente  exp(  i2tti  ,  used  in 
llu  pie.eul  papei  i 

HM  I  W  Maul, and  .01.1  I  W.,||  i  I||,|  Soe  Am.  64  (1974  I 
1  273; 

A  Wall  l.e  i  J  Op,  .s,  64(|974)|17S 

[I6|M  II.oii  .,i,.|  I  W.,l(  i'.uniples  of  optics  (Perganmn 
Pfevv, Ijslout  and  New  York,  hlh  oti  ,  1980),  Appendix 

HI 

1 1 7 1  l:  Woll,  optics  (.11  .  .1  .Ml  >87.-4  ( 2 . 1 4  >. 


References 

1 1  |  A  Walther.  J  Opt.  Soe.  Am.  58  (1968)  1256. 


Radiance  functions  that 
depend  nonlinearly  on  the 
cross-spectral  density 


John  T.  Foley  and  M.  Nieto-Vesperinas 


a  reprint  from  Journal  ol  the  Optical  Society  of  America  A 

volume  2,  number  9,  September  1985 


85 


The  US.  Government  Is  authorized  to  reproduce  and  sell  this  report. 
Permission  foi  further  reproduction  hy  others  must  be  obtained  from 
the  copyright  owner. 


Radiance  functions  that  depend  nonlinearly  on  the 

cross-spectral  density 


John  T.  Foley*  and  M.  Nieto-Vesperinas^ 


1 1>  '|Kir!  nu'ii{  nf  I'hvsii  s  unit  Astronomy.  I  'tuvrtsity  nf  Itnrhi-sti-r,  tlnrhcslrr.  JVetr  York  14827 


Received  March  20,  1085;  accepted  May  8.  1085 

Recently  a  new  definition  of  radiance  was  proposed  (•).  Opt.  Soc.  Am.  A  1, 558  ( 1984)]  that  depends  nonlinearly  on 
the  cross  spectral  density  id  the  field  and  satisfies  the  three  major  postulates  of  traditional  radiometry.  We  show 
that  there  are  an  infinite  number  of  such  radiance  functions  Their  utility  is  discussed. 


It  is  well  known  that  one  of  the  main  problems  encountered 
in  the  attempt  to  connect  the  theory  of  partial  coherence  with 
traditional  radiometry1  ( is  that  there  is  no  radiance  function 
that  depends  linearly  on  the  cross-spectral  density  of  the  field 
and  satisfies  the  three  major  postulates  of  traditional  radi 
ometry  for  planar  sources  of  any  state  of  coherence.4-  ’  More 
specifically,  consider  a  planar  source  of  finite  area  I)  located 
in  the  plane  2  =  0  that  emits  light  into  the  half-space  ?  >  0. 
I.et  r  lie  a  two-dimensional  position  vector  in  the  planet  = 
0,  Win.  r-j,  e)  he  the  cross-spectral  density  in  that  plane,  and 
s  lie  a  three-dimensional  unit  vector  whose  z  component  is 
nonnegative.  Kribergr  showed  that  there  is  no  radiance 
function  /fir,  s, e)  that  satisfies  the  following  four  conditions 
for  planar  sources  of  any  state  of  coherence: 

111  /fir,  s,  i>)  depends  linearly  on  W(r,,  r  ,,  cl, 

(111  /fir,  s.  c)  £  o  for  all  r  and  s, 

(III)  /fir,  s,  e)  =  I)  when  r  <  I), 

(IV)  cos  II  ,)  a  /fir.  s,  c)d -V  =  ,/(s.  c), 


r2,  e)<M r)p  cld-r,  =  \,Ai>)<t>n{r2.  c). 


(2) 


The  eigenfunctions  are  orthonormal  over  the  domain  D, 
i.e., 


*(r,  e)c/>„,(r,  eld 


=  /> 


nm » 


Cl) 


and  the  eigenvalues  are  real  and  nonnegative.  Expansion  ( 1 ) 
holds  irrespective  of  whet  her  the  set  of  functions  |</>„(r,  e)j  is 
complete  in  the  Hilbert  space  of  functions  that  are  square 
integrable  over  /-l.*-*' 

In  Ref.  7  the  following  definition  of  radiance  was  pro¬ 
posed: 


H( r,  8,  i') 


cos  II 


expfi/i-s  ■  r'Ki’tr',  r,  c)dV 


(41 

where  (Ur\  r,  r),  the  generating  function,  was  given  by 


where  cos  0  =  s  •  i  and  J(s,  e)  is  the  radiant  intensity  of 
physical  optics.  For  example,  Walther’s  two  definitions1-2  fail 
to  satisfy  requirement  (II)  for  certain  types  of  sources.  '-'’ 

In  an  interesting  recent  paper7  a  new  definition  of  radiance 
was  introduced  lhal  la)  depends  nonlinearly  on  W  (r,.  r-j,  c) 
and  (h)  satisfies  requirements  111)  (IV)  for  sources  o'  any  state 
of  coherence.  The  purpose  < it'  this  ( 'ominimhatinn  is  to. show 
t  hat  I  he  radiance  ( mu  l  ion  ot  Ref  7  is  not  inique  in  t fu  se  re¬ 
spects  By  using  the  methods  of  Ref  7  we  w  ill  show  that  there 
are  an  infinite  number  oi  radiance  ion  -ii.ats  that  sali-.tv 
conditions  (al  and  tl>)  above. 

The  <  ross  spectral  density  ot  lln  ource  ean  be  represented 
by  the  Mercer  expansion’’ 

Wir,.  r  ,.  rl  =  ,\„  i  e  (r ,  rl.y.tr  ..  r).  ID 

where  the  </>,,  Ir.  el  .‘ml  N.li)  r, livi-lv.  tin-  eigen 
functions  ami  the  eigem  aloe-  ol  lie  l-’redhohn  integral 
equation 


/»’(«■',  r,  el  =  x/i(r)L  r,  id  (ft) 

n 

and 


Kqualion  (•))  lories  /fir,  s,  e)  to  lie  noiinegat ive;  therefore 
condition  (II)  is  fulfilled.  Since 

(i(r-.  r.  r)  =  I).  r  ./  I),  (7i 

condition  Mil  is  I  oil  died.  Hv  using  Kqs.  ( I ),  fit),  and  (a)  it 
is  a  straightforward  matter  to  show  that 

W'(ri,  rJt  e)  =  J  (i'*(r ,,  r.  rlfitr-.,.  r.  eld  r,  (8) 

and  il  follows  front  Kqs.  t-l)  and  (8)  that  condition  (IV)  is 
fulfilled  (see  Ret.  7  lor  details). 

We  will  now  show  that  lor  a  given  cross  spectral  density- 
function  there  are  an  inti  nit  e  number  of  generating  functions 
f<t r  .  r  i-i  t bat  obey  Kqs.  1 7 )  and  i  hi.  The  corresponding  ra 
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dtonrc  functions  obtained  by  using  K(|.  (4)  will  tiiorefore  de- 
1 10  in  I  nmilineurly  on  U’t  r  i .  r,  cl.  ulioy  conditions  (III  < IV) 
anti,  in  general,  i«‘  different  Irmn  I  no  railianoo  tiinoimn  nl  Hot. 

Consider  the  expansion 

(ftr-.  *■.  o)  =  \i,{r)  V  ( J'tt/i,. t r f  )</>„, * ( r ,  el,  (9) 

IM  'I 

where1" 


certain  typos  of  sources  Walt  tier's  radiance  functions  ulioy. 
approximately.  I  tie  equation  of  radial  ive  transfer  when  they 
propagate  into  the  half  space  c  >  0. 1  1  -  Since  W ( r  i ,  r,,  rl 

propagates  into  the  halt  space  c  >  u  according  to  two  linear 
partial  differential  equations  (the  Helmholtz  equations  on  the 
variables  r,  and  r  .  respectively),  the  nonlinear  radiance 
functions  may  not  propagate  in  this  simple  manner.  Noser 
theless,  the  nonlinear  radiance  functions  are  interesting  and 
deserve  further  study. 


iIL:K1„<eH-<  (HI) 

m  it 

i  r,  cl  satisfies  Kq  <7l.  By  using  Kqs.  <9)  and  01)  one 
obtains 

(i*(r . ,  r.  cUi'l  r ..  r,  tdd-r 

1 1 

-  HI  V  «m„*ii’iti„„,,ii'i«n,(r|.  ir2,  r).  (11) 

ni  n  ‘ 

Kquation  (1)  can  he  rewritten  as 

VVi r i .  r j,  id  =  v  v  A„ t «■>«>„„  </•„ * t r i ,  tdi'i,,  (r..,e)  (12) 

'i  n ' 

li  lollows  from  Kqs.  (Ill,  112).  and  <41  that  the  generating 
turn  lion  defined  tiy  Kq.  d*>  obeys  Kq.  (Ml  if  and  only  if  the 
expansion  coefficient*  satisfy  the  sealed  unitarilv  condition 

)  ftnifi*(l'|ftmn  (t'l  "  A„(C)A„,|  (Id) 

‘li 

Therefore  tiny  radiant v  junction  of  tin  form  (  I),  where  the 
Lien, ■'■til me  function  is  of  the  form  (ft)  and  the  coefficients 
i  or-,  mi  i/uality  ill))  and  Eq.  iMi,  depends  non  linearly  on 
the  i  r.  ss- s [>cct rnl  density  in  the  source  plain  and  satisfies 
i ■oiidiiions  l// (  i  / V ' )  The  radiance  function  ot  lief.  7  cor¬ 
responds  to  the  choice 

CtmaO’  \  AfTll/')Am,|  1 14) 

\nother  simple  choice  would  tie 

a  ,,,'CI  \  ,\„  it  ,rxp|in...(r)|a.  I  I  fit 

>  e  re  e.e  li  o (el  is  real 

Die  above  result  brings  to  mind  two  questions.  First,  of 
II  tat  | sc.,- 1 hie  nonlinear  radiance  functions  that  arc  possible. 
■'  i  h.-re  one  ‘O.g.,  the  Hermit  inn  one  of  Kef.  7 1  that  is  preter 
.die  '  nless  one  imposes  additional  physical  restrictions  on 
i  he  problem,  t  no  answer  to  this  question  is  clearly  no. 

Second,  are  l  Imse  nonlinear  radiance  fund  ions  preferable 
1  e  the  two  definitions  of  Walt  her"  This  is  an  open  question 
a  tins  turn  .  however  the  following  point*  are  relevant  Kach 
4  Wnlllicr's  radiance  Iiiih  turns  is  nomiegative  when  the 
our,  i  is  quasi  homogeneous11  Since  quasi  homogeneous 
■  nun  i  '  are  globally  spat  ml  lx  incoherent,  VS  all  tier’s  radiance 
•uni  tu ill-  behave  prnperK  tor  the  tvpes  ol  sources  (incoherent ) 
ior  whiih  tr.niilioii.il  radioinetr,  was  developed.  Also,  for 
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